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FOREWORD 


Volume I of this report presents the results of a study 
effort conducted during the 1973-1975 time period. The initial 
study effort had a very broad scope and considered a number of 
different system approaches. As a result of this effort, a 
more detailed study was conducted to define a Shuttle borne 
laser ranging experiment. Volume II of this report presents 
the results of the second study effort. 
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1.0 INTRODUCTION AND SUMMARY 

The purpose of the study was to synthesize and evaluate laser measurement 
systems which showed a capability to meet planned space-based earth physics 
investigations; to define critical system and component performance areas; and 
to develop optimum system approaches based on parametric evaluation of perform- 
ance, cost, and risk. 

The measurement accuracy goals were 2 cm rms range estimation error and 
0.003 cm/ sec rms range rate estimation error, with no more than 1 cm (range) 
static bias error. Further, each selected measurement concept was to be capable 
of meeting these objectives with any of the specified targets. 

The specified links included links between a high altitude, transmit/ 
receive satellite (HATRS), and one of three targets, (1) a low altitude target 
satellite, passive (LATS), (2) an active low altitude target, and (3) a ground- 
based target. The number of links was subsequently expanded to include operations 
with a primary transmit/receive terminal intended to be carried as a Shuttle pay- 
load, in conjunction with the Spacelab program. The primary emphasis of the 
Shuttle links was ranging to passive retro-reflectors on the ground, although 
ranging to other orbital satellites was a desired capability. 

The early part of the study was oriented toward defining the constraints 
and parameters associated with these links. As with most radar links, the 
principle constraints on feasibility were the target signature, the minimum 
allowable beamwidths, and the required received signal energy to achieve the 
desired ranging accuracy. The most critical link is the link between the HATRS 
and the passive LATS, and we compared the various measurement concepts for this 
application. We investigated CW and pulse-burst mode-locked laser measurement 
concepts, to satisfy the range rate measurement goals, These measurement concepts 
were found to be significantly short of meeting either the range or range rate 
measurement accuracy goals for the HATRS to passive LATS link. The major problem 
for the CW mode-locked laser measurement concept was due to the large required 
signal energy to overcome the background from a sunlit earth. A pulse-burst, 
mode-locked laser measurement concept was evolved to reduce the duty cycle, and 
hence reduce the effect of background. This concept was initially considered 
feasible, even though significant development would be required. This concept, 
however, was found to be considerably inferior to the more conventional single 
pulse measurement concepts for range measurements. Since both range and range 
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rate (or range-sum- rate) measurements would be made on the same fundamental 
waveform, the pulse-burst, mode-locked measurement concept provided no fundamental 
measurement advantage and was, therefore, discarded. 

A micro-wave frequency modulated, single frequeney CW laser was considered 
as a remotely possible alternate technique for obtaining range-sum- rate measure- 
ments. However, since the passive LATS considered in this study was based on the 
LAGEOS, the dispersive effect of multiple cube-corners distributed on a spherical 
surface would preclude meaningful measurements at side-band frequencies materially 
greater than achievable with the single pulse, mode-locked laser measurement 
concept, thus this concept was considered inappropriate. 

The remaining two measurement concepts were based on a presently available 
technique (a Q-switched, cavity dumped laser) and a reasonable extrapolation of 
this laser to add a mode-locker within the cavity to reduce the emitted pulse- 
width from the 4 to 6 ns range to values on the order of 0.2 ns. These two 
measurement concepts were retained throughout the balance of the study; no con- 
clusion as to an optimum choice was reached. The Q-switched, cavity dumped 
laser, in ruggedized form, exists today. Although this measurement concept can 
only approach the accuracy goals, development of a space qualified system, using 
this technology, would be relatively straight forward. On the other hand, the 
Q-switched, mode-locked cavity dumped laser can readily meet the measurement 
accuracy goals (with suitable targets), but the laser does not exist, even in 
laboratory form at this time. An oscillator/amplifier combination to achieve 
equivalent performance has been built and is currently in operational evaluation, 
however this system is not readily adapted to a space environment. We concluded 
that the short pulse technique is decidedly superior, and has the long range 
potential to meet the science objectives; however, for more immediate applica- 
tions, such as a Shuttle/Spacelab experiment, the readily available Q-switched, 
cavity dumped laser was a more reliable concept and has the potential to meet 
most, if not all, of the projected experiment goals, 

Detailed evaluation of the target signature of the LAGEOS configuration 
revealed a number of problems which will limit the achievable measurement accuracy. 
The impulse response of the LAGEOS was found to vary significantly with optical 
frequency and velocity aberation, resulting in a varying mean response, which is 
also potentially thermally sensitive. In addition, when the short pulse technique 
is employed, coherent interference effects cause the detected return pulse to vary 
in shape and amplitude on a pulse to pulse basis. The effect of this variation 
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in detected pulse shape and amplitude is dependant on the measurement technique, 
but will certainly increase the rms measurement error. Although the variation 
in pulse shape for the longer pulse laser is less visible, the effect is similar, 
since the centroid of the effective return pulse will vary in a similar manner. 

For targets with smaller effective dimensions, such as planar arrays at a near 
normal incidence, the dispersive effect is virtually eliminated; the pulse-to- 
pulse mean amplitude variations will of course remain. These variations in 
signal amplitude can be countered by averaging or thresholding techniques. 

When the target is active, i.e. , equipped with a laser receiver and transmitter 
the signature problems are virtually eliminated, and viable, high margin links 
can be configured. The selected target terminal configuration employed an 
asynchronous receive/transmit system to simulate an ideal target. The key to 
the operation of the system is precise recording of the time-of-arrival of the 
received pulses and of the time-of- departure of the transmitted pulses. Sub- 
sequent processing can then extract the propagation time estimates, even though 
the clocks at the HATRS and LATS (or ground) terminals are not perfectly 
synchronized. 

A major concern during the study was the accuracy with which the transmit 
beam and receiver fields-of-view could be oriented. For the link between the 
HATRS and the passive LATS, very narrow transmit beamwidths and small receiver 
fields-of-view will be required. The net pointing accuracy requi red will be on 
the order of one microradian. We determined that this accuracy is theoretically 
achievable, based on current and projected state-of-the-art in guidance and 
navigation technology, although many problems remain to be solved. 

The pointing accuracy required for the Shuttle experiment is considerably 
less stringent, more on the order of 50 to 100 microradian. However, the naviga- 
tion uncertainties and attitude reference errors in the Shuttle systems results 
in pointing uncertainties on the order of milliradians to degrees unless auxiliary 
subsystems are added to improve the references. The acquisition and pointing 
control problem for a minimum cost experiment requires further study to ascertain 
feasibility and interface requirements. 

The net conclusion is that both the Shuttle experiment and the long range 
earth physics mission objectives can be met with existing and projected technology; 
no extraordinary development problems are envisioned. The concepts are con- 
sidered fundamentally sound, and reflect the application of current ground-based 
laser radar technology to space-borne system applications. 
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2.0 LINK ANALYSIS 


2.1 SUMMARY. One of the early activities in the study was to select a link 
analysis summary format which would serve, with minor modification, to document 
both active and passive link margin summaries. The format shown in Table 1 
evolved from this activity, and is similar to an optical communication type of 
link margin analysis format. For links with passive targets, all entries shown 
are used. When analyzing a link with an active target, lines 7, 8, and 9 are 
deleted. 

The format shown in Table 1, for ranging with passive targets, is identical 
in content to the classical radar range equation, with a few (4tt/a ) factors 
added to convert from range to free space loss and from radar scattering cross- 
section to target receive/transmit antenna gain product. The classical radar 
range equation can be written,^^^ 

Pr = PjGjAeCrS.^S^/(4u)2 (1) 


= received signal power (watts) 

Pj = transmitted signal power (watts) 

Gy = transmit antenna gain (relative to an isotropic antenna) 
Ag = receiver effective antenna area (m ) 
a = radar scattering cross-section of the target (m ) 

Sy,S^ = transmittance of the transmit and receive paths 
R = range to the target (m). 


The target radar scattering cross-section is defined as the product of the 
projected target area, the transmit antenna gain of the target, and the target 
reflectivity. 


^^^Skolnik, M. I., "Introduction to Radar Systems," McGraw-Hill, New York (1962), 
pg 37-5. 
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TABLE 1 


HATPf; TO LA TP (PASS I VE) 

Cl PC HIT MARGIN PUMMARY 

OPTICAL WAVELFNCTH .53 HIGRra^iETERP, RANGE 30000. KM 
?!!L PEv^TDTH = .2 0 N A NOP EC O MR S 



PARAMETER 

U N I TP 

VALUE 



! 

TPANPriTTFR POWER 

DBW 

- 20.00 

10.0 

MILL I WAT TP 

P 

TRANPriTTER LOPPEP 

D3 

-.97 



3 

TRANPMIT ANTENNA GAIN 

D3 

121.07 

5. 

M IC RORAD I A NP 

4 

POINT IMG LOPP 

DB 

- 1-.39 

I. 

M ICRORADIANP 

5 

FREE PPACE LOPP 

D13 

-297.04 



fS 

ATMOPPHERIC LOPP C DOWN) 

DB 

.00 



7 

TARGET GAIN PRODUCT 

D3 

208.50 




A TM 0 P P H F R I C 1 ,OPP (UP) 

DB 

.00 



Q 

FREE PPACE LOPP 

DB 

-297; 04 



1 0 

RECEiyE ANTENNA GAIN 

D3 

131.00 

.PI 

METERP DIA. 

1 1 

RFCEIVFP LOPPEP 

DB 

- 3 . 98 



12 

RECEIVED PIGNAL POWER 

DBW 

- 15 9.85 



13 

TRAMPflT PULPF RATE 

D3/PEG 

10.00 

1 0 . 

pulpep/pec. 

14 

ENERGY PEP PUL PE 

• DBJ 

CO 

• 

1 



15 

ENERGY PER PHOTON 

DBJ 

-iS4.2C 



1 5 

RECEI VED PHOTOMP/PULPE 

DB 

14.4 1 



1 7 

QUA NT DM EEFICIEMCY 

Dil 

- C .02 

25.0 

PERCENT 

10 

RECEIVED P-E/PULPE 

DB 

S.39 

P.91 

P-E 


19 

BACKOROIIND RADIANCE 

DBW 

-17.70 

.0170 

WATTS/Mi2-A-P 

20 

PEC El VER FOV( PTERRAD . ) 

1)3 

-87.07 

5 0. 

M ICRORADIANP 

21 

OPTICAL FILTER BANDWIDTH 

D3 

P.99 



22 

RFC FIVE ANTENNA AREA 

• DBM2 

-5 .5 1 



23 

RECEIVER LOPPEP 

DB 

-3.98 




RFCEIVED BACK. POWER 

D3W 

-107. 2P 



C'S 

ENERGY PER PHOTON 

DBJ 

- 184.2P 



2P 

R EC E I VED PHOT 0 NP / SEC 

DB/SFG 

77.00 



27 

QUAMTHM EFFICIENCY 

D3 

-P.02 



2.0 

BACKGROUND P -E/PEG 

D3/PEC 

70.98 



2.9 

RECEIVER GATE WIDTH 

DH-PEC 

-8P.99 

2. GO 

NANnPECnMDP 

3 0 

BACKOROUND P-E/GATE 

DB 

- IP. 01 

. 03 

P-E 



LI NK f ARGIM 




18 

RECEIVED P-E/PULPE 

DB 

8.39 

P.91 P-E 

31 

PEQUIRED P-E/PULSE 

D3 

.00 

1.00 P-E 

32 

LI Ml- l ARGi .M 

DB 

8.39 

P.91 
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a = r G-J.J 


( 2 ) 


r = target reflectivity (assumed = 1) 

Gyj = target transmit gain = 4ir/fi 
Aj = target projected area (m ) 

= effective target transmit beamwidth (steradians) 


The radar range equation can then be rearranged to the form shown in Equation (3), 

2 


Pr ^ ^I'^T'^R^T^r'^TT^TR^FS 


Gr = Ag = radar receive antenna gain 


(3) 


TR 


YJ Aj = target receive antenna gain 


Lps = (x/4irR)'^ = free space loss 
X = wavelength (m) 


These terms are shown in the first twelve lines of Table 1. Provisions were 
made to divide the loss terms into significant components, and the Gyj'G^I^ product 
was shown rather than the individual terms. The parameter values are shown in dB, 
thus line 12 is the sum of the first eleven lines. For easy referenee, the major 
parameters are also printed in natural units. 

Line 1 shows the transmit power, and is chosen to provide adequate signal 
margin within the constraints of feasibility for space-qualifiable hardware. 

The second line shows the optical transmittance of the transmit optics path. 

This term is a function of the optical design of the transmitter and the optical 
elements necessary to form and steer the transmit beam. These losses are 
discussed in Section 2.2. 

The radar transmit antenna gain, relative to an isotropic (4-?r steradian) 
radiator, is shown in line 3, and is discussed in Section 2.3, as is the transmit 
antenna pointing loss term, shown on line 4. 
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The free space loss terms are computed for the range and wavelength shown 
in the heading, and appear (for passive target ranging) in lines 5 and 9, 
Atmospheric losses are also shown in two places, lines 6 and 8. These losses 
are discussed in Section 2.4. 

The target gain product is shown in line 7, and discussed in Appendix F. 

The radar receive antenna gain, line 10, is discussed in Section 2.3., and the 
receiver optical losses, line 11, are discussed in Section 2,2, 

The received signal power, line 12, in dBW, is simply the sum of the first 
eleven lines. Most of the systems considered in this study employed pulses of 
optical energy for ranging. Line 13 shows the pulse repetition rate, and line 
12 minus line 13 defines the energy per pulse, shown in line 14. 

( 2 ) 

The energy per photon, is defined as hv, where h is Plank s constant 
-34 

(6,6256 X 10 Js) and v is the optical frequency (c/x). This is shown in 
line 15 (also line 25), and line 14 minus line 15 yields the mean number of 
received signal photons per pulse, line 16, delivered to the photodetector. 

Line 17 shows the detector quantum efficiency, and the sum of lines 16 and 17 
yields the mean number of detected signal photoelectrons per pulse, shown 
on line 18. 

Many of the potential laser ranging concepts are severely background 
limited, and all are affected to some extent. Lines 19 through 30 show the 

background level calculations. The background source, for most of the links, 

is the sunlit Earth. The background calculation is similar to the signal 
calculation. The background radiance, in watts/m -angstrom-steradian, is 
shown in line 19. This is documented in Section 2,4. The receiver field-of- 
view, steradians, is shown in line 20. Note that the field-of-view in ste- 
radians is approximately tt/ 4 times the square of the full plane FOV (in radians). 
This approximation is virtually exact (within 0.001 dB) for fields-of-view less 
than 0.1 radians. 

An optical filter is used to reduce background il Tumi nation, and the filter 
bandwidth, in Angstroms, is shown in line 21. The receiver antenna area (m ) 
is shown in line 22, and the receiver optical losses in line 23 (same as line 

11). The received background power, line 24, is simply the sum of lines 19 

through 23. 

^^Ubramowitz and Stegun, "Handbook of Mathematical Functions," Dover Publications, 
New York, 5th Edition, pg 7. 
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Conversion from background power to photoelectrons per pulse is accom- 
plished by converting first to photons/sec (line 26), then, considering quantum 
efficiency (line 27), to detected photoelectrons per second (line 28). The 
detection techniques for the candidate laser measurement concepts all employ 
some form of time gating of the detector output. The effective receiver gate 
width (line 29), summed with the background rate yields the detected background 
photoelectrons per gate (or pulse), shown in line 30, 

Following line 30, the final margin calculation is shown. Line 18 is 
repeated, for convenience, followed by the required photoelectrons per pulse, 
line 31, and the signal margin in line 32. The required signal photoelectrons 
per pulse are discussed in Appendix F, and are either input as a specified 
parameter or calculated based on the detection technique and observed back- 
ground level. 

Appendix D presents the link margin analysis summaries for the various 
links, the computer program used to perform the summary analyses, and the 
specific parameters used. 

2.2 OPTICAL EFFICIENCIES. The transmittance of the receiving and transmitting 
optics is a function of the number and types of optical elements employed, 

and the coatings used. A conservative estimate of the number and type of elements 
for a typical optics system was made, and that figure was used for all terminals. 
Table 2 summarizes the element counts, the efficiency per device, and the total 
efficiency computation for the typical terminal receive and transmit optics. 

The most significant loss element is the narrow band optical filter. In the 

o 

range of optical pass-bands from 1 to 100 A, the filter efficiency is not 
grossly a function of bandwidth, or wavelength for that matter, but rather a 
strong function of the materials, techniques, and the supplier. The available 
data shows peak transmissibilities ranging from '^40% to as high as 85% in the 

G 

vicinity of 5A bandwidths in the visible and near IR frequencies. Sixty percent 
transmittance is used as a reasonable compromise value. 

2.3 ANTENNA GAIN AND POINTING LOSS CALCULATIONS 

2.3.1 HATRS . The transmit antenna gains for the candidate HATRS laser radar 
systems are required to be relatively large, to permit operation with passive 
targets at long range. The receive antenna gains, similarly, are required to 
be relatively large. The receive antenna size varies from 0.6m to 1.0m dia- 
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TABLE 2 

TYPICAL RADAR TERMINAL OPTICAL EFFICIENCIES 


INTERFACE 

TRANSMIT 

RECEIVE 

Reflecting Surfaces 

9 @ 0.99 

11 @ 0.99 

T ransmi ssi ve El ements 

5 @ 0.98 

5 @ 0.98 

Bandpass Filter 

NA 

0.6 

Beamsplitter 

2 @ 0.98 

2 0 0.98 

TOTAL 

0.793 

0.466 

Allocated 

0.80(-0.98 dB) 

0.40 (-3.98 dB) 


meter, and is assumed to be a centrally obscured (Cassegrain) configuration. 

The receive antenna gain is simply, 

Gp = (D^ - d|) (4) 

Dp = primary diameter 

Dg = secondary diameter 

The transmit antenna gain is a function of the beam intensity profile, 
dimensions, and truncation. The nominal transniit beam has a Gaussian intensity 
distribution. We have assumed an offset feed arrangement is used, as shown in 
Figure 1, in order to use the same primary optics as the receiver. The transmit 
beam diameter is truncated by the beam combining mirror, and the beam diameter 
shown in the primary aperture is the truncation diameter expanded to final 
size. The far-field gain of a truncated Gaussian beam (on-axis) is given by 
Equation (5). 


2(7T 


V^)‘ 




= truncation diameter (m) 

2 

= l/e power diameter of the beam (m) 


( 5 ) 
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OPTICS ARRANGEMENT 



FIGURE 1 


The truncation at the beam combining mirror also results in a loss of 
receive antenna gain. The significance of the losses to the transmit and 
receive gain is dependant on the desired transmit beamwidth (and gain). If 
the image of the transmit beam truncation aperture does not overlap the 
secondary or fall outside of the primary, then the product of the receive and 
transmit antenna gains is given in Equation (6). 

2 2 ^ 

The two Toss terms, on the right, can reduce the net signal significantly 
unless considerable care is used to choose the parameters. Since the transmit 
beam is truncated, the far field antenna pattern is no longer exactly 
Gaussian. For the systems considered here, the approximation is relatively 
good. This can be seen by computing the far field pattern of the truncated 
beam.^^^ 

^^^Buck, A. L., "The Radiation Pattern of a Truncated Gaussian Aperture Distribu- 
tion," Proceedings of the IEEE, March 1967, pg 448. 
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( 7 ) 


e|K 6} = J ^ ° Jo(Ker) rdr 

0 

Jq(x) = Bessel function of the first kind 
K = 2tt/X 

0 = angle off bores ight, radians 
_2 

r^ = e power beam radius 
r-j. = truncation radius 


Clearly, 




( 8 ) 


Equation (7) can be reduced to a Bessel function series as shown in Appendix G. 
The result is, 

n 

, . 1 0 ^ ' 

E 


:|K6) 4 t 0 „(Ker,) /'■t 

n=1 \Ke kJ " ' 




(9) 


Thus , 


G(e 


GTo 


E 

n=l 


(A) ^ 

177 


2 




(10) 


Equation (10) is thus the gain loss for the truncated Gaussian beam as 
a function of angle off boresight. Combining (10) and (5) results in 
Equation (11). 


G^(0) = 2(uDq/a)' 


oo 


E 

(^) ''"‘“'■t’J 



2 2 
2 -2r^/r^ 


( 11 ) 


By expanding the Bessel functions in series form, and rearranging terms, 
the expression in Equation (12) is obtained, and was used for the computation of 
the radiation pattern plotted in Figure 2. 
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where 


ej(x) 


^ -^ = the incomplete exponential 
K=0 


Figure 2 shows that the effect of truncation is to decrease the on axis 
gain, and slightly broaden the far field pattern. The pattern calculated for 
the truncated Gaussian beam can be approximated quite well for pointing errors 
to about 2-8 micro radians as a nontruncated Gaussian beam with an on-axis gain 
equal to the gain of the truncated Gaussian beam. 


FAR FIELD ANTENNA PATTERN 



POINTING ERROR >^RAD FIGURE 2 
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2.4 MISCELLANEOUS LINK PARAMETERS 


2.4.1 Background Noise and Sources . Most of the links considered in this study 
operate predominantly with the Earth in the field-of-view of the receiver. The 

background radiance of the sunlit Earth, considering only diffuse refi action, 

9 ° 2 ° 

is '^^0.017 watts/m -A-steradian at 0.53 yra wavelength, and 0.00875 watts/m -A- 

steradian at 1.06 ym wavelength. 

An active LATS receiver would normally view deep space, with a reasonably 
large field-of-view. The moon, under these conditions, is the brightest body 
in the sky other than the sun, and has a spectral radiance of 4,7 x 10 ^ 

W/m -A-steradian at 0.53 ym and 1.6 x 10" WjAn -A-steradian at 1.06 ym . 

An active ground-based terminal would see atmospheric scattered sunlight 
during daylight operation, which has a spectral radiance^^^ of about 3.5 x 10 ^ 

p o -3 2^ 

w/m -A-steradian at 0.53 ym and 2-5 x lO" W/m -A-steradian at 1.06 ym. 

These values have been used in the link margin analyses. In all cases, 
operation with the sun (or strong specular solar reflections) in the field-of- 
view has not been considered either feasible or necessary. 

2.4.2 Atmospheric Transmissibility - For links involving ground targets, two 
loss mechanisms have been considered, nominal atmospheric attenuation due to 
molecular and aerosol absorption and scattering, and degradations due to scintil- 
lation effects. The nominal atmospheric transmission model is shown in Figure 3 
for 0.53 ym and 1.06 ym wavelengths, as a function of ground terminal altitude 

for zenith angles of 20°, 40° and 60°. The data is for a clear atmosphere, and 
f5l 

Elterman's^ ' model was used. At modest zenith angles, an atmospherie trans- 
missibility (each way) of '\>60-70% is reasonable at 0.53 ym and 80%-90% at 1.06ym. 

For link analysis purposes, we have arbitrarily allocated a 3 dB loss term 
to account for atmospheric scintillation. The basis for this allocation is an 
assumption that each ground target will be composed of a number of retro- 
reflectors. We further assumed that the retro- reflector spacing and arrangement 
would be such that each retro-reflector effectively contributing to the return 
signal energy could be considered as scintillating independently. Consequently, 

( 4 ) 

' ^Pratt, W. K. , "Laser Communication Systems," John Wiley & Sons, New York (1969), 
pg 121-125. 

^^^Elterman, L. , "UV, Visible, and IR Attenuation for Altitudes to 50 km, 1968," 
AFCRL, ERP NO. 285 (1968). 
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even under severe scintillation conditions, the worst case fade depth should be 
relatively modest if the number of retro-reflectors contributing to the return 
is large. 


ATMOSPHERIC TRANSMISSION VS GROUND 
TERMINAL ALTITUDE FOR ZENITH ANGLE 



1^| I I I ^ i 

2000 4000 6000 8000 10,000 12,000 



GROUND TERMINAL ALTITUDE, hj - km 


FIGURE 3 
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2.5 TARGET SIGNATURE AND GAIN PRODUCT. The passive targets considered in this 
study are all augmented with one or more high quality, optical retro-reflectors , 
also referred to as cube-corners or corner reflectors. In order to determine 
the amount of signal energy returned to the radar, it is necessary to determine 
the energy returned from each reflector on the target, and to combine the 
reflections from each reflector at the radar. The return from each reflector 

is a function of the angle of incidence of the incoming wave, and the magnitude 
and direction of the velocity aberration vector, also referred to as point-ahead 
angle. The individual reflector gains were computed as discussed in Appendix F. 

2.6 SIGNAL POWER REQUIREMENTS. The required signal photoelectrons (pe) per 
pulse, to achieve the desired range estimation accuracy goals, is a function of 
the laser pulse shape, the target characteristics, and the range estimation 
process employed. The candidate range estimation processes may be categorized 
as estimation processes which are optimum, in some sense, and those which are 
sub-optimum but offer implementation advantages. These range estimation processes 
depend on sensing the time when a laser pulse is emitted and when the reflected 
pulse is detected. The time of arrival (TOA) of the reflected pulse can be 
estimated in many ways. 

A TOA estimation process, which is optimum in a maximum likelihood sense, 
is described in Appendix E. Also discussed and analyzed in Appendix E are 
several other estimation processes which offer implementation advantages, and 
approach the performance of the maximum likelihood estimator. These techniques 
and their analytical treatments are appropriate for shot-noise limited applica- ' 
tions where the target is a single cube-corner retro- reflector in a nondispersivd 
channel . 

However, when the target is composed of multiple retro-reflectors, which 
is characteristic of targets which have been synthesized for the HATRS links, 
scintillation is encountered, due to coherent interference, at the detectors, 
resulting from summing the returns from each retro-reflector. This type of 
scintillation is discussed in Appendix F, and the mean detected signal amplitude is 
shown to have an exponential distribution in the limiting case. 

The remaining system functions which are covered in the link margin analysis 
program result from the postulated use of a ground based beacon, of some sort, to 
aid in the initial acquisition (spacially) of a set of ground-based retro- 
reflectors clustered in a small geographical area on the ground. Two types of 
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beacon systems were hypothesized to span the spectrum of feasible types for 
active ground beacons. 

First, a relatively high PRF, low energy per pulse system, such as might be 
obtained from a GaAs laser diode beacon, was postulated. For this technique, the 
angle tracking loop bandwidth and the beacon PRF determine the required beacon 
energy in pe/pulse. 

The other type of beacon considered was a low PRF, modest energy per pulse 
beacon, similar to existing ruby lasers. In this case, the spatial acquisition 
and tracking function is better described as a bang-bang system, where the primary 
considerations are false alarm and missed detection probabilities. The computa- 
tion algorithm which is employed for this purpose is based on estimating the 
threshold setting, and subsequently the required signal, to satisfy the Chernov 
bounds.' ^ This technique is modestly conservative, and was chosen to evade 
numerical difficulties for large signal and background level conditions. These 
bounds are determined, using transform techniques, as shown in Equations (12) 
through (17). 


The second moment generating function is defined as. 


^^(s) = In 




sx 


dx 


For a Poisson process, with mean intensity m. 


i'x(s) “ In 


^ nL g“m+sk 
k=o 


= m(e^ - T) 

The Chernov bounds are then, for false alarms. 


PfA = > y} i exp|if.^(s) - si|j][(s)| 

provided s > o 


'y " 


( 12 ) 


(13) 


(14) 


(6) Van Trees, H.Lo, "Detection, Estimation, and Modulation, Vol . I", 
Wiley and Sons, New York (1967), pages 121-122 
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and, for missed detection, 

Pun ' PP)* ‘ Ti iexp(t|> (s) - s*l(s)> , 

“ ' ’ * ^ '' ' 'v = <(s) (15) 

provided s < o 

In these expressions, y is the threshold setting in photo-electrons. 

Then, with m = nj^T, equation (14) bounds the false alarm probability. Thus, 

PpA 1 exp{Y - m - Y ln(y/m)| , m = n^^T , (16) 

Y > m 

The first step, therefore is to find the threshold setting, y, which just 
satisfies this criterion. Then, with m = n^ + nj^T, we iteratively solve to 
find the value of m, hence which just satisfies Equation (17), derived from 
Equation (15). 

Pmd iexp|Y - m - Y ln(Y/m)j- , 

Y < m 

These equations and algorithms have been implemented in the link margin analysis 
program presented in Appendix D. 


2-14 



3.0 LASER RANGING SYSTEM CONFIGURATIONS 

The laser ranging system configurations developed during the study were 
selected to maximize the achievable ranging accuracy and to be compatible with a 
variety of targets. Two major active terminals were defined. First, a terminal 
suitable for use in a high altitude, transmit/ receive satellite (HATRS), was 
defined to operate with low altitude target satellites (LATS), either active or 
passive, and ground based retro-reflectors. The other terminal was configured 
for use as an experiment carried on a Spacelab mission, and is primarily intended 
for operation with ground-based retro-reflectors. Target configurations were 
defined for both the active and the passive targets, either orbital or ground 
based. 

System performance is limited by the type of laser selected. From a ranging 
accuracy standpoint, the most attractive laser is a mode- locked, cavity dumped 
Nd:YAG laser, which outputs very narrow (<0.2 ns) pulses of modest energy. How- 
ever, this laser does not exist in space qualified form at this time, although 
all of the (known) necessary elements have been developed for different applica- 
tions. Thus, it was deemed advisable to have an alternate system configuration 
based on the current state-of-the-art, a Q-switched, cavity dumped NdrYAG laser, 
with a characteristic pulse-width of 4 to 6 ns. The choice between the two laser 
types is a matter of risk assessment. In the narrow pulse case, the risk is 
principally that unknown factors may delay the development of the new laser con- 
figuration. The risk, for the more conventional laser techniques is on develop- 
ing a space qualified laser with sufficient amplitude and pulse shape stability 
to enable meeting an over-all rms ranging accuracy equivalent to l/30th to 
l/50th of the pulse width, which is a nontrival goal under the best of circum- 
stances. There are, of course, other components which present development risks. 
These components, however, are common to both concepts, and limit, if at all, 
only the time scale for implementation. 

3.1 HATRS TERMINAL. The HATRS terminal is configured to operate with a variety 
of targets at very long range. The key terminal parameters postulated for this 
study are summarized in Table 3 for both the short pulse and the long pulse laser 
configurations. 
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TABLE 3 

HATRS TERMINAL PARAMETERS 


PARAMETER 

SHORT PULSE 
CONFIGURATION 

LONG PULSE 
CONFIGURATION 

Wavelength 

0. 53 ym 

0.53 ym 

Pul sewidth 

0.2 ns 

4.0 ns 

Transmit Energy 

l.OmJ/pulse 

10.0 mJ/ pulse 

PRF 

10 pps 

10 pps 

Average Transmit Power 

10.0 mW 

100.0 mW 

Transmit Beamwidth 

5 . 0 y rad 

5.0 yrad 

Telescope Primary Diameter 

0.61m 

1.0m 

Detector Quantum Efficiency 

25% 

25% 

Detector Field of View 

50. yrad 

50. yrad 

Optical Filter Bandwidth 

A 

5A 


The physical configuration which was adopted is inustrated in Figure 4. 

The interface with the host spacecraft may be either hard or gimbals, depending 
on the host spacecraft design and functional requirements. The basic assumption 
is provision of a precision optical bench, which supports all alignment critical 
components. The alignment critical components include the transmit and receive 
optics, the laser and the detectors, and the attitude reference sensors. The 
attitude reference system is a steller-inertial system which employs strapdown 
precision rate-integrating gyros, strapdown star sensors, and a radar telescope 
focal plane mounted star tracker to establish and maintain a precision attitude 
reference of the radar receiver boresight vector. 

A typical ranging operation begins with establishing an inertial attitude 
reference update, accomplished by slewing the receiver telescope to one or more 
known stars. The star position with respect to the receiver boresight is 
established, using either an offset star sensor located on the telescope focal 
plane or the radar system quadrant PMT detectors with the narrow bandwidth optical 
filters temporarily removed. Once the attitude reference is established, the 
telescope can be open-loop pointed at the expected target position. The radar 
system employs a narrow transmit beam and a relatively wide receiver field-of- 
view, thus the acquisition process is principally a matter of scanning the 
transmit beam until the target is illuminated. 
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HATRS TERMINAL PHYSICAL CONFIGURATION 



FIGURE 4 

The nominal region of pointing uncertainty includes components due to un- 
certainty in both the HATRS and the target ephemeris data, attitude reference 
errors, and transmit/ receive alignment errors. Since the transmit and receive 
boresight vectors must be offset to account for velocity aberration (point-ahead 
angle) effects, there is also a transmit beam deflection drift and scale factor 
error included in the pointing uncertainty. 

Once the target is illuminated, the receiver detects the angle of arrival 
of the reflected signal, and slews to zero the mean receiver pointing error. The 
alignment error between the transmit beam and the receiver boresight is detected 
and corrected by dithering the transmit beam. The system maintains pointing 
control using the inertial reference system; the corrections are filtered and 
used to eliminate bias errors. Range measurements are obtained in a similar 
manner. For spatial acquisition, the range gate is widened to account for range 
prediction uncertainties. Once the target is reliably detected, a split-gate 
range tracker is used to extract a range prediction error estimate which is 
filtered and used to bias the range prediction to zero the mean range tracking 
error. 

One implementation approach is shown summarized in Figure 5 . A central 
computer is used to perform the necessary data processing and prediction functions. 
A typical range measurement begins with selection of a nominal transmit pulse 
departure time. This departure time is chosen to ensure that the receiver will 
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not be blanked by the TR shutter at a time when a return pulse is expected. The 
desired departure time is delivered to the pulse generation electronics unit 
(PGE). The PGE compares the commanded departure time with the central clock 
time, and at a programmed time prior to coincidence, the PGE initiates the pulse 
generation sequence. The sequence begins by energizing the lamp pump. Subse- 
quently, the cavity is Q-switched, and lasing begins. For the mode- locked laser 
technique, the mode-locking RF signal is also applied to the intracavity mode- 
locking element. After the laser has reached stable conditions, the cavity is 
dumped and the high energy pulse proceeds through the optics system. Prior to 
dumping the cavity, the receiver is temporarily blanked by a switchable attenua- 
tor or shutter, and a sizable bias voltage is applied to the quadrant PMTs 
photocathodes to gate them off. 

The laser output pulse passes through a frequency doubler. Then a small 
amount of the optical pulse energy is extracted and delivered to the pulse 
departure time detector. The output from this detector is delivered to a 
threshold detector, which strobes the departure time estimator at threshold 
crossing. The departure time detector output is also integrated and it's magni- 
tude is converted to digital format for subsequent processing in the central 
computer to correct the departure time estimate for actual pulse amplitude. The 
optical pulse proceeds through the transmit beam steering mirror assembly, which 
establishes the desired transmit offset angles. The pulse then passes through 
a hole in the beam combining mirror, and exits the telescope assembly in an off- 
set position, as shown in Figure 5. Subsequently, the switchable attenuator in 
the receiver is opened. 

The next pulse arrival time, estimated from the vehicle ephemeris data, 
the pulse departure time, and the filtered range tracking error, is delivered to 
the receiver gate timer. At a programmed time prior to coincidence with the 
central clock time, the receive gate timer issues a PMT enable command, and 
subsequently issues gate sample commands as shown schematically in the figure. 
Four sampling modules, one for each PMT, are used to measure the intensity of 
the received signal at three sampling times, the early gate, the center gate, 
and the late gate sample times. The sampling modules employ a sliding window 
integrator technique to measure the PMT output during the three overlapping time 
periods. The sampling module outputs are then converted to digital data, and ' 
delivered to the central computer for subsequent processing. 
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The system can be most readily converted to operate with active targets 
by assuming two color operation, i.e., the HATRS radar transmits at 0.53 ym, and 
the target transmits at 1.06 ym. This would require inserting a dichroic beam 
splitter in the receive optics, and adding a detector, probably a silicon photo- 
diode, or semiconductor photocathode PMT. The active target receiver signal 
processing will differ significantly frcm the passive target receiver signal 
processing techniques. The operational concept is based on recording pulse 
departure times and pulse arrival times at both ends of the link, with no attempt 
made to synchronize or provide a transponder function in the target. This free 
running operation provides all of the data necessary for subsequent estimation 
of range, assuming the data from both terminals is available. Both the HATRS 
and LATS transmit beamwidths and receiver-fields of view are selected to allow 
totally open loop pointing in this mode. 

The selected configuration for the HATRS active LATS receiver is shown in 
Figure 5 and uses a PMT detector with a weighted counter approximation of a 
maximum likelihood time-of-arrival estimator. If a silicon photodiode detector 
were chosen, instead of the PMT, the time-of-arrlval estimator module would 
be modified slightly to replace the weighted counter with an analog matched 
filter, since the signal would be preamplifier noise limited. In either case, 
the AGC amplifier maximum gain is clamped to minimize false alarms during high 
background level conditions. 

3.2 SPACELAB LASER RADAR EXPERIMENT. The Spacelab radar terminal configuration 
employs many of the concepts evolved for the HATRS terminal, but is significantly 
different in several important areas. Table 4 shows the key terminal parameters 
assumed for each of the two candidate laser techniques. 

One major difference between the HATRS and Spacelab terminal is inclusion of a 
two-color ranging technique in the Spacelab terminal to permit assessment of atmo- 
spheric delay by simultaneous ranging at 0.53 ym and 1.06 ym wavelength. A second 
significant difference is that the Spacelab terminal does not track the returned 
signal angle of arrival. The pointing technique is principally open-loop, with 
provisions for a ground beacon angle of arrival measureraent to initialize the point- 
ing system at the start of a ranging pass over a cluster of ground-based retro- 
reflectors. Third, the range estimation technique is based on a pulse-by-pulse 
detection process , rather than by tracking the average signal over a number of 
pulses. Finally, less stringent pointing accuracy requirements and Shuttle attitude 
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TABLE 4 

SPACE LAB TERMINAL PARAMETERS 


PARAMETER 

SHORT PULSE 
CONFIGURATION 

LONG PULSE 
CONFIGURATION 

Wavelength 

0.53 urn 

1 . 06 pm 

0.53 ym 

1.06 pm 

Pulsewidth 

0.2 ns 

0.2 ns 

4. 0 ns 

4.0 ns 

Transmit Energy 

1.0 md/pulse 

4.0 mJ /pulse 

10.0 mJ/pulse 

40.0 mJ/pulse 

PRF 

1 pps 

1 pps 

1 pps 

1 pps 

Transmi t Power 

1.0 mW 

4.0 mW 

10 mW 

40 mW 

Transmi t Beamwi dth 

375 yrad 

750 prad 

375 prad 

750 prad 

Telescope Primary Diameter 

0.3m 

0.3m 

0.46m 

0.46m 

Detector Quantum Efficiency 

25% 

1.0% 

25% 

1.0% 

Detector Field of View 

375 yrad 

750 prad 

375 prad 

750 prad 

Optical Filter Bandwidth 

5A 

5A 

5A 

0 

5A 


motion disturbanGes are satisfied with a gimballed Cassegrain telescope in a 
Coude^ configuration. Alternatively, a fixed telescope with a gimballed flat 
mirror could be employed, subject to installation limitations. 

Figure 7 is a functional block diagram of the selected Spacelab laser con- 
figuration, and shows representative implementation techniques for these func- 
tions. The attitude reference sensors are considerably less critical than for 
the HAIRS terminal , however they must still be tightly coupled to the laser 
radar optical reference bench. 

Since the system functions on a pulse-by-pulse basis, transfer from one 
ground target to another can be accomplished readily between pulses, particularly 
if the angular separation is within the dynamic range of the tracking benders. 

The term, benders, is used generically, since it appears that torsion-bar-mounted 
torque motors may be more rugged and reliable than piezo-electric bender bimorphs 
for beam steering functions. 

Since the Spacelab is in a low altitude orbit, for a relatively short mission 
duration, and is subject to potentially sizable ephemeris changing perturbations, 
difficulties in converting ranging data into data of geodedic significance are 
anticipated. One possible solution is to define a Spacelab laser radar experiment 
which provides "simultaneous" ranging nteasurements with two or more ground 
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targets. The most obvious technique is to increase the transmit beamwidth and 
receiver field -of -view, and provide some sort of dual return pulse detection 
technique. However, this approach appears to result in serious link margin problems 
even for targets relatively closely clustered. A more attractive technique from 
energy management considerations is to divide the transmit pulse into two equal 
parts, and form two individually targeted transmit beams. Individual beam steering 
for sizable angular offsets ('v +1®) requires a relatively good linearity device 
0.2%). The most attractive device for this application appears to be preeision 
Risley prisms with 10 bit optical encoders for deflection and orientation feedback 
control . 

The selected system concept can, however, roughly approximate this performance, 
by jumping from target- to- target between pulses, resulting in a tightly inter- 
laced set of data. It was beyond the scope of this study to determine the impact 
of these competing concepts on the value of the measurements, thus these alterna- 
tives are simply presented for consideration. 

3.3 ACTIVE LATS TERMINAL CONFIGURATION. The active LATS terminal is very similar, 

functionally, to the HATRS terminal, and most of the same support elements are re- 
quired, although the performance requirements can be relaxed to some extent. 

Based on our evaluations, the receive aperature weight and cost are not a significant 
part of the subsystem weight and cost when the diameter is less than several inches 
A one-inch diameter (0.0254m) aperature was chosen as a reasonable value for the 
receiver. For the narrow pulse HATRS radar, with the beam spoiled to TOO urad, 
this choice resulted in a 31.6 dB link margin, based on requiring '\-20 pe/pulse to 
ensure detection of every pulse. The rms time of arrival estimation error, for 20 
pe/pulse is ^^32 ps. 

For the transmitter, the required product of energy per pulse and transmit 
antenna gain, to provide a 6 dB margin above the 20 pe/pulse level, is '^/23 dBJ. 

If the transmitter outputs 0.1 mJ/pulse, the transmit beamwidth can be -^4 mrad. 

We allowed about 0.5 mrad for pointing error, which is almost totally allocatable 
to attitude reference errors. This pointing accuracy (0.03®, la) is well within 
the capabilities of relatively inexpensive stellar-inertial attitude reference 
systems, but considerably more stringent than normally achievable for simpler 
attitude reference techniques. 

Figure 8 summarizes the selected terminal configuration. A fixed telescope 
with a gimballed flat was chosen for this terminal, to minimize complexity. An 
alternate design, with three additional folding mirrors, can be configured to 
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FIGURE 8 FUNCTIONAL BLOCK DIAGRAM - ACTIVE LATS TERMINAL 
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provide full hemispheric coverage for a minimal weight penalty. 

An alternate functional concept is shown in Figure 9 > which also shows the 
special HAIRS receiver configuration to extract one-way Doppler frequency measure- 
ments. If a precision clock, with suitable frequency accuracy and stability, is 
provided, at both ends of the link, one-way Doppler frequency measurements can 
provide reasonably accurate range sum rate measurements. In this concept, 
the active LAIS terminal is composed of a gimballed array of retro-reflectors, 
for ranging purposes, and a modest beamwidth (1.06 vm wavelength), 200 Mpps mode- 
locked laser, whose mode-lock frequency is phase-locked to a standard frequency 
derived from a precision Rubidium clock. 



FUNCTIONAL BLOCK DIAGRAM - ALTERNATE ACTIVE LATS TERMINAL 
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4.0 LINK DYNAMIC ANALYSES 

Two of the ranging links evaluated in the Laser Ranging Systems Study 
involve a High Altitude Transmit- Receive Satellite (HATRS) in a 22,000 km circular 
orbit ranging from ground based retro-reflectors or retro-reflectors on a Low 
Altitude Target Satellite (LATS) in a 1000 km circular orbit. These passive 
ranging links are severely power limited due to the extremely long ranges, and 
hence very narrow laser transmit beamwidths (5 microradians) must be employed. 

The laser pulse rate for the selected HATRS terminal configuration is 10 pulses/ 
second. The resulting laser telescope pointing requirements are very stringent 
for these long, passive links because of the narrow transmit beamwidth. Further- 
more, the high accuracy pointing during target acquisition and tracking must be 
achieved in open-loop fashion to some extent since no wideband pointing error 
information is available with a 10 pps pulse rate and the low energy in the return 
signal collected at the HATRS. 

At the beginning of the acquisition phase, laser telescope pointing angle 
commands are computed on the basis of HATRS attitude and position estimates and 
estimated target location. The pointing angle commands then are transformed 
into spacecraft attitude or gimbal position commands, depending upon the system 
mechanization, and the laser telescope is slewed to provide nominal alignment 
of the telescope boresight axis along the HATRS-to-target line of sight (LOS). 

The pointing error of the laser telescope boresight axis relative to the true 
HATRS-to-target LOS at initiation of acquisition is a function of HATRS and 
target ephemeris errors, HATRS attitude reference error, and misalignment between 
the HATRS attitude reference and laser ranging systems. 

If the expected pointing error is larger than some specified fraction of 
a transmit beamwidth, a sequential search of the initial pointing uncertainty 
area must be performed until a return signal is received by detectors in the 
HATRS laser terminal. The acquisition strategy that is devised will be 
strongly influenced by short-term stability of the pointing error since it 
defines the overlap required between adjoining subareas during the sequential 
scan process and hence it affects acquisition time. Thus the open-loop pointing 
requirement for acquisition in the long, passive links is a function of the 
allowable acquisition time and the acquisition strategy employed. However, a 
preliminary requirement of 10 yrad or better is felt to be a reasonable perfor- 
mance objective. 
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Although wideband pointing error information is not available in the long, 
passive links, it is conceptually feasible to determine static and very low- 
frequency (i.e., less than 0.02 Hz) boresight errors since the reflected signal 
strength from a passive target is directly related to HAIRS pointing error at the 
instant the laser pulse was transmitted. Two concepts for measuring the static 
and low-frequency boresight errors consist of (1) a quadrant array of ranging 
detectors for angle of arrival estimation augmented by a conical scan of the 
transmit beam to detect receiver/transmitter boresight errors and (2) a single 
ranging detector which employs a slow conical scan for transmit beam pointing 
error estimation from return signal amplitude measurements. Therefore, after 
the target is acquired, one of these techniques could be used to estimate the 
static and low-frequency boresight error and appropriate compensation for these 
errors could be provided during tracking. The primary pointing error contributors 
which are calibrated with these techniques are ephemeris errors and static 
misalignment between the attitude reference and the laser ranging system. 

The most stringent pointing requirement for the long, passive links occurs 
during the tracking phase when the target must be tracked in an open-loop 
manner within 1.0 yrad pointing error. The achievable tracking accuracy depends 
upon spacecraft and environmental disturbances, attitude reference stability 
(i.e., short-term drift and noise), and HATRS-to-target LOS kinematics. 

In summary, extremely accurate HATRS and target location prediction, pre- 
cise HATRS attitude determination, and precision laser telescope slewing are 
required for acquisition and tracking in the long, passive links. The 10 yrad 
or better open-loop pointing requirement for acquisition is treated as a some- 
what "soft" requireraent since target acquisition can be achieved with larger 
pointing errors. The primary impact of larger pointing error at the beginning 
of acquisition is longer acquisition time. The 1 yrad open-loop pointing 
requirement for tracking is a firm one for a 5 yrad transmit beamwidth. 

4.1 EPHEMERIS ERROR EFFECTS. The positions of the HATRS and the target must 
be known very accurately in order to define the HATRS-to-target line-of-sight 
vector with an accuracy compatible with open-loop pointing within five or ten 
microradians. 

4.1.1 HATRS and LATS Ephemeris Prediction . The LATS and possibly the HATRS 
currently are envisioned to be drag-compensated spacecraft. An accelerometer 
is used to measure the effect of all nongravitational forces acting on this 
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type spacecraft such as the forces due to solar radiation or meterorite impact. 
Then a "compensation" system is exercised to cancel the sensed nongravitational 
forces. The important result is that orbit prediction errors with drag- 
compensated spacecraft are thought to result primarily from ground station track- 
ing errors and errors in the geopotential field model used for orbit prediction. 

GSFC performed an in-house study to assess ephemeris prediction capability 
for the HATRS and a drag-compensated satellite in a 300 KM circular orbit. 

Their computer program estimates 100 parameters including position and velocity 
components for both satellites, a tracking bias error for each of six tracking 
stations, a drag residual for the Tow altitude spacecraft, and 81 geopotential 
terms. Their conclusion is that for a drag-compensated satellite that has been 
tracked by six stations (2 passes per day) for ten days, the uncertainty 
associated with predicting the position of the satellite 24 hours ahead is 
approximately 3 meters (la) provided all terms are solved for. This total 
error can be broken down into a 2.5 meter cross track component, a one meter 
along track component and a two centimeter radial component. These projections 
serve as the basis for estimating effects of ephemeris errors on achievable 
open-loop pointing system accuracy. 

4.1.2 Pointing Angle Error Sensitivities. Sensitivities of HATRS pointing angle 
error to HATRS and LATS position errors were derived for coplanar HAIRS and LATS 
orbits (1) with the line-of-sight tangent to the low altitude orbit and (2) with 
the low altitude satellite directly under the HATRS in order to bound the effects 
of ephemeris errors. Similarly, pointing angle error sensitivities for the HATRS- 
to-ground links were derived with the ground site (1) on the HATRS horizon and 
(2) directly under the satellite. Analytical expressions for the pointing angle 
error sensitivities and numerical sensitivities for the HATRS-to-LATS link are 
presented in Table 5. The numerical sensitivities to ephemeris errors for the 
HATRS-to-ground link are somewhat smaller since the range is larger. 

The numerical sensitivity values for the HATRS-to-LATS link show that the 
pointing angle error per unit position error is approximately 0.05 prad/meter 
if acquisition is performed with the LATS directly below the HATRS. Assuming 
that the along track and cross track positions errors for both satellites are 
the same, the root-sum-square (RSS) HATRS pointing error is approximately O.TO 
prad/meter. The numerical sensitivities of pointing error to HATRS and LATS along 
track and cross track position errors are about 33 percent lower if acquisition 
is performed when the HATRS-to-LATS line of sight is tangent to the LATS orbit. 



TABLE 5 

HATRS POINTING ANGLE ERROR SENSITIVITIES TO HATRS AND LATS POSITION ERRORS 

(COPLANAR ORBITS) 



C|^ are along track, radial, and cross track errors for HATRS. 
Cj^ are along track, radial, and cross track errors for LATS. 

R are HATRS radius, LATS radius, and HATRS-to-LATS range. 








The effects of radial position errors are negligible in an RSS sense since the 
radial errors are much smaller than the along track and cross track errors, and 
the pointing angle error sensitivities also are much lower. The total RSS 
pointing error for acquisition with the LATS on the "horizon" is approximately 
0.072 yrad/meter assuming equal errors. 

4.2 ATTITUDE REFERENCE SYSTEMS. The sub-are second open-loop pointing require- 
ment for the long, passive links dictates stellar sensing, and hence only stellar- 
inertial attitude reference systems (ARS) are considered. High-quality reference 
gyros provide relatively high-bandwidth, short-term data and star sensors provide 
precise low-bandwidth attitude measurements at discrete times. Optimal filtering 
generally is employed to process these measurements in order to estimate and 
compensate for important measurement errors and to estimate spacecraft attitude. 
Spacecraft angular rates can be derived from incremental attitude data or can be 
measured with another set of gyros if higher bandwidth rate data is required. 

Gyro drift rate bias must be estimated and appropriate compensation must be 
provided in the attitude estimation algorithms to achieve high-accuracy attitude 
determination. Sub-arc second attitude determination might also require periodic 
calibration of additional attitude reference sensor parameters such as gyro scale 
factor and gyro input axis alignment. Parameters which require infrequent 
calibration, say every few days, could be estimated using ground-based computer 
algorithms to minimize onboard computer requirements. 

Gyro performance should not be a limiting factor in achieving sub-arc 
second attitude determination in the 1980's. Projected performance of the third 
generation gyros (TGG) is compatible with sub-arc second pointing, and current 
high-quality gyros such as the Honeywell GG334 are close to being adequate for 
some attitude reference system mechanizations. A comparison of current gyro 
performance and performance goals of the Draper Lab TGG is shown in Table 6. 

The Draper Laboratory TGG has undergone significant laboratory “Resting, and 
results of these tests demonstrate the feasibility of meeting the specified 
performance goals. The frequency response of the Draper Lab TGG is essentially 
flat out to 30 Hz which is significantly better than the frequency response of 
current high-quality reference gyros. 

Gyro drift rate bias is a very large error contributor unless appro- 
priately compensated. Random drift also is an important error source in any 
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TABLE 6 

GYRO PERFORMANCE DATA 
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*TOTAL DRIFT UNCERTAINTY 
♦♦UNSPECIFIED 






high-accuracy stellar-inertial reference system. The Honeywell GG334 gyro 
has a 0.01 deg/hr (la) random drift specification, but much better performance 
has been observed during laboratory testing. The drift uncertainty performance 
goal for the Draper Lab TGG is less than 0.004 sec/hour. This assumes that 
the absolute drift of the gyro can be calibrated and compensated with resolu- 
tion that is cormensurate with the 0.004 sec/hour specification. 

Other important gyro performance parameters are scale factor stability and 
input axis alignment stability. Scale factor stability, which is a measure of 
the angular uncertainty associated with torquing the gyro through a specified 
angle, is extremely important if the spacecraft is slewed through large angles 
or if input rates to the gyro are high. A scale factor stability specification 
of 100 to 200 parts per million (PPM) represents current state-of-the-art gyro 
technology. However, the scale factor stability can be calibrated to 10 PPM 
or so. The performance objective of the TGG is quoted as 1 PPM with 0.1 PPM 
potentially achievable. 

The current state-of-the-art in strapdown stellar-inertial attitude reference 
systems is exemplified by the Honeywell Space Precision Attitude Reference System 
(SPARS) to which 3 to 4 s^, la, 3-axis performance capability is attributed. 

The SPARS employs two strapdown starsensors which are solid state devices 
utilizing a silicon photodetector array as the sensing element, and three strapdown 
pulse-rebalanced, rate integrating gyros (Honeywell GG334 series). Primary SPARS 
attitude error contributors are star sensor random error (2 sec, la) and random 
gyro drift (0.01 deg/hr, la). No full -up SPARS has flown to date, but a modified 
version with relaxed accuracy requirements (15 sec, 3a per axis) is scheduled to 
be flown early in 1975. 

The SPARS initially was designed for low altitude orbit missions, but Honey- 
well has performed a study for the Air Force in which they evaluated SPARS perfor- 
mance in high altitude orbits— up to synchronous altitude. The study results 
indicate that a SPARS type of attitude reference system warrants consideration 
for the HATRS application. Higher accuracy star mappers and higher quality gyros 
are required for the laser ranging system than for the SPARS, but the necessary 
components could be available by the early 1980's if projected improvements in 
sensor performance are realized. 

The primary impact of high altitude orbits (hence lower orbit rates) on SPARS 
performance is increased acquisition time since it takes much longer for filter 
convergence with the less frequent observable star transits due to slower star 
field scan. The frequency of observable star transits also defines the attitude 


4-7 



reference update interval during normal operation, but the longer update intervals 
should not be a significant problem by the 1980's with the extremely low drift 
gyros currently being developed. 

The state-of-the-art in stellar-inertial attitude reference systems using 
gimballed star trackers is typified by the TRW Precision Attitude Determination 
System (PADS) which was designed to provide a 3.6 £ec, lo, per axis attitude 
reference in an equatorial, synchronous orbit. The PADS design activity is an 
outgrowth of the Precision Pointing Control System (PPCS) study performed by 
TRW for Goddard Space Flight Center since early 1970. The PADS employs a two- 
axis gimballed star sensor using a photomultiplier tube detector and precision 
strapdown, pulse- rebalanced, rate integrating gyros. Telescope gimbal angle is 
measured by an Inductosyn. Kalman filtering is used to incorporate the star 
sensor measurements and to provide periodic estimates of spacecraft attitude 
and gyro drift rate bias. 

Higher accuracy stellar sensing is achievable on spacecraft such as the 
HATRS which have large optical systems if the star sensor is boresighted to 
the experiment optics. Integration of the star sensor and the large optical 
systems results in larger aperture diameter and focal length than normally is 
feasible within size and weight constraints of separate star sensor assemblies. 

4.3 OPEN-LOOP POINTING SYSTEMS. Two generic approaches can be employed for 
pointing a space-based payload toward a desired target. One method of implement- 
ing the pointing requirement is to mount the payload on a set of gimbals and to 
slew the gimballed payload assembly in order to achieve and maintain alignment 
of the payload sensor boresight axis with the host spacecraft-to- target line- 
of-sight. The second approach, body pointing, is to rigidly attach the payload 
to the spacecraft structure and to perform the pointing function using the 
spacecraft attitude control system. 

Brief descriptions and typical operational scenarios for three open-loop 
pointing system mechanization concepts are discussed in this section. Two of 
the pointing systems use the gimballed payload approach and the third employs 
body pointing. The objective of this section is to identify major system com- 
ponents and primary pointing error contributors in order to assess the 
feasibility of achieving open-loop pointing capability between 5 to 10 yrad 
by the early 1980's in the long passive links. All three mechanization concepts 
require significant improvements in component performance relative to current 
state-of-the-art in order to realize this open-loop pointing accuracy goal. 
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However, only components which currently are being developed or components 
which are reasonable extropolations from current state-of-the-art are accept- 
able since full implementation of the laser ranging system is planned for the 
early 1980' s. 

4.3.1 Gimbaled Telescope Pointing Systems . Numerous gimbaled telescope point- 
ing system implementation concepts can be contrived, with the system differences 
being primarily the type and/or location of the attitude reference sensors. Two 
implementation concepts which characterize the basic options available with the 
gimbaled pointing system approach are discussed in this section. 

One of the gimbaled telescope pointing systems evaluated in this study 
employs a SPARS type of ARS which consists of strapdown star sensors and strap- 
down reference gyros, all mounted on a precision optical assembly (POA) which 
is a precision mounting structure. The POA provides a stable thermomechanical 
interface to the spacecraft for establishing and maintaining precision alignment 
between pointing system components. Another star sensor is located on the focal 
plane of the laser telescope and is used solely for alignment calibration. The 
second gimbaled system employs a star sensor located on the laser telescope 
focal plane as the primary attitude reference sensor and strapdown gyros mounted 
inside the laser telescope gimbals. The telescope gimbals are rigidly attached 
to the POA in both systems. 

Gimbaled System with SPARS Type of ARS (System I) . The basic elements of 
this open-loop pointing system are a SPARS type of attitude reference system and 
a gimbaled laser telescope assembly. The ARS consists of two or more star 
mappers and high-quality strapdown reference gyros, probably third generation 
gyros. The star mappers are wide field-of-view sensors with no moving parts; 
they depend upon spacecraft motion for star field scan. They contain no servo 
loops or angle encoders, and hence these star sensors are mechanically simpler 
than gimbaled star trackers. 

The gimbaled telescope assembly contains the laser telescope, gimbals and 
associated drive electronics, and gimbal angle encoders. Also, an image dissector 
type of star sensor is located on the laser telescope focal plane for use in 
calibrating alignment between the laser telescope and the attitude reference 
system. 
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The HAIRS is controlled about a local vertical /orbit plane reference. Con- 
tinuous star field scan is provided by the nominal spacecraft pitch rate that is 
required to maintain a local vertical orientation, and thus attitude reference 
updates can be provided throughout acquisition and tracking. 

Whenever the laser ranging system is to be exercised, laser telescope point- 
ing angles are computed using spacecraft attitude and position estimates and 
estimated target location, These desired pointing angles are transformed into 
laser telescope gimbal position commands, and nominal target acquisition is 
achieved by slewing the telescope through the required gimbal angles. The initial 
pointing uncertainty area then is searched to acquire the target, if necessary, 
and the static and low-frequency boresight errors are estimated. Then the track- 
ing phase is initiated, and the gimbals are driven in accordance with measured or 
derived spacecraft angular rates during tracking. 

The achievable pointing accuracy with this gimbaled system probably is 
limited by gimbal related errors. The two types of gimbal related errors which 
are of primary concern are gimbal position reference error due to angle encoder 
inaccuracy and gimbal control loop tracking error due to gimbal servo nonlineari- 
ties such as friction and torque motor hysteresis. Sub-arc second attitude 
determination should not be a limiting factor by the I980's if projected improve- 
ments in sensor performance, particularly gyros, are realized. Gyro drift rate 
uncertainty is the most important gyro error contributor with this approach. 

Scale factor stability and input axis alignment stability are less important gyro 
errors since the spacecraft angular rates are fairly small with the spacecraft 
closely controlled about a local vertical/orbit plane reference. The largest 
spacecraft angular rate is orbit rate which is essentially constant in the 
nominally circular HAIRS orbit. 

Only the short-period rates, which are small, produce attitude reference error 
due to scale factor instability and gyro input axis misalignments. The gyro 
measurement error due to scale factor and input axis alignment stability effects 
with a constant input rate is compensated for in the state estimate of gyro drift 
rate since this measurement error is a bias error in the short term. 

The current state-of-the-art in gimbal angle encoder performance is inadequate 
for the Laser Ranging System application with system mechanizations which rely 
upon these devices for gimbal position determination during target acquisition 
and tracking. No angle encoders with readout accuracy specifications in the sub- 
arc second region have been space qualified to date. Both Inductosyns and optical 
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shaft encoders have fractional sub-arc second measurement potential, but these 
high-accuracy angle encoders probably will be fairly large with diameters of 
twelve inches or greater since increased resolution is generally obtained by 
increasing encoder dimensions. An angle encoder with readout accuracy of 0.1 
prad or better is required with this gimbaled system in order to provide 1 prad 
pointing during tracking, and this is well beyond the performance of any angle 
encoder currently being developed, to our knowledge. 

Farrand Controls, Inc. has a twelve inch diameter Inductosyn with a >0.5 s^ 
accuracy specification if the discs are directly mounted to the shaft. The 
accuracy specification denotes the precision to which angular position can be 
measured with respect to an absolute external reference. This should not be 
confused with the resolution specification which defines the smallest movement 
which can be measured or reproduced by the angle encoder. Thus the accuracy 
specification includes mechanical misalignments and electronic noise and biases 
in addition to the resolution of the device. 

Baldwin Electronics has a ten inch optical shaft encoder with 21 bit resolu- 
tion (0.62 sec) and a full circle accuracy quoted as one bit RMS plus resolution. 
Baldwin currently is working on a 16-inch diameter optical encoder with a 24 bit 
resolution but they have not determined its readout accuracy capability so far. 
Their goal is to calibrate the 24 bit encoder v/ithin jjD.l ^ full circle accuracy. 

ARS Sensors Inside Laser Telescope Gimbals (System II). One approach to 
improving open-loop pointing accuracy with a gimbaled laser telescope system 
is to locate the attitude reference sensors (i.e. star sensor and gyros) inside 
the laser telescope gimbals. The star sensor is located on the focal plane of 
the laser telescope and hence misalignment between the attitude reference axis 
and the laser telescope boresight axis is minimized. Also, integration of the 
star sensor with laser telescope optics results in larger aperture diameter and 
focal length than normally is feasible within size and weight constraints imposed 
upon a separate star sensor assembly, and hence utilization of laser telescope 
optics provides greater stellar sensing accuracy potential. The two primary 
advantages of placing the gyros on the gimbals is a significant relaxation in 
the gimbal angle encoder readout accuracy requirement and the availability of 
higher bandwidth gimbal motion data relative to System I. 

The acquisition scenario with this gimbaled pointing system concept differs 
substantially from the acquisition scenario for the gimbaled pointing system 
with the attitude reference sensors strapped to the precision optical assembly. 
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When a command is received to perform a ranging experiment, the laser teleseope 
assembly is slewed to acquire a specified guide star. The inertial reference 
is zeroed and a second guide star is acquired. Large angular separation of the 
two guide stars (ideally 90 degrees) is desired in order to derive accurate 
three-axis attitude information, but smaller angular separation is acceptable in 
order to maintain reasonable gimbal angle range requirements. The gimbal angle 
range required during acquisition depends primarily upon the number of star 
measurements which can be obtained within the allowable acquisition period and 
star measurement noise. Finally, after an attitude reference has been established, 
the laser telescope pointing angles are computed and the laser telescope is 
slewed through the commanded gimbal angles. Then the initial pointing uncertainty 
area is searched, the target is acquired, static and low-frequency boresight 
errors are estimated, and tracking is initiated. 

The gimbal angle encoder provides gimbal position information during acquisi- 
tion of the first guide star. The gyros, which are mounted on the telescope 
gimbals, are employed to provide a reference for precision slewing during sub- 
sequent guide star acquisition and target acquisition maneuvers. Therefore scale 
factor stability and input axis stability are very important error parameters 
with the high gyro input rates and the large angles through which the gyros must 
be torqued during the slewing maneuvers. Current gyros are inadequate but 
acceptable position measurement accuracy is achievable over fairly large angles 
using third generation gyros if their quoted performance goals are met. For 
example, the telescope boresight alignment error following a 30 degree maneuver 
is approximately 0.1 sec with a scale factor stability of 1 PPM. 

Effects of gimbal dynamics such as friction probably are the largest tracking 
error contributors. However, with the gyros mounted within the gimbals, higher 
bandwidth gimbal jitter measurement is available than with the gyros strapped 
to the POA, since, for example, the frequency response of the Draper Lab TGG is 
essentially flat out to 30 Hz. 

A disadvantage of boresighting the star sensor to laser telescope optics is 
that no star measurements can be obtained during tracking when the laser telescope 
is Earth oriented. C(xnplete reliance is placed upon the gyros for maintaining 
an attitude reference of sufficient accuracy. However, attitude error propaga- 
tion due to gyro drift uncertainty during tracking should not be a problem with 
the third generation gyros currently being developed. Attitude error due to a 
0.10 sec/hour drift rate is a tolerable 0.05 sec for a 30 minute tracking period 
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assuming a linear relationship between attitude error buildup and time, and 
this is an order of magnitude larger drift uncertainty than is projected for 
the third generation gyros. 

Comparative Evaluation of Glmbaled Pointing Systems . A qualitative compar- 
ison of the two gimbaled telescope pointing systems is presented in Table 7. 

There are many implementation options available within the class of gimbaled 
pointing systems, but the two implementation concepts compared in Table 7 ade- 
quately serve as a basis for discussing the important tradeoffs which should be 
considered in selecting a gimbaled pointing system for the HATRS. 

The basic conclusion is that the attitude reference sensors should be 
located inside the payload gimbals in order to minimize the effects of gimbal 
angle encoder error and gimbal servo nonlinearities on pointing system perfor- 
mance. This conclusion is based upon projected gimbal angle encoder readout 
accuracy attainable with Inductosyns or optical encoders of reasonable size. 

If the ARS sensors are mounted within the gimbals, gimbal position information 
is required only during acquisition of the initial guide star, and hence a 
moderate accuracy angle encoder (say 18 to 20 bit full circle accuracy) is 
sufficient. Also, low frequency gimbal dynamics can be measured if the gyros 
are strapped to the gimbals. 

The gimbaled pointing system with the ARS sensors mounted within the 
gimbals has other comparative advantages in addition to making achievable pointing 
accuracy essentially independent of gimbal angle encoder errors. Boresighting 
the primary star sensor to the laser telescope optics eliminates the need for 
a separate alignment calibration sensor, provides higher stellar sensing accuracy 
potential, and minimizes the misalignment between the laser telescope and star 
sensor boresight axes. However, star measurements cannot be performed while 
tracking ground based or near Earth targets, and large gimbal angle ranges are 
required for guide star acquisition. An additional star sensor mounted orthog- 
onally to the primary sensor could be used to decrease gimbal range requirements 
and to provide periodic attitude reference updates during tracking, if required. 
This additional star sensor could be a lower accuracy device than the one bore- 
sighted with the laser telescope optics. 

A comparative disadvantage of gyro placement inside the payload gimbals is 
increased sensitivity of laser telescope pointing error to gyro scale factor 
stability and input axis alignment stability errors due to higher gyro input 
rates and the larger angles through which the gyros are torqued during guide 
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TABLE 7 

GIMBALED TELESCOPE POINTING SYSTEM COMPARISON 
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star acquisition maneuvers. However, the pointing error attributable to these 
gyro erros is compatible with sub-arc second pointing if projected TGG performance 
goals are met. Another adverse effect of mounting the gyros inside the gimbals 
is increased power dissipation in the gimbaled assembly and attendant thermal 
control problems. 

Some advantages of the system in which the ARS sensors are strapped to the 
precision optical assembly (System I) are less stringent gyro requirements, 
"continuous" attitude reference update availability, and a smaller gimbal angle 
range requirement (== +25 degrees). However, it is felt that these advantages do 
not outweigh the primary comparative disadvantage of the system which is the 
requirement for an extremely accurate angle encoder. 

4.3.2 Body Pointing System . Body pointing is an attractive alternative to 
gimbaled telescope pointing for the Laser Ranging System because it would 
eliminate the potential problems directly attributable to the gimbals. A brief 
description of a candidate body pointing system and a typical operational scenario 
are presented in this section. 

A SPARS type of stellar-inertial attitude reference system is employed with 
two or more strapdown star sensors and strapdown, rate integrating, high quality 
reference gyros (third generation gyros) mounted on the precision optical assembly. 
The star sensors would be oriented to enable periodic attitude updates to be pro- 
vided during laser ranging experimentation. The time between attitude reference 
updates depends primarily upon spacecraft angular rates and detector sensitivity. 
The laser telescope also is rigidly attached to the POA, and another star sensor, 
probably an image dissector type, is located on the laser telescope focal plane. 
This star sensor is used for initial on-orbit alignment and periodic calibration 
to compensate for alignment shifts due to thermal gradients. 

A typical scenario for the acquisition and tracking phases with body pointing 
is as follows. The target acquisition phase begins with the HATRS in some orienta- 
tion which enables a precise attitude reference to be maintained between ranging 
experiment periods. The desired inertial orientation of the laser telescope is 
determined on the basis of HATRS and target position estimates. Then the space- 
eraft attitude maneuvers required to point the laser telescope along the HATRS- 
to-target LOS are determined using the estimated inertial attitude of the HATRS, 
and the HATRS is slewed accordingly. The resulting orientation of the laser 
telescope boresight axis relative to the true HATRS-to-target LOS will be within 
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specified error limits which allows acquisition to be completed within a prescribed 
time interval using a sequential search technique if necessary. After the target 
is acquired, a conical scan boresighting technique is employed to estimate static 
and low-frequency boresight error. Then the tracking phase is initiated, and the 
spacecraft precisely tracks the HATRS-to-target LOS. 

Precision spacecraft slewing throughout the acquisition and tracking phases 
is performed by the HAIRS attitude control system using reaction wheels, control 
moment gyros, or some hybrid combination of the two control actuators. Attitude 
control accuracy requirements are tighter with body pointing than with gimbaled 
telescope pointing since spacecraft angular rates are a direct measure of laser 
telescope pointing stability with the body pointing approach. The factors affect- 
ing achievable pointing stability with a momentum exchange attitude control sys- 
tem include induced mechanical vibrations due to rotor imbalance and spin bearing 
imperfections, servo nonlinearities such as friction torques and CMG torquers 
ripple and compliancei and sensor noise. All these problems have been studied 
in depth within the industry for other space-based optical system applications 
(e.g.. Large Space Telescope). 

4.3.3 Body Pointing versus Gimbaled Telescope Pointing . A comparative evalua- 
tion of the body pointing system and the gimbaled telescope pointing system 
with the attitude reference sensors located inside the payload gimbals is pre- 
sented in Table 8. Body pointing has the basic advantage of being the simpler 
to mechanize. There are no gimbals, and hence no gimbal servos and no gimbal 
angle encoders. However, more precise attitude control is required with body 
pointing than with the gimbaled systems since the gimbal control system can 
compensate for observable, spacecraft motion. 

Both pointing systems require high accuracy stellar sensing. Star sensors 
located on the focal plane of the laser telescope provide higher accuracy stellar 
sensing potential than is achievable with separate star sensor assemblies of 
reasonable weight and size, but the use of laser telescope optics Drecludes 
attitude reference updates when tracking ground based targets. 

When body pointing is used, it is possible to locate star sensors on the pre- 
cision optical assembly oriented to allow star measurements while the telescope 
is earth oriented (actiye ranging). The boresight star tracker is used only to 
periodically calibrate between the normal mode star sensors and the telescope 
boresight. For the gimbal led system, however, it is not at all certain that 
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TABLE 8 

BODY POINTING VERSUS GIMBALED TELESCOPE POINTING 
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star sensors can be located for star measurements while the telescope is earth 
oriented (host vehicle interference), thus we elected to rely on the boresight 
star sensor for periodic ARS updates between ranging operations. 

A very "quiet" spacecraft is required for both the body pointing and gim- 
baled telescope pointing systems. This is an inherent requirement with body 
pointing since laser telescope pointing stability is directly related to attitude 
control system performance and attitude control bandwidth is constrained by the 
lower structural modes. The "quiet" spacecraft requirement also is applicable to 
gimbaled pointing system for the long, passive ranging links since wideband 
pointing error information is not available due to severe power limitations and 
a low laser pulse rate, and hence no compensation for high-frequency pointing 
jitter can be provided. 

Body pointing is the recommended approach for implementing the very stringent 
laser ranging system pointing requirements. Although attitude control accuracy 
requirements are tighter with body pointing, it is felt that this comparative 
disadvantage is heavily outweighed by the simpler mechanization and the absence 
of the performance-limiting gimbal related problems due to angle encoder readout 
error and gimbal dynamical errors. The body pointing approach contains fewer 
components, probably is lighter, and alleviates thermal design problems attendant 
to high energy dissipation rates inside the laser telescope gimbals. It also 
is easy to implement redundancy with body pointing, whereas, an equivalent level 
of redundancy is not feasible with the gimbaled system. 

4.4 POINTING ERROR BUDGET FOR LONG, PASSIVE LINK. Preliminary pointing error 
budgets for the recommended body pointing system are presented in Tables 9 and 
10 for the long, passive ranging links. The open-loop pointing requirements 

for the acquisition and tracking phases are 5.0 yrad and 1.0 yrad respectively. 

All pointing error contributors are assumed to be independent and normally 
distributed. 

A 1.5 yrad pointing error allocation is specified for the effects of HATRS 
and target position errors during acquisition. This corresponds to approximately 
20 meter (3a) along track and cross track position errors for the HATRS and the 
target if acquisition is performed with the target on the horizon. This ephemeris 
error allocation is about three times higher than required if the 3 meter (la) 
prediction capability projected by Goddard is realized. 
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TABLE 9 

ACQUISITION POINTING ERROR BUDGET FOR LONG, PASSIVE LINKS 
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TABLE 10 

POINTING ERROR BUDGET FOR TRACKING PHASE 
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On-orbit calibration of the aligntnent between the laser telescope boresight 
axis and attitude reference axes is required after the operational orbit has been 
achieved in order to remove static misalignments attributable to manufacturing 
tolerances and boost environment effects. Then, periodic on-orbit calibration 
must be performed during the mission to measure alignment shifts due to thermal 
gradients and structural creep. The pointing error budget for the acquisition 
phase contains a 2.5 yrad allocation for residual calibration error, low-frequency 
alignment shifts between calibration periods, and high-frequency dynamics such 
as structural vibrations. 

A 1.0 yrad error allocation is specified for attitude control errors. This 
includes the effects of spacecraft and environmental disturbances plus actuator 
nonlinearities and asymmetries. The remaining allowable pointing error during 
acquisition is allocated to the attitude reference system. The ARS error com- 
ponents include star sensor error, gyro errors and computational error, with the 
primary gyro error sources being drift uncertainty, scale factor stability, and 
input axis alignment stability, 

A preliminary pointing error budget for the tracking phase is presented in 
Table 10 for a 1.0 yrad tracking requirement. It is assumed that static and 
very low-frequency boresight errors are estimated using a conical scan boresight- 
ing technique following target acquisition and prior to the start of the tracking 
phase. Hence it is assumed that the only contribution of static boresight mis- 
alignment and low-frequency boresight alignment shifts to tracking error is a 
calibration or measurement residual. The boresight error estimator also measures 
the pointing error due to HATRS and target location errors. There are short- 
period ephemeris propagation errors but their effects on pointing error should 
be small in a "root-sum-square" sense compared to the 1.0 yrad tracking accuracy 
requirement. Thus, a preliminary allocation of 0.33 yrad is made for the residual 
calibration error (i.e. estimation error) plus short-term ephemeris effects. 

A 0.33 yrad error allocation also is specified for attitude reference system 
to laser ranging system misalignment due to short-period thermal effects and 
structural vibrations. This imposes stringent requirements on the precision 
optical assembly which is the thermo-mechanical interface between the ARS sensors 
and the laser telescope. The HATRS attitude control accuracy requirement is 
defined as 0.50 yrad in the tracking error budget. This defines the allowable 
pointing error due to angular deviation of the spacecraft relative to a desired 
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reference as defined by the attitude reference system. The achievable attitude 
control accuracy is influenced by HATRS and environmental disturbances, control 
actuator performance, and the degree of interaction of the attitude control 
system with spacecraft structure. 

The remaining approximately 0.70 yrad allowable pointing error is allocated 
to the attitude reference system. The 0.70 prad error allocation includes effects 
of star sensor error, gyro error, and computational errors. In order to satisfy 
this requirement, third generation gyros and star sensors with an overall accuracy 
of 0.25 prad or so are required. The availability of adequate gyros should not 
be a problem, but the availability of star sensors with the accuracy just quoted 
is less certain. 

The pointing error budgets for the Laser Ranging System acquisition and 
tracking phases impose very stringent requirements for ephemeris prediction, 
attitude determination, attitude control, and spacecraft structural and thermal 
design. However, it is felt that the 5 prad open-loop pointing requirement can 
be met by the early 1980's with reasonable extrapolations from the current state- 
of-the-art in ephemeris prediction capability, attitude sensor technology, and 
attitude control system performance. Achievement of a 1 prad pointing capability 
for tracking also appears to be feasible if an acceptable boresighting technique 
can be developed to estimate static boresight errors following target acquisition. 


4-22 



5.0 PERFORMANCE ANALYSES 

Many types of system configurations were examined, and performance analyses 
made, prior to selecting the preferred systems configurations, described in 
Section 3. In this section, the performance analyses for these alternate con- 
figurations are described and compared to the performance predictions for the 
preferred systems. Additional performance analyses for the preferred system con- 
figurations are also presented. 

5.1 HAIRS RANGING MEASUREMENTS. The first, and most fundamental, concern is to 
determine the amount of received signal energy necessary to meet the rras range 
measurement accuracy requirements. This can be expressed in several ways. The 
most useful measures are in terms of detected photoelectrons per pulse and 
detected photoelectrons per second. The former is used for systems which make 
range estimates on a single pulse or average over a few pulses. The latter 
description is more appropriate for systems which average over many pulies 

to obtain range estimates. 

All of the systems considered in this study employed side-band ranging 
techniques, due to unresolved problems of coherent carrier detection at or near 
visible frequencies. Side-band ranging techniques rely on estimating the range 
by measuring the phase delay of the modulation waveform which is imposed on the 
(optical frequency) carrier. When the modulation waveform is periodic, and the 
period is less than the maximum possible propagation time, a range ambiguity con- 
dition exists, which must be resolved in some manner for true range measurements. 
When the modulation waveform period is significantly larger than the range un- 
certainty, the ambiguity can be readily resolved. However, when the waveform 
period is much less than the range uncertainty, the ambiguity cannot be easily 
resolved, and the system is limited to estimating the rate-of-change of the phase 
path length, commonly referred to as Doppler frequency measurements or range-sum- 
rate measurements. 

These techniques can be used in either a closed-loop system or an open-loop 
system. In this context, an open-loop system implies that the receiver and the 
transmitter are physically separated, and the reference (transmitter) phase is 
estimated at the receiver (as opposed to measured). An open-loop system, by 
way of illustration, can be configured by providing precision clocks at both 
the receiver and the transmitter, and requiring that the transmitted modulation 
waveform reflect the state of the transmitter clock at the time of transmission. 
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The receiver then estimates the phase path length by subtracting the estimate of 
the transmitted clock state from the receiver clock state to obtain a propogation 
time estimate. The preceding arguments on waveform period apply also to this 
system, thus, either open-loop range estimates or range-rate estimates can be 
made. 

5.1.1 Ranging with Passive Targets. Ranging measurement accuracy is fundamen- 
tally a function of the signal strength, the magnitude of the system noise, the 
net system bandwidth, and various uncompensated error sources within the system. 

At optical frequencies, with high quality detectors and optical systems, the 
received signal is said to be shot-noise limited, i.e., individual photoelectron 
emissions are detectable. Under these conditions, the detector output can be 
considered as a series of impulses, each signifying a photoelectron event, 
which are passed through a filter representing the impulse response of the 
detector and the amplifiers that precede the signal processing electronics. 

Under ideal circumstances, the time of arrival of each detected photoel ectron 
event can be determined, and used to estimate the phase of the received signal 
modulation waveform. If the estimation process is optimum, in some sense, the 
expected accuracy (as a function of received signal energy) of this system is a 
lower bound on the achievable range estimation process in a realizable system. 

In Appendix E, a maximum likelihood estimation process is defined for the gener- 
alized modulation waveform, and evaluated for the specific case where the wave- 
form is a series of pulses. Several sub-optimum estimation algorithms are also 
synthesized and compared to the maximum likelihood estimation algorithm for the 
pulse modulation format. 

Two basic approaches to range measurements using pulsed systems are the 
single pulse estimates, and estimates derived from the average of multiple 
pulses. For the single pulse case, an estimate of the range or propagation 
time is made for each received pulse which is detected. For the multiple pulse 
case it is necessary to average over several received pulses prior to 
making an estimate, hence we must include a predictor module, since the pulses 
are spaced in time. In the simplest form, the signal waveform has a periodic 
component, and the predictor module is a simple phase-locked loop. The loop 
bandwidth is chosen on the basis of the link dynamics. However, when the pulse 
rate or the usable loop bandwidths must be reduced below the level necessary 
for tracking with a simple (second order) phase-locked loop, a more powerful 
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predictor module must be employed. This is possible if the link dynamics are 
governed by a known set of laws. In general, this is the case for ranging 
between two satellites in orbit about the earth; the laws of motion are well 
known, only the initial conditions and the exact magnitude of the gravitational 
harmonic terms and the disturbances are unknown. The results of these uncertainties 
are what we refer to as prediction noise, i.e. , errors in predicting the pro- 
pagation time. In general, the effective bandwidth of the prediction noise is 
much less than the effective bandwidth of the link dynamics. 

CW Mode-Locked Laser Analysis - A cw mode-locked laser was initially considered a 
viable approach for ranging applications. Conceptually, either dual apertures or 
an interrupted a«/ mode could be used. However, the link was found to be seriously 
background limited, which resulted in requiring excessive transmitter power, a 
very large aperture, or a very large cross-section target. Further, if the target 
has a finite dimension in range, the return pulse is broadened, and the effective 
background level increased still further. These considerations led to removal of 
the cw mode-locked laser from the viable candidate list. 

Pulse-burst Mode-locked Laser Concept - Initially, we considered the pulse- 
burst, mode-locked laser technique as an alternative to the CW mode-locked laser, 
which could, conceivably, provide the advantages of narrow pulses, but with a 
reduced duty cycle to reduce the effects of background. Conceptually, the laser 
would be very similar to the CW mode-locked laser, with the exception that the 
laser would be pulse pumped with a low duty cycle, high energy lamp. The 
resulting output would be a series of mode-locked pulses, with a duration of 
about one microsecond (-^200 pulses/burst). The pulse bursts would be repeated 
periodically, effectively resulting in a sampled system approach to maintaining 
the phase-lock- loop in operation. A relatively high pulse burst frequency is 
necessary to maintain phase-lock, since it would be very difficult to acquire 
phase-lock during a pulse burst. Since the pulse-burst laser does not now 
exist, and significant development would be required to ascertain feasibility, 
and no significant advantages (over the single pulse techniques) could be found, 
this measurement concept was also judged inappropriate. 

Single Pulse Laser Measurement Concepts - Two types of lasers were included 
in this category, a Q-switched, cavity dumped laser, with a typical pulse- 
width of -^4 to 6 ns FWHM, and a Q-switched, mode-locked, cavity dumped laser with 
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a typical pulsewidth of 0.1 ns FWHM. The two lasers are very similar, the basic 
difference is the addition of a mode-locking crystal in the cavity for the short 
pulse laser. Clearly, addition of an intra-cavity element will reduce the output 
power, thus we assumed 10 dB less output power for the short pulse laser than the 
long pulse laser to account for this difference. Since the short pulse 
laser does not exist even in breadboard form, the conservatism of this assumption 
is unknown. 

Two basic measurement concepts were investigated for these two low PRF lasers. 
First, we investigated the performance of single pulse estimation techniques, where 
each received pulse is processed to make a range estimate. Subsequently, we 
analyzed a split gate rate tracking technique. This analysis was also applicable 
to angle tracking when using a quadrant detection technique. 

The single pulse estimation techniques were based on the concept of measuring 
the time at which a transmit pulse is emitted, and estimating the time of arrival 
of this pulse after reflection from the target. Even with an ideal target, an 
infinite resolution clock, an infinite bandwidth detector, and a perfectly stable 
laser, errors in estimating the time-of -arrival of the reflected pulse will be 
encountered due to the statics of the detection process (shot noise). Appendix E 
describes an analysis conducted to determine the optimum detection accuracy, 
and to evaluate the performance of several suboptimum, but potentially more 
readily implemented, time-of-arrival estimation techniques. The fundamental 
concept was that a signal pulse was contained within a finite time "window", and 
the problem was to estimate its location within this window. A number of 
approximations were necessary to arrive at usable answers, so a simulation was 
performed to verify the validity of the results. The simulation results matched 
the theoretical predictions within the expected error of the simulation under the 
conditions examined (N^ >12 pe/pulse). 

At lower signal levels, however, the accuracy is expected to degrade signifi- 
cantly, since there is an increasing probability that few or no signal photo- 
electrons will be detected. In a practical system, we would impose a threshold 
criterion, such that no time of arrival estimate would be made unless the detected 
signal were above this threshold. In order to determine the effect of a threshold 
limit on the variance of the estimation error, we will begin with Equation (8) of 
Appendix E, If the expectations are evaluated conditioned on the occurance of L 
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or more photoelectron events in the pulsewidth, we find, for the maximum likelihood 
filter and a raised cosine pulse shape. 


ML 
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This is exactly the same expression obtained previously (Equation 31, Appendix 
E) with a multiplicative threshold function. For the mismatched filter case 
(correlator). 
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This is also exactly the same expression obtained previously (Equation 32, 
Appendix E) with the multiplicative threshold function. If we choose L=1 , the 
threshold function reduces to. 


g(T) = 1 



+ n^T) 


( 20 ) 


Then, for nj^ = 0, and small, 
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Simply stated, these expressions show that if at least one photoelectron 
event occurs, the variance of the time-of-arrival estimation error will not exceed 
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2 7 

(T/2tt) for a maximum likelihood estimator with zero background, or 2(T/2 tt)^ for 

a correlator or matched filter. Thus, for a 4 ns FWHM pulse, the limiting rms 
error, assuming a ML estimator, would be 1.275 ns rms (19 cm rms); for an 0.2 ns 
FWHM pulse, the error would be .064 ns rms (.95 cm rms). Clearly, for the short 
pulse laser under these conditions, detection probabilities will govern the require- 
ments for signal energy. 

For the long pulse laser, however, " (4 /tt(. 133)) = 91.2 pe/pulse is 

required to achieve the accuracy requirement. This implies, of course, an ultra- 
stable pulse shape, optimum processing, and an ideal target. 

As noted in Appendix F, the targets considered in this study are far from 
ideal. Two significant effects are noted, pulse stretching and scintillation 
resulting from coherent interference. When the target has finite range dimensions, 
such as the LAGEOS, coherent interference may result in a random variation of the 
apparent mean of the return pulse. When the transmit pulse is short with respect 
to the effective target dimension, the individual pulse shape can be grossly 
distorted by the coherent interference, and the mean pulse shape will be distorted 
by the impulse response of the target. Distortion of the mean pulse shape can be 
countered, to some extent, by using optimal filtering in the receiver, however, 
note that this pulse distortion results in an increase in the effective pulsewidth 
and time-of-arrival estimation error, requiring an increase in mean signal level 
to maintain acceptable performance levels. 

Scintillation resulting from coherent interference can have a much more profound 
effect on the performance of the system. These effects arise from two distinct 
causes. First, consider the case where the target is a planar array of cube- 
corners at, or near, normal incidence. In Appendix F we show that the optical 
pulse energy (N^) of the return signal from such a target has a random variation, 
with an approximately exponential probability distribution. The expected value of 
this pulse energy is simply the mean energy computed assuming noncoherent summing 
of the returns from each cube-corner. Since a maximum likelihood estimator requires 
knowledge (or estimation) of the signal and background levels, this technique 
appears to be inappropriate when the signal energy is a pulse to pulse random 
variable. The correlator approach, however, does not require estimating these 
parameters, and is, therefore, a usable process. 

The effect of coherent interference on the variance of the TOA estimation error 
can be determined for the zero background, unity threshold case by averaging 
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equation (19) over the probability space of N^, resulting in. 



an (1 + Ng) 


( 23 ) 


For comparison purposes, consider a 4 ns FWHM pulse. Inverting Equation (23) 
numerically to determine the required to achieve 0.133 ns rms estimation error 
results in = 1309.3, an 9.6 dB increase in required mean signal energy compared 
to a nonscintillating target. 

Next, we consider the case where a larger threshold criterion (L > 1) is used. 
To accomplish this objective, we return to equation (8) of Appendix E, and evaluate 
the expectations conditioned on the occurance of exactly K photoelectrons in the 
pulse width. The result is. 



Now, suppose we set the threshold for accepting a measurement large enough 
such that the probability of obtaining a measurement is negligible unless N^ is 
much greater than n^^T. Then, the variance of the error, given K, is effectively 
dependent only on K and T. We could therefore choose L such that if a measure- 
ment is accepted, the variance will not exceed the accuracy goal . For the 4 ns 
pulse, this results in L = 183. If we assume that the background is negligible, 
the probability of obtaining K pe events from the scintillating target was shown 
in Appendix F to be given by. 



Hence, the detection probability is. 
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Now, if we assume the system can emit 10 pps, and we will accept an average 
of one measurement per second, Pp = 0.1, and = 79. However, since it requires 
10 pulses, on the average, to obtain one measurement, the average required energy 
per measurement is 790 pe. Using the concept of minimum energy per measurement as 
a criterion, we select Pp = 0.366872 and = 182 pe/pulse, which results in 
'^496 pe/measurement, on the average. Thus, at best, we find the scintillation has 
cost '^>4 dB in link margin for the long pulse laser. The short pulse laser concept 
is also affected, although the criterion for selecting a threshold is more likely 
to be minimization of false alarms due to background. 

The next phase of the analysis will consider a split-gate tracking loop 
as an alternative to single pulse detection techniques. Since the link dynamics 
are such that a 10 pps system cannot track unaided, a range prediction module 
is required to drive the loop. The phase-1 ocked-loop is used to correct low 
frequency errors and static biases in the range predictions. Consider a phase- 
locked-loop constructed using a split-gate integrator to derive an error 
signal. The split-gate integrator simply counts the photoelectron events in two 
disjoint segments of the signal field, and uses the difference in the two 
counts as an error signal to guide the repositioning of the partition boundary 
to the center of the signal field. Since the events are not weighted by time 
of arrival, the technique is suboptimum. Background photoelectrons do not bias 
the estimate, provided the gates are balanced. When the loop is locked up, the 
center of the window, on the average, will coincide with the mean centroid of the 
received pulse, thus the time of transition from early gate to late gate is the 
estimate of the time-of-arrival of the centroid of the received pulse. 

The analysis begins by noting, for a fixed mean number of photoelectrons 
detected per pulse, that the number of photoelectrons detected in each half of 
the gate are random variables, governed by the Poisson distribution. For simplic- 
ity, let us define the observed count (for one pulse) in the "early" gate as KE, 
and KL for the late gate. Then, assume that the pulse is positioned in the 
vicinity of the partition, such that Z(1 +e)/2 is the pulse energy in the late 
gate, and e is a positioning error term, -l>e<l. Define V = KE-KL. The next 
step is to determine the mean and variance of V. These are most easily determined 
in the transform domain. It can be shown that 0^(w) and 0|_(w), the characteristic 
functions for the early and late gates are: 

= exp Z^(e'^“-1) 4 >l(o))= exp Z^Ce'^’"-!) 
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Then, the characteristic function of V = KE + (-KL) is the product of (u) and 
((>l_ (-U)), since, 

£ {gjco(-KL)J, ^ |e-ja)(KL)| 


Thus, 


/ 

• \ 


<f>y(t»)) 

= exp /z^(e’^“-l ) ■ 

f Zg (e"J“-l )j 

Since, 

E Ivl^ 

_ 1 rd^vCu)] 

j L J ^ “ 0 


and 


L d^9v(a)) 
1 do3^ 


we find. 

E ivi. 

= Z, - Z2 = ze 




•= n,jT + Z + (Ze)^ 

(27) 

Then, if 

Z is governed by an exponential 

distribution, as shown in Appendix E, 

fz(z) = 

1/m exp (-Z/m), 



e{v| = / e{v/z} 

. f^(Z) dZ = me 

(28) 




;{v^} = »l,T * m f 

2(me)^ 

(29) 

= n. T + m + 

(me )^ 

(30) 


Then, the expected signal -to-noise ratio in the tracking loop is simply, 
SNR = Nm /( 0 |^T + m), where N is the number of pulses used to form the estimate 
(N = 1/loop bandwidth). Then, the rms tracking error, in the normalized form, 
is simply ep |\/|5 = 1 / ANR. For large SNR conditions (6-9 dB or more), and a 
temporally Gaussian pulse, = t/ 2Tr/SNR, where 2x is the pulse width at the 
60% power points. For zero background, this becomes t/2tt/NM. Assume 

T = 0.2 ns. Then, we see that Nm > 3.53 is required to meet epMo > 0.133 ns. 
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The rms error for a maxitnutn likelihood detector, under zero background conditions, 
never becomes as large as 0.133 ns for a 0.2 ns pulse ~ 0.061 ns) by 

equation 21, if an estimate can be made at all 1 pe detected). However, the 
assumption of zero, or at least negligible, background level is meaningful only 
as long as the observation time, i.e., the period of time over which the return 
pulse is expected is comparable to the pulsewidth. 

In order to determine realistic estimates of required mean signal energy 
under dynamic conditions, a series of simulations were performed for the split- 
gate range tracking system concept. The first series of runs were made to 
determine optimum system parameters, and employed a first order range tracking 
loop, with no range prediction errors. 

The simulation results indicated that an AGC function of the form, ^^(KE+KL) , 
yielded the best results, provided the AGC amplifier gain was clamped to a maximum 
value of unity. 

The simple first order loop was then replaced with a second order loop as 
shown in figure 10. Conceptually, the transmit pulse departure time is used to 
estimate the arrival time of the reflected pulse. This estimate is summed with 
the output of the range tracking loop, and delivered to the gate timer. The 
sequence is that the early gate is opened prior to the expected arrival time. 

At the expected arrival time, the early gate is closed and the late gate opened. 
Subsequently, the late gate is closed. Then, the sum and difference values of 
the two integrator outputs are used to update the range tracking loop. 

The return pulse arrival time prediction will be in error by some amount, 
due to factors such as quantization noise, interpolation errors, unmodeled (or 
incorrectly modeled) perturbations of the orbit path(s), and errors in the initial 
conditions. The unmodeled perturbations were considered the most difficult error 
source to control, since hard data on the higher ordered harmonics of the 
Earth's gravitational field simply do not exist at this time, and, in fact, 
one of the scientific objectives for this ranging system is to obtain this 
data. In view of the lack of hard data, we elected to estimate the unmodeled 
perturbations and then apply a suitable safety factor to the model and examine 
the effects on ranging system performance. 

The first step was to compute the frequency spectrum of gravitation forces 
on a low altitude spacecraft due to a small point mass concentration on or near 
the surface of the Earth. This was done by computing the radial acceleration 
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force, as a function of time in the orbit plane, and then taking the Fourier 
transform of the data. We used the fast Fourier transform (FFT) technique, so 
the output is a line spectrum, in terms of harmonics of the orbit period, i.e. , 

a(t) = f; D (31) 

n= -00 


The n = 0 term was dropped, and the remaining terms normalized by dividing 
by the magnitude of the first harmonic term. The resulting data, plotted in 
Figure 13, shows that the magnitude of the coefficients (D^) decreases in a nearly 
exponential manner with increasing frequency. 

If we assume the radial displacement function can be adequately modeled by 
simply integrating twice, we obtain, 


d(t) = C t “n 

11= n ^ 

n^O 

The next step in this process is to sum the contributions of many coplanar 
mass concentrations at various locations. Thus, each term in the harmonic 
expansion is, 

^ j27rnt/T V jn<l>0 

n2 ® ^ \ % /e 


Where m.. is the mass surface concentration and f. is the angular location with 

vJ \3 

respect to the orbit plane. If we assume m. and <j>., are independent random 

J J 

variables, we can express the autocorrelation function as. 


where 


R(T) = 



j2T7nT/t 

e 
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Thus, assuming independence, 


R{T) = f; ea2TrnT/x (33) 

n= -» \n V 

The next step is to assess the magnitude of the constant C. The simplest 
approach is to compare the differences in position for a spacecraft under the 
influence of two slightly differing gravitation models. Figure 12 presents the 
results of one such internal study, performed in 1971. The traces show 
reasonably periodic deviations, with an estimated fundamental (one orbit period) 
component of about + 60 feet. Improvements in models and increased number of 
gravitational harmonics considered should significantly reduce this deviation, 
perhaps to more on the order of a few feet by the 1980 time period. Thus, if 
we assign a value of "C" equivalent to 60 feet, we are probably conservative by 
an order of magnitude. 

Figure 13 shows the simulation results, for a 200 ps laser system, obtained 
for three values of C, equivalent to 0, 60 feet, and 120 feet, with mean signal 
energies of 1 , 2, and 5 pe/pulse, a scintillating target, and loop bandwidths on 
the order of 0.5 Hz. This bandwidth setting appeared to be near optimum, as 
wider loop bandwidths had markedly increased variances, while narrower loop 
bandwidths experienced increasing bias terms. At 1 pe/pulse signal levels, 
increasing "C" to 120 feet resulted in one loss of lock in ten trials, indicating 
near marginal operation. Loss of lock was not observed at any other condition 
examined. Thus, we can confidently predict that the system will function 
satisfactorily at 1 pe/pulse, provided the prediction error does not materially 
exceed the nominal 60 foot peak (120 foot peak-to-peak) level. 

5.1.2 Ranging With Active Targets : The target was assumed to have a receiver 

and a transmitter, and a precision clock. Both the ranging system and the target 
would emit pulses repetitively, and record in local clock time, the transmitted 
pulse departure time and the times when received pulses were detected. These 
disjoint event time recordings are used for subsequent processing to extract 
range data. 
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When the active target contains a receiver and a transmitter, the detection 
problem is much the same as for passive target ranging, except that our 
ability to predict the arrival time of the next received pulse is degraded 
considerably. First, since we are not obtaining real time ranging data, 
ephemeris errors are not reduced and the range will be changing in a different 
manner than expected. Second, the receiver does not have knowledge of when 
the pulse was actually emitted. If the transmit pulses could be emitted at 
very precise intervals, it would be possible to reduce the uncertainty by 
extrapolating from the last received pulse. However, this is the only known 
situation where precise control of emission time is advantageous. Consequently, 
we elected not to require precise emission time control, and to use a time of 
arrival estimation technique which does not require tight range gating. When 
operating with an active target, the links are not seriously power limited, 
thus a suboptimum TOA algorithm does not impose a serious penalty. A modified 
version of the weighted counter with an adaptive threshold, discussed 
in Appendix E, is shown in figure 14. The modifications are intended to adapt 
this concept to ungated operation. When a signal is encountered, the operation 
is identical to the circuit of Appendix E. When only background is present, 
however, an AGC clamp is applied to reduce the false alarm probabilities to 
acceptable levels. Also, the simple integrator has been replaced with a 
sliding window integrator to permit continuous scan. 

5.2 POINTING CONTROL SYSTEM ANALYSIS. The fundamental problem in the pointing 
control system for the HATRS terminal is to control the direction of propagation 
of the very narrow (5 yrad) transmit beam to ensure that the transmit energy 
illuminates the target. The sources of error in this control process include 
prediction errors (ephemeris uncertainties) , attitude reference uncertainties, 
and thermal and mechanical errors. Many of these uncertainties have a low 
frequency content, i.e., they change in magnitude very slowly, and thus can be 
combatted with active signal sensing techniques. The energy of signal reflected 
from the target can be used to form an estimate of the pointing error if the 
transmit beam is dithered (systematically varied) about the nominal expected 
direction. These estimates can then be used to bias the pointing commands to 
center the transmit beam on the target. The accuracy of this centering 
process is dependent on the mean reflected signal energy, the type and magnitude 
of the dither process, and the bandwidth of both the disturbance spectrum and 
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the filter used to smooth the error estimates. The pointing control concept 
requires very precise open-loop pointing, which is augmented with heavily fil- 
tered pointing error estimates to reduce static pointing errors. The attitude 
reference system inputs the estimated boresight direction vectpr to the ephemeris 
computer, which computes the desired boresight direction vector, and outputs a 
pointing error estimate to the pointing control loop filter. The loop filter sums 
the computer error estimate the dither command, and the filtered bias error 
signal and outputs acceleration commands to the pointing control actuators. 

The open-loop error can be expected to contain very low frequency error 
sources (thermal drift), moderate frequency error sources (ephemeris prediction 
errors), and wide-bandwidth error sources (vibration). The boundaries of these 
domains are related to the transmit pulse rate. The very low frequency error 
sources are highly correlated over many transmit pulse periods. The moderate 
frequency error sources are characterized by auto-correlation functions which 
decay significantly in a finite number of pulse periods. The wide bandwidth 
error sources result in errors which are virtually uncorrelated from pulse-to- 
pulse. We note that the pointing error estimation process is noisy, and that 
this noise is virtually uncorrelated from pulse-to-pulse. The pointing control 


5-18 






system filter acts as a low pass filter to the estimation error noise, and 
as a high pass filter to the pointing error noise. The optimum control system 
filter is the filter which minimizes the residual mean square pointing error, 
which is simply the sum of the estimation error noise which passes through the 
system filter and the pointing error noise which is outside the system filter. 

The first step in determining the performance of the system is to select 
a dither process. We evaluated two techniques and found virtually identical 
performance. The first dither process used a conical scan about the nominal 
pointing direction. The second technique employed a bang-bang type of dither, 
where the beam was deflected by a small angle (Ae) is each of four direc- 
tions (+x, -X, +y, -y). The latter approach was chosen as the best approach 
since only four pulses were required to form a pointing error estimate. 

The transmit beam has a nearly Gaussian spatial intensity distribution, 
thus the mean reflected signal energy (m) is related to the pointing error at 
the instant a pulse is transmitted by 

m = m^ exp -(x^ + y^)/R^ 


where m^ is the received signal energy when x = y = 0, x and y are the coordinates 
of the pointing error, and R is the transmit beam radius at the e"^ power point. 

X, y, and R may be expressed in linear dimensions (meters) or angular dimensions 
(radians), or simply normalized in terms of the beam radius. The latter approach 
simplifies the notation and is used subsequently. 

Now, assume that the mean prediction pointing error (xq, yo) is constant for 
four successive pulses. Denote the receiver output for these four pulses as 
Z^v, i =1, 2, 3, 4. We will assume the receiver is an ideal photoelectron 
counter. Thus, the Z-j are independent, Poisson random variables. Next, assume 
that the actual pointing direction is offset from the nominal direction by an 
amount Ax for the first pulse, -Ax for the second pulse. Ay for the third pulse, 
and -Ay for the fourth pulse. Then, let gx and gy be defined. 
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Consequently, g and g are random variables with means n and n , and variances 

rt X X X y 

a)C and given by 

n = m„ le'^^o Ax)2-yo2 _ g-(Xo-Ax)2-yo2 

X 

aZ = 2nbT + 

for X, and similarly for y, and nj^T is the background photo-electrons/gate. 

These equations can be rewritten, 

n^^ = 2 m^j sinh (ZX^ax) e'^^o^ 

Ox^ = 2n|^T + Zrn^j cosh (2 Xq AX)e"^^o^ ■*■ Xo^ 

When x«» y„ and Ax are small, n^ :i -4m^ x^ Ax, thus the estimate of x is 

u \j X 0 0 0 

simply, 

Xq = -n^/Am^Ax 

and the variance of the estimate of x^ is, 

®xo^ - ax2/(4moAx)^ 

Thus far, however, we have not employed all of the information available from 
the return signal. If the receiver field-of-view is divided into four or more 
non-overlapping detection areas, we can obtain an estimate, from each received 
pulse, of the received signal angle-of -arrival with respect to the receiver bore- 
sight. If the error between the transmit boresight and receive boresight is 
small, the angle-of-arrival estimates are a good approximation of the transmitter 
pointing error. Ignoring signal to noise ratio considerations, the usable band- 
width of the angle-of-arrival estimation process is four times the bandwidth 
of the transmitter pointing error estimation process. 

The simplest solution to the angle-of-arrival estimation process is to 
divide the receiver field-of-view into four quadrants. The angle-of-arrival 
estimates are based on the detected signals from each quadrant, f(y) = 

A+B-C-D, and f(x) = A-B-C+D, where the quadrants are lettered counter-clockwise 
from the conventional first quadrant in an x, y coordinate system. As noted 
previously, each signal is a Poisson random variable, whose mean is determined 
by the amount of the target image energy on that quadrant. Now, A, B, C, and D 


-{xQ+Ax)2-yo^ ^ g -(xq-aX)^ - yo^j 



5-20 



are independent (Poisson distributed) random variables, thus the sum of any two 
is also a Poisson distributed random variable, with a mean rate equal to the 
sum of the mean rates. Consequently, we can use the analytical form developed 
for the split-gate range tracker, where e, the division parameter, now refers to 
the spatial (or angular error) division parameter. 

As discussed previously, the split gate tracker was evaluated by simulation, 
and we found under typical conditions that at one signal pe/pulse, it was 
possible to achieve 'v 0.03 ns rms tracking error for a temporarily Gaussian 
pulse of 0.2 ns width at the e"2 power points. Thus, the achieved performance 
is 0.15 ns rms per ns pulse width, which would scale to 0.15 yrad rms per 
urad image diameter at the detectors. For a one meter diameter collector, 
the diffraction limited image diameter is 2.44 X/D = 1.3 yrad, (at the first null). 
Thus, with diffraction limited optics the rms angle of arrival estimation error 
would be less than 0.2 yrad, which is well within the 0.5 yrad budget. 

In the split gate tracker simulation, prediction noise was a major 
consideration in choosing the tracking loop bandwidth. For angle of arrival 
estimation, prediction noise is several orders of magnitude smaller in angular 
units and, hence, virtually negligible. The other pointing error sources, such 
as mechanical vibration, must simply be controlled to acceptable limits. Linder 
these assumptions, therefore, we can confidently predict that if the range 
tracking functions can be accomplished, angle of arrival tracking to 0.4 yrad 
rms can also be accomplished. 
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6.0 CONCLUSIONS 


Laser ranging from a space-borne terminal to active and passive targets in 
space or on the ground has been shown to be feasible, and that 2-cm rms range 
measurement accuracy is achievable. Two laser ranging systems have been synthe- 
sized, one to meet the requirements for very long range measurements, and the 
other for relatively short range measurements. The very long range links 
('\-30,000 km) are severely power limited, requiring precision optical systems, 
narrow transmit beamwidths (5 yrad), and sophisticated signal processing 
techniques if the ranging accuracy goals are to be met. These requirements can 
be met with implementable systems provided the target is suitably configured. 

A LAGEOS type of target was evaluated in detail, and found to contribute 
significantly to the uncertainties in range measurement accuracy. Both narrow 
pulse ('v/200 ps) and wide pulse (4 ns) lasers were considered for the terminals. 

If the narrow pulse laser can be developed for this application, the inherent 
ranging accuracy capability is well within the sub-centimeter region, and will 
be most seriously limited by the target characteristics and signal processing 
implementations. These limitations do not appear insurmountable, although 
development is required. The wide pulse laser concept theoretically capable 
of meeting the range estimation accuracy goals, provided the output pulses 
are very stable and repeatable. The required precision pointing control and 
alignment control functions are well within the state-of-the-art projected for 
the 1980's time frame. 

A laser ranging system concept for the relatively short range (300-500 km) 
measurements was evolved for a Space Shuttle or Spacelab experimental mission, 
with ranging to ground based retro-reflectors the fundamental requirement. For 
this application, limitations in ephemeris prediction accuracy and available 
target viewing time resulted in significant differences from the long range 
terminal concept, although many of the components are identical. This concept 
is considered feasible and implementable, but several unresolved questions 
remain to be studied before the requirements for the experiment can be 
established. 

Measurement strategy was discovered to be a major factor affecting both the 
experiment results and the design of the laser ranging system. Detailed evaluation 
of the effect of measurement strategy on the target grid relative position 
reconstruction accuracy is needed to determine whether simultaneous (multiple 
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beam) ranging is needed, or if sequential ranging measurements can meet the 
experiment objectives. 

Another unresolved question of potential significance concerns the integration 
of the experiment and the Shuttle/Spacelab. The Shuttle requirements and capabilities 
are not ideally matched to the laser ranging experiment. Factors such as c.g. 
uncertainty, structural flexure, empheris uncertainties and vibration spectrum 
may be very significant to the design and feasibility of the laser ranging 
experiment, and deserve more attention than could be allocated in this study. 
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APPENDIX A 


FAR FIELD RADIATION PATTERNS - LAGEOS CUBE-CORNERS 

1.0 INTRODUCTION 

The far field radiation patterns and the programs used to generate this 
data are included in this appendix. The data file generated in this activity 
was also used by the program, discussed in Appendix B, to compute the return 
pulse shapes. 

2.0 RADIATION PATTERN ANALYSIS PROGRAM 

The program used to compute the radiation pattern and gain was written in 
Fortran IV, for MDAC's Sigma 9 time sharing computer system, and was used in 
the batch (off-line) mode. The program is listed in the following pages. 


RAY5 07: 59 MAR 06, ’7:3 


100$PS=2000 

110 COHPLEX AXC512),AY(100,1000,ELEMT 

120 COMPLEX ALEMT 

130 DIMENSION DPC60,TRANC116) 

140 COMMON AIN,AROT,FRAD, FR2,ADROT, ALAM, AGRAD, IC0N,3MSPRD 

150 NAMELIST FDl ACM, ALAM, AIND, AD IHED, IREMD,MMAX 

160 Pl=3. 1415926535898 

170 NMAX=256 

180 NDIM=50 

190NDIM2=100 

200 FDIACM=3.8 

210 ALAM=.53 

220 AIND-0 

230 ADIHED=7.3 

240 CALL OPENFIC’RAYDAT‘,1) 

245 CALL 0PENF2C’RAY0UT’ ,2) 

250 IREND=0 
260 1 CONTINUE 
270 INPUT CD 

280 IFClREND.EO,0)GOTO210 
290 CALL CLOSFl 
295 CALL CLOSF2 
300 STOP 

310 210 CONTINUE 
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320 PRINT301^NMAX,FDIACM, ALAM^AIND, ADIHED 

330 301 FORMATC// 'SAMPLE SPACING, FACE DI A, LAMBDA, INC I DENCE, DIM . ' 
3405/5F12.2) 

345 WRITE C2^)FDI ACM, ALAM,AIND^ ADI HED 
350 PRAOa-5«-FDIACM 
360 FR2=FRAD”"2 
370 DELX=FDIACM/5. 

380 DELTIIE=ALAM"1000./NMAX/DELX 
400 AIN=AIND”PI/180, 

410 delz=delx/coscain) 

420 DELYsDELX 

430 BMSPRD=4,75351’*ADIHED 

440 ADR0T=ADIHED/2. 

450 A1NDF=ASIN(.6849«SINCAIN)) 

460 ALP=COS(AINDF) 

470Dl=l./snRTC3.) 

480D2=1./SortC2.) 

490 D3-1./SORTC6.) 

500 AN=cos(.955317+aindf9 
510 AM=Sort(C1.-AN«”2)/2.) 

520 AL=AM 

530 '==FRAD/D2 

540 X0=P«D1 

550 DSTOP=2”P/ALP 

560 DO 2 1=1,100 

570 2P=DELZ"Cl-NDIM+.5) 

580 DO 2 J=l,100 

590 yp=dely"C<J-ndim+,5) 

600RAD=ZP«”2+YP”"2 
610 ELEMT=0 

620 I F C R AD . GT . F R 2 ) G O TO 2 01 

630 X1=X0-YP»D2--ZP»D3 

640 Y1=X0+YP”D2-Z'’”D3 

650 Zl=Xd+2”ZP»D3 

660 X2=AL”DSTOP-XO-X1 

670 Y2=AM”DSTOP-XO-Yl 

680 22=AN”DSTOP-XO-Z1 

690 RAD2=X2'!5V2+Y2«”2+Z2""2 

700 IFCRAD2.GT.FR2)GOT0201 

710 PKA1=6 .2332E3/ALAM"BMSPRD 

720 ifCyp.eo.0)yp=i.e-9 
730 DX=X 1 /AL 
740 DY=Y1/AM 
750 DZ=Z1/AM 

760 IFCCDX. LE. DY).AND.CDY. LE.DZ))PMA2=pHA1“CZP«.86603+YP/2.);GOTO303 

770 I FCCDX. LE.DZO . and. CDZ. LE .D''0)pHA2=PHA1"vP; GOTO 303 

780 IFCCDZ. LE.DXO.AND.CDX. LE.DY))PHA2=-PHAl”r7P«.86603-Yo/2.9;G0T03n3 

790 I FCCDZ. LE . D'^) .AND. CDY . LE . DX) ')PHA2=pHA1"C-Z'’*' . 866 O 3-YP / 2 . 9; 001030 3 

300 IFCCDY. LE.DZ).AND.CDZ.LE.DX9 9PHA2=pHA1«C-Y'’);G0T03O3 

810 PHA2=-phA1“C-Zp«.86603+YP/2.) 

320 303 PHA2=PHA2/ALP 

830 PflA=°HA2”l.E-6 

340 ELEMT = CMpLXCC0SCPHA),SINCPHA9') 



850 201 CONTINUE 

860 AYCI,U)=ELEMT 

870 2 CONTINUE 

880 I SIGN- 

-1 

S'; 0 DO 3 I 

=1^100 

900 DO 20 

J=101,NMAX 

910 20 AXCU>=0 

920' DO ^ U 

=1,190 

930 4 AXCUO-AYa ^ 

940 CALL FFTCNHAX.AX, TSTGN) 

950 DO 5 >) 

= 1,51 

960 5 AYCl 

,.0=AXCJ9 

970 DO 21 

U=52,100 

'930 21 AYn, O = AXC.J+NMAX-100') 

990 3 CONTINUE 

1000 DO 6 

1=1,100 

ion DO 22 

0=101, NMAX 

1020 22 AX<,J).=0 

1030 DO 7 

.) = l,lO0 

1040 7 AXCJ9 = AYC.J, I) 

1050 CALL 

fftcnmax, AX, I S into 

1060 DO 8 

0=1,51 

n70 ELEMT 

= AX CvJ ) “CON JG CAX Cu ) ) 

1080 8 AYC-J, n = ELEMT 

1090 DO 23 

J=52,100 

1100 ELEMT 

=AX CU +NMAX- 100 ) "COM JG C AX C-J+NM AX -10 0 ) ) 

11 n 23 Avr.l . n=ELEMT 

1120 6 COM 

TIMUE 

11 30 FACTORS C 4”^ I ” C 1000 . ” DE LX/ ALAfI) 2 ) 2 

1150 CALL 

STATSCFACTOR, AY, DELTME, TP.AN) 

1155 O' RITE 

C 2,) TRAM 

1160 GOTO 

1' 

1179 END 


1180 

SUBROUTINE FFTCME, X, MASCM) 

II 9 OC 

MMAX= LARGEST VALUE OF fJ TO BE ''’nOCFSSED 

120 9 C 

.‘JONDUMMY DIMENSION M(NMAX0 

121^G 

FOR EXAMPLE, IF NMAX=25 TMFN 

1220 

DIMENSION MC9) 

1230 c 

DIMENSION X(2"''NO 

1240 

DIMENSION XCNE) 

1250 

COMf’LEX X,V/K,N0LD,0 

1260 

ASIGN=NASGN 

1270 

IF (ME- 1)10, 10, 11 

1230 11 

DO 12 1 = 1,15 

1290 

IT=NE/C25^:M) 

1300 

1 11= I 

1310 

I F( 1 1 -1)12, 13 , 12 

1320 12 

CONTINUE 

1330 

GO TO 10 

1340 13 

N=III 

1350 

LX=NE 

1360 

FLX=LX 

1370 

DO 1 1=1, N 


A-,3, 



1380 

1 

MCT)=2-«»CN-I) 

1390 


DO 4 L=1,M 

1400 


NBL0CK=2»"CL-1) 

1410 


LBLOCK=LX/NBLOCK 

1420 


LBHALF= LB LOCK/2 

1430 


K=0 

1440 


DO 4 IBLOCK=l,NBLOCK 

1450 


FK=K 

1460 


V= AST Gr J « 6 .2831 85 3- FK / FL X 

1470 


VJ K = C M P L X ( C 0 5 C V ) , S I N ( V ) ) 

1480 


I S TART=LBLOCK5« C I BLOCK-1) 

1490 


DO 2 I=1,LBNALF 

1500 


J=ISTART+I 

1510 


JH=U+LBHALF 

1520 


n=XCUH)"WK 

1530 


X(JH)=X(J)-0 

1540 


xC'j)=xrj)+o 

1550 

2 

CONTINUE 

1560 


DO 3 1=2, N 

1570 


1 1 = I 

1580 


ifCk.lt.mCD) go to 4 

1590 

3 

K=K-MCI) 

1600 

4 

K=K+M(II) 

1610 


K=0 

1620 


DO 7 U=1^LX 

1630 


IFCK.LT.J) GO TO ^ 

1640 


UOLD=XCU) 

1650 


XC.))=X(K+1) 

1660 


XCK+l)=NOLD 

1670 

5 

DO 6 I=1,N 

1680 


1.1 = I 

1690 

* 

IFCK.LT.MCD) GO TO 7 

1700 

6 

K=K-MCI) 

1710 

7 

■K=K+MCII) 

1720 


I FCAS I GN . LT, 0 , 0) RETURN 

1730 


DO 8 1=1, LX 

1740 

8 

xCi)=xci)/flx 


1750 10 CONTINUE 
1760 RETURN 

1770 END 

2370 SUBROUTINE STATS CF, A, DEL, TRAN) 
2380 COMPLEX AC100,100) 

2390 DIMENSION VALCSO, 3)/TRANCllfi) 

2400 ASP=2.5/DEL 

2410 AMIN=ASP/2 

2420 AMAX=69.5/DEL+AMIN 

2430 DO 1 1 = 1,50 

2440 DO 1 U=l,3 

2450 1 VALCI,J)=0 

2460 DO 2 1=1,100 

2470 II=I 

2480 IFCII.OT. 51)11=102-1 
2490 DO 2 J=l,100 
2500 JJ=J 





2510 IFCJJ.GT.51)Jd=102-J 

2520 ARG2SFLOATCCI I-1>««2+CJ J-l)«»2)+l . E-6 

2530 ARG5S0RTCARG2) 

2540 I FCARG.lt, AMI N)G0T04 
2550 I FC ARG. GT , AMAX+ASP>G0T04 
2560 ID=1+INTCCARG-AMIN)/ASP) 

2570 IFCID.GT.50)GOT04 
2580 X=REALCAa,J))”F 
2590 X2=X5«X 

2600 VALC I D , 19 =VALCT D, 1)+1 

2610 VALaD,2)-VAUaD,2>+X 

2620 VALClD,3>=VALaD,3)+X2 

2630 4 CONTINUE 

2640 2 CONTINUE 

2650 DO 3 1=1/29 

2660 XI=VALCI/1) 

2670 X2=VALCI^2) 

2680 X3=VALCI,3') 

2690 XM=X2/XI 

2700VAR=X3/XI-XM”XM 

2710 STD=SORTCVAR+1.E-20) 

2720 ANG=AM I N+ AS P«c FLOATC O - . 5 9 

2740 ANG=ANG«DEL 

2741 II=4«I 

2742 TRANCII-39=ANG 

2743 TRANCIT-2)=XI 

2744 TRANCM“1)=AL0GCXM9 • 

2745 TRANCM9=AL0GCSTD) 

2750 PRINT 100,ANG,XI^XM,STD 
2760 3 CONTINUE 

2770 RETURN 

2780 100 FORMATCG12.4,f6,0,2G12.4) 

2790 END 

3.0 RADIATION PATTERN DATA 

Two programs were used to process the data, from the preceding program, 
for presentation. These programs and the resulting data are presented in the 
following pages. The first set of data, program GAINPLOT, shows the mean gain 
as a function of point-ahead angle, for incidence angles from 0 (normal) to 
40 degrees. The left-hand number is the point-ahead angle in micro-radians. 

The number on the right-hand side is the gain product in dB. The plotting 
field has 50 increments, with the scale factor shown in the heading. The first 
half of the data is for 0.53 ym wavelength illumination, the second half is for 
1 .06 ym wavelength. 

The second set of data, program CROSPLOT, shows the mean gain as a func- 
tion of incidence angle, for point-ahead angles from 5 to 60 microradians. 

The same conventions and scale factors were used, except that the left most 
column of numbers are the incidence angle in degrees. 
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100 DIMENSION TRANC8), ALISTC18^2'I) 

110 DATA STAR, X/‘"S’X*/,DB/i|, 3^129^/ 

120 CALL OPENFlC't£DAT',l) 

130 DO 1 1=1,18 

140 READC1,)D1,D2,D3,D4 

150 oRINT,D1,D2,D3/D4 

160 D02 il=l,12 

170 READ(1,)TRAN 

180 ALISTa,25«J-l)=TRANC3) 

19 0 A L I S T C 1 , 2 «.J ) = TRAN C 7 ) 

200 2 CONTINUE 
210 D03 J=1,3 
220 3 RE AD Cl, ODUM 
230 1 CONTINUE 
240 CALL C LOS FI 
250 AMAX=0. 

260 AMIN=1.E4 
270 D04 1=1,18 
280 D04 J=1,24 
290 DUM=ALISTCI,.J) 

300 I FCDUM. GT. AMAX) AMAX=DUM 
310 IFCDUM. LT.AMIN>AMIN=DUM 
320 4 CONTINUE 
330 VO =170 

340 VX=INTCDB”AMAX0+1. 

350 DV=VX-VO 
360 SCALE=DV/50. 

370 PRINT100,VX,VO, SCALE 
380 DPA=2,5 
390D010 1=1,18 

391 IFCI .EO.IOOUTPUT’ .53 MICROMETERS WAVELENGTH* 

392 IFCI .Ed.lO)OUTPUT*1.06 MICROMETERS WAVELENGTH* 
400 AINC=5.”CI-1) 

410 I FCl .Gt.9) AINC=A1NC-45 . 

420 PRINT101,AINC 
430 DOIO J=1,24 
440 DUM=DB5JALI STCI , J) 

450 DS=CDUM-VO)/SCALE 
460 1D=INTCDS)+10 
465 I F c I D . LT . 10 ) I D= 10 
470 PA=J"DPA 

480 PRINT102,PA,STAR, ID,X,STAR,DUM 
490 10 CONTINUE 
500 STOP 

510 100 FORMAT ( * MAX , M I N * , 2 F8 . 2/ ' SCALE- ’ , F6 . 2, * / D1 V * 0 
520 101 FORMATC/* INCIDENCE ANGLE *,F4 .0, » DEGREES* 
521g,/T4,*PA*,T65, ’GAIN*) 

530 102 FORMATCF6.1,T10,A,TN,A,T60,A,F6.2) 

540 END 



GAIN'’LOT 


08:08 HAR 06,»75 


3.80000 

.530000 

.000000 

7.30000 

3.80000 

.530000 

5.00000 

7.30000 

3.80000 

.530000 

10,0000 

7.30000 

3 .- 3 ( 000 ^ 

i , 5 ^ i > eeo - 

15 iOOOO - 


3.80000 

.530000 

20.0000 

7.30000 

3.80000 

.530000 

25.0000 

7.30000 

3.80000 

.530000 

30.0000 

7. 30000 

3.80000 

.530000 

35.0000 

7,30000 

3.80000 

.530000 

40.9000 

7,30000 

3.80000 

1.06000 

.000000 

7.30000 

3 . 80000 

1.06000 

5.00000 

7.30000 

3.80000 

1.06000 

10.0000 

7,30000 

3.80000 

1.06000 

15,0000 

7.30000 

3.80000 

1.06000 

20.0000 

7.30000 

3.80000 

1.06000 

25.0000 

7.30000 

3. 30000 

1.06000 

30.0000 

7.30000 

3. 80000 

1.06000 

35.0000 

7.30000 

3.80000 

1.06000 

40.0000 

7.30000 

MAX,MIN : 

201.00 170 

.00 


SCALE = 

. 62 /DIV 




.53 MICRnMETERS ''/AVELENGTH 


INCIDErJCE ANGLE 0 . DEGREES 
'’A 


GAIN 


• J 

5.0 " 

7.5 ” 

10.0 ” 

12.5 " 

15.0 « 

17.5 '* 

20.0 « 

22.5 " 

25.0 « 

27.5 " 

30.0 :: 

32.5 ” 

35.0 

37.5 " 

40.0 « 

42.5 ” 

45.0 « 

47.5 ” 

50.0 « 

52.5 ” 

55.0 « 

57.5 « 

60.0 ” 




20 H .19 

197.27 



19^.71 

191.61 

196.25 

198.89 

200.22 

290.48 

199.93 

199.39 

199.38 
200. ni 
200.46 

200.33 
199.50 
198.08 

196.55 

195.33 
194.63 
193.99 
192.84 

190.38 
187.57 
185.13 
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INCIDENCE ANGLE 5. DEG REE <5 
PA 

2.5 

5.0 

7.5 " 

10.0 « 

12.5 

15.0 ■* 

17.5 « 

20.0 « 

22.5 "• 

25.0 « 

27.5 » 

30.0 « 

32.5 

35.0 ” 

37.5 •' 

40. n ” 

42.5 '* 

45.0 ” 

47.5 « 

50.0 « 

52.5 ” 

55.0 " 

57.5 ■* 

60.0 ” 



GAT N 

19B.85 

197.02 

194.29 

194.00 

195.71 

197.54 

198.90 
199.53 
199.50 
199.16 

198.91 

198.87 
199.05 
199.04 

198.68 

197.88 
106.84 
195.63 
194.39 
193.33 
192.26 

100.68 
iBo . 05 
186,96 


INCIDENCE ANGLE 10. DEGREES 
DA 

2.5 “ 

5.0 " . 

7.5 " 

10.0 « 

12.5 " 

15.0 » 

17.5 ” 

20.0 » 

22.5 « 

25.0 « 

27. 5 " 

30.0 « 

32.5 " 

35.0 " 

37.5 « 

40.0 « 

42.5 '• 

45.0 « 

47.5 *' 

50.0 « 

52.5 « 

55.0 « 

57.5 « 

60.0 « 



GAIN 

194.81 

195.03 

195.52 
196,00 

196.11 

196.53 
196.76 

196.82 
197.07 

197.37 

197.83 
198.14 

198.27 

198.03 
197.50 

196.73 

195.93 

195,25 

194.66 

194.03 
193.20 
19i.71 

190.27 
188.36 



X 



INCIDENCE ANGLE 15, DEGREES 


PA 

2.5 


5.0 

5? 

7.5 

K 

10.0 

«« 

12.5 

n 

15.0 


17.5 

« 

20.0 


22.5 

«* 

. > 

25.0 


27.5 


30.0 


32.5 


35.0 


37.5 


40.0 


42.5 


45.0 


47.5 


50.0 

** 

52.5 


55.0 


57.5 

>> 

60.0 

IT a 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


GAIN 

19 ^ 1.21 

195,31 

195,33 

195.52 

195.52 
13'!. 92 
19'). 37 
13'). 39 
195.03 
195.93 

196.55 
196.75 

196.53 

196.09 
195.30 
19 '). ')B 

19'). 03 

193.55 

193.97 

192.37 

191.53 

190.42 


INCIDENCE ANGLE 20. DEGREES 
RA 

2.5 ” 

5.0 « 

7.5 « • 

lOvO 

12.5 « 

15.0 

17.5 ” 

20.0 « : 

22.5 ” 

25.0 « 

27.5 " 

30.0 « 

32 . 5 •' 

35.0 « 

37.5 " 

40.0 « 

42.5 « 

45.0 " 

47.5 « 

50.0 ” 

52.5 •* 

55.0 « 

57.5 " 

60.0 « 


X 

X 

X 


X 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


-AIN 

'5.34 

14.77 

'3.99 

)2.4l 

12.36 

) 2.66 


193.11 

193.42 

133.58 

lon.n? 

104.09 

193.96 

193.89 

193.48 

193.02 

192.811 

192.74 

192.48 
192.01 


191.70 

191.03 



INCIDENCE 


ANGLE 

25. 

DEGREES 




PA 







GAIN 

2.5 

«• 

0% 




' X 

^0 

197.56 

5.0 

*■ 




X 

!? 

196.78 

7.5 





X 

%0 

195.39 

10.0 

%$ 

0% 




X 

'ft# 

•fit 

193.87 

12.5 

0'* 




X 

ft# 

191.79 

15.0 

%M 

0M 




X 

ft# 

190.02 

17.5 





X 

ft# 

#» 

183.77 

20,0 





X 

ft# 

188.51 

22.5 





X 

ft# 

'.r> 

188.87 

25.0 

•4 




X 


190. p8 

27.5 





X 


191.35 

30.0 





X 

ft# 
.tf H 

191.71 

32.5 





X 

ft# 

192.48 

35.0 

«4 




X 

ft# 

192.53 

37. 5 





X 

ft# 

192.23 

HO.O 

0\ 




X 

ft# 

# V 

192.35 

^42. 5 

V 




X 

ft# 

191.71 

M5.0 





X 

ft# 

191. PO 

47.5 

V 




X 

.ft# 

IQP . 86 

50.0 

5 




X 

ft# 

190.72 

52,5 





X 


190 . 4 n 

55,0 





X 

ft# 

no .pi' 

57.5 

r « 



X 


190 . IP 

60,0 




X 


180.97 

n^cinENCF 


ANGLE 

30. 

DEGREES 




PA 







GAIN 

2.5 





X 


195.53 

5.0 





X 


19/[.0 8 

7.5 


• 





194.00 

10.0 





/\ 


193.11 

12.5 





X 

ft# 

191.81 

15.0 





X 

ft# 

19^.83 

17. 5 





X 

ft# 

190.03 

20.0 





X 

ft# 

180.38 

22.5 





X 

.ft# 

188.54 

25.0 





X 

ft# 

188.79 

27.5 





X 

0.0 

188.86 

30.0 





X 


183,77 

32.5 





X 

.ft# 

189.23 

35.0 





X 

ft# 

4 ^ 

189.39 

37.5 





X 

ft# 

189,41 

40,0- 





X 


189.78 

42.5 





X 

ft# 

vv 

189.61 

45.0 





X 


189.24 

47.5 





X 


189.30 

50.0 





X 


189.32 

52.5 





X 

ft# 

#> 

188.98 

55.0 





X 

ft# 

188.57 

57. 5 





X 

ft# 

188.53 

60,0 





X 

ft# 

188,30 
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TNCIDENCE ANGLE 35. DEGREES 
PA 

2 . 5 » 

5.0 « 

7.5 « 

10.0 « 

12.5 ” 

15.0 « 

17.5 « 

20.0 « 

22.5 » 

25.0 « 

27.5 « 

30.0 « 

32.5 « 

35.0 " 

37.5 « 

40.0 '• 

42.5 " 

45.0 " 

47.5 •* 

50.0 « 

52.5 " 

55.0 « 

57.5 « 

50.0 » 



GAIN 

« 185.53 

« 1P5.44 

« 185.45 

185.55 
« 185.52 

186.95 
187.54 
” 188,15 

5* 188.16 

« 188.58 

” 188,81 

" 133,52 

” 188.75 

” 138.52 

188.46 

•• 188.57 

« 188.30 

” 187.98 

« 188.02 

137.81 
“ 187.41 

« 187.01 

” 186.57 

« 185.29 


INCIDENCE ANGLE 40. DEGREES 


PA 

2.5 " X 

5.0 « X 

7.5 »' ‘ X 

10.0 « X 

12 . 5 X 

15.0 •» , X 

17.5 " X 

20.0 ” X 

22.5 " X 

25.0 « X 

27.5 ” X 

30.0 « X 

32.5 •' X 

35.0 « X 

37.5 » . X 

40.0 « X 

42.5 ” X 

45.0 ■* X 

47.5 " X 

50.0 « X 

52.5 « X 

55.0 - X 

57.5 " X . 

5n n X 


GAIN 

” 185.40 

•' 185.49 

” 185.64 

185.81 
“ 185.93 

« 186.11 

« 186.32 

'• 186.50 

” 185.55 

5' 185.5 8 

185,71 
" 185.50 

“ 185,50 

« 186.49 

185.34 

" 186.30 

« 185.05 

» 185.75 

” 185.50 

« 185.56 

" 185.22 

“ 134.87 

” 134.60 

” 184,20 


A -11 ' 



1 . CT(j micrometers VMVELENGTil 


irCIDEMCE AtriLE 0, DEGREES 
- °A ' 


2.5 X 

5 . 0 « X 

7.5 •« X 

10.0 « • X 

12.5 •' X 

15.0 » X 

17.5 « X 

20.0 « X 

22.5 « X 

25.0 '« X 

27.5 « , X 

30.0 « X ■ 

32.5 “ X 

35.0 « X 

37.5 « X 

40.0 X 

42.5 « X 

45.0 « X 

47.5 « X 

50.0 « X 

52.5 « X 

55.0 « X 

57.5 « ■ X 

60.0 « X 


GAIN 

5' ,^0 

« 198.09 

'• 197.63 

•« 196.69 

" 196.22 
194.81 
" 194.05 

” 192.60 

» 191.17 

« 190 M^O 

” 189.84 

« 189.96 

5! 190,34 

190.68 
" 190.95 

•' 190,98 

" 190.75 

« 190.35 

5- 189,54 

5- 188,57 

5' 187,38 

" 185,77 

" 184.04 

« 182.87 


INCIDENCE ANGLE 
PA 

2.5 X 

5 . 0 ” ' 

7.5 ” 

10.0 « 

12.5 

15.0 « 

17.5 « 

20 . 0 ” 

22.5 « 

25.0 

27.5 '• 

30.0 » 

32 . 5 ” 

35.0 ” 

37.5 " 

40.0 « 

42.5 « 

45.0 '-s 

47.5 « 

50.0 

52.5 « 

55.0 « 

57.5 *' 

60.0 « 


DEGREES 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


GAIN 

■' 00 
19 7.* 21 

" 196.82 

" 196.04 

« 195,64 

" 194.45 

" 193.31 

" 192.55 

" 191.22 

" 190.36 

" 189.37 

" 189.05 

" 189.01 

" 189.26 
" 189.48 

” 189.64 

” 189.64 

” 189.40 

" 189.05 

“ 188.28 

” 187.59 

” 186.36 

5? 184.98 

” 183.75 
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INCIDENCE ANGLE 10. DEGREES 

PA 



2.5 X 

5.0 « 


X 

7.5 •' 


X 

10,0 « 


X 

12.5 


X 

15.0 « 


X 

17.5 •* 


X 

20.0 w 


X 

22.5 '' 


X 

25.0 « 


X 

27.5 » 


X 

30.-0 » 


X 

32.5 « 


X 

35.0 « 


X 

37.5 « 


X 

40,0 


X 

42.5 " 


X 

45.0 ” 


X 

47.5 " 


X 

50,0 


X 

52.5 " 


X 

55.0 « 

X 


57.5 

X 


60.0 « 

X 



GAIN 

« . 00 
» 195.5^ 

« 195.23 

» 19i|,63 

« 19iJ.32 

« 193.^0 

” 192.97 

« 191. 9B 

5' 191.06 

•' 190.33 

" I89v^»^ 

” 133 .'3* 

» 183.53 

« 188.32 

'• 188.21 

” 188.23 

" 188.18 

» 187. 9 

« 187.86 

» 187.13 

« 186.92 

" 186.07 

« 185.35 

" l8iJ.39 


INCIDENCE ANGLE 15. DEGREES 
PA 

2 . 5 X 

5.0 « 

7.5 " ' 

10.0 

12.5 " 

15.0 « 

17.5 « 

20.0 » 

22. 5 « ■ 

25.0 « 

27.5 « 

30.0 « 

32.5 " 

35.0 « 

37.5 « 

40.0 « 

42.5 '■ 

45.0 « 

47.5 " 

50.0 « 

52.5 ” 

55.0 « 

57.5 « 

60.0 


X 

X 

X 

X 

X 

X 

X 




GAIN 

" .00 
« 192.96 

« 192.74 

« 192.35 

•' 192.13 

'• 191.49 

« 191.30 

« 190.58 

« 19"^. 17 

« 189.62 
“ 189,20 
« 188.58 
" 188.53 

" 187.91 

” 187.82 
*« 187.40 

« 187.17 

" 186.96 
” 186.57 

« 185.98 

« 185.69 

« 185.00 

5' 184.60 

« 183.61 
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INCIDENCE ANGLE 20, DEGREES 
PA 

2.5 X 

5.0 « 

7.5 " 

10.0 

12.5 ” 

15.0 « 

17.5 *' 

20.0 “ 

22.5 -J 

25.0 « 

27.5 “ 

30.0 ” 

32.5 “ 

35.0 ” 

37.5 " 

i»0.0 ” 

42.5 ” 

45.0 “ 

47. 5 " 

50.0 « 

52.5 “ 

55.0 « 

57.5 “ 

60-^0 ■ « 



GAIN 

00 

IBP ! 85 
” 188.74 

” 188.6^ 
« 188.47 

” 188.14 

” 188.25 

« 187.79 

" 187.95 

187*60 
« 187.76 

187.25 
” 187.75 

55 186.96 

" 187.35 

“ 186.81 
" 186.46 

55 186.58 

” 185.57 

55 185.90 

55 184.79 

184.52 
55 183.68 

55 183 ;18 


INCIDENCE 

ANGLE 25. DEGREES 




PA 




GAIN 

2.5 

X 



no 

t •' 

5.0 


X 


182.62 

7.5 


X 

MM 

182.70 

10.0 


X 

%M 

183.00 

12.5 

M0 

X 


183.'^^ 

15.0 


X 


183.11 

17.5 


X 


183.66 

20.0 ■ 


X 

%M 

183.49 

22.5 


X 

T > 

134. 33 

25.0 

«« 

X 

mm 

184.18 

27.5 


X 

MM 

. 0 - 

184,00 

30.0 

f\ 

X 

MM 
^ > 

134.57 

32.5 


X 

ft# 

185.56 

35.0 


X 

.ft# 

184.80 

37.5 


X 


185.65 

40.0 


X 

ft# 

185.12 

42.5 

% 

X 

ft# 

184.95 

45.0 


X 

ft# 

✓ ft 

185.38 

47.5 

A* 

X 

ft# 

■#> 

184.20 

50.0 


X 

ft# 

># 

185.30 

52.5 


X 

ft# 

183.81 

55.0 


X 

ft# 

eft 

184.12 

57.5 


X 

ft# 

132.93 

60.0 


X 

ft# 

183.47 
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INCIDENCE ANGLE 
PA 

2.5 X 

5.0 « 

7.5 ” 

10.0 " 

12.5 ” 

15.0 « 

17.5 « 

20.0 « 

22.5 ” 

25.0 « 

27.5 « 

30.0 « 

32.5 ” 

35.0 « 

37.5 " 

40.0 « 

42.5 '* 

45.0 » 

47.5 ” 

50.0 « 

52.5 ” 

55.0 ”■ 

57.5 ” 

60.0 ” 


30, DEGREES 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


GAIN 

« .00 

“ 176.45 

« 176,67 

« 177, 30 

« 177.39 

« 177.77 

178.60 

*= 178,57 

« 179.77 

5! 179.72 

» 180,7^ 

” 180,51 

« 131.78 

181,06 
” 182.22 
” 181.73 

” 131.75 

« 182,40 

5' 181.27 

” 182.83 

- 181.34 

;; T Pp 07 

iSoiSO 

" 182.11 


INCIDENCE ANGLE 
PA 

2.5 X 

5.0 “ 

7.5 " ' 

10.0 

12.5 ” 

15.0 " 

17.5 ” 

20.0 ” 

22.: 5 ‘ » 

25.0 ” 

27.5 ” 

30.0 •'* 

32.5 " 

35.0 « 

37.5 " 

40.0 

42.5 “ 

45.0 « 

47.5 « 

50.0 « 

52.5 ” 

55.0 « 

57 . 5 " 

60.0 ” 


35. DEGREES 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


GAIN 

" ,00 

” 178,83 

” 178.77 

178.68 
178.61 
” 178.41 

« 178,48 

" 178,20 

" 17 8. '31 

178.09 
178.24 
“ 177.90 

178.39 
« 177,78 

" 178,35 

•* 177,92 

- 177,81 

« 178.23 

” 177,29 

178,52 
" 177,27 

” 177,90 

” 176. OP 

178.08 
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INCIDENCE 

ANGLE 40. DEGREES 

PA 


2.5 

< 

5.0 

X 

7.5 

X 

10.0 

X 

12.5 

X 

15.0 

X 

17.5 

X 

20 .0 

X 

22.5 

X 

25.0 

X 

27.5 

X 

30.0 

X 

32.5 

X 

35.0 

X 

37.5 

X 

40.0 

X 

42.5 

X 

45.0 

X 

47.5 

X 

50.0 

X 

52.5 

X 

55.0 

X 

57.5 

X 

60.0 

”STO”5' 

X 


r,ATN 

« . 00 
'* 179 Ji5 

” 179. 3 B 

” 179.25 

” 179.10 

« 178,9G 

" 178.87 

178.62 

” 178.41 

*' 178,18 

" 177.90 

" 177.56 

» 177.38 

" 176.94 

" 176,71 

« 176.34 

" 175.91 

" 175.78 

" 175.01 

« 175.24 

” 174.33 

" 174,25 

» 173.37 

« 173.77 


CROS^LOT 08:50 MAR 06/75 


mo DIMENSION TRANC8), ALISTC18,24) 
110 DATA STAR//’”' /XV, DB/4. 34294/ 
120 CALL OPENFIC’TEDAT* ,1) 

130 DO 1 1=1,18 

140 READC1,>D1,D2,D3/4 

150 PRINT, D1,D2,D3,d4 

160 D02 J=l,12 

170 READ(1,)TRAN 

180 ALISTCI,2«J-1)=TRANC39 

190 ALIST(I,2«J>=TRANC7) 

200 2 CONTINUE 
210 D03 J=1.3 
2 20 3 • RE Af>f-1-/)DUM 
230 1 CONTINUE 
240 CALL CLOSFl 
250 AMAX=0. 

260 AMIN=1.E4 
270 D04 1=1,18 
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280 D04 J=1,24 
290 DUM=AL1ST(I,J) 

300 IF(DUM.GT,AMAX)AMAX=DUM 
310 I FCDUM.lt , AMIN)AMIN=DUM 
320 4 CONTINUE 
330 VO=170 

340 vx=intcdb«amax)+i, 

350 DV=vx-yo 
360 SCALE=DV/50. 

370 PRINT100,VX,VO,SCALE 
380 DPA=2.5 
390 DOll 117=1,2 

400 IADD=9*'CIW-1) 

401 IFCl'/.En,i:)f5RINTl04 

402 IFCIW.EO,2)PRINT105 

403 104 formatC'.53 micrometers WAVELENGTHM 

404 105 FORMATC’1.06 micrometers wavelength’') 

410 DOll J=2,24,2 

420 oA=iJ»DPA 
430 PRINTIO 3 /PA 

440 103 FORMATC/' POINT-AHEAD ANGLE F4 , 0 , ’M I CROP, AD IANS ' 

4415, /T4, ' INC. ' /T 65 , 'GAIN ' 9 

450 DOll 1=1,9 

450 DUM=DB«ALISTCI+IADD, J) 

470 ID=10+INTCCDUM-V0)/SCALE) 

480 IFCID.LT.10)ID=10 
490 AINC=5.“CI-1) 

500 PRINT102, AINC, STAR, ID, X, STAR, DUM 
510 11 CONTINUE 
•520 STOP 

530 100 F0RMATC'MAX,MIN' ,2F8,2/'SCALE = * ,f6.2, '/DIV ') 

540 101 'FORMATC/’ INCIDENCE ANGLE* , f4 ,0,’ DEGREES’) 

550 102 FORMATCF6.1,T10,A,TN,A,T60,A,F6.2) 

560 END 


#RUN CROSPLOT 

CROSf^LOT 

03:52 MAR 

3 . 80000 

.530000 

3.80000 

.530000 

3.80000 

.530000 

3.80000 

.530000 

3.80000 

.530000 

3.80000 

.530000 

3 . 80000 

.530000 

3.80000 

.530000 

3.30000 

.530000 

3.80000 

1.06000 


06, ’75 


.000000 

7.30000 

5.00000 

7,30000 

10.0000 

7.30000 

15.0000 

7.30000 

20,0000 

7,30000 

25.0000 

7.30000 

30.0000 

7.30000 

35.0000 

7.30000 

40.0000 

7.30000 

.000000 

7,30000 
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3.80000 

1.06000 

5.00000 

7.30000 

3.80000 

1.06000 

10.0000 

7.30000 

3.80000 

1.06000 

15.0000 

7.30000 

3.80000 

1.06000 

20.0000 

7.30000 

3.80000 

1.06000 

25,0000 

7.30000 

3.80000 

1.06000 

30.0000 

7.30000 

3.80000 

1,06000 

35.0000 

7.30000 

3 . 80000 

1.06000 

40.0000 

7.30000 


“1AX,MIN 201.00 170.00 

SCALE= . 62/di V 

. 53 MICROMETfeRS WAVELENGTH 


>^01 NT-AHEAD ANGLE 
INC . 

,0 

5.0 « 

10.0 - 

15.0 « 

20 . 0 ” 

25.0 " 

30.0 “ 

35.0 " 

40.0 


5.MICRORADIANG 




X 


X 

X 


X 


GAIN 

197.27 

197.02 

195.03 
192.01 

194.77 

196.78 
194.98 
186.44 

IB5.49 


POINT-AHEAD ANGLE 10. MICRORADIANS 
INC. 

.0 « 

5.0 •' 

10.0 " 

15.0 " 

20.0 

25.0 * 

30.0 ” 

35.0 " • 

40.0 


X 


X 


X 


X 

X 


X 

X 

X 


GAIN 
101.61 
194.00 
196, on 
195.31 

193.43 

193,87 

193.11 

136.66 

135,81 


■^OI NT-ANEAD ANGLE 15 ,M I CRORAD I ANS 
INC. 

.0 ” 

5.0 » 

10.0 ” 

15.0 w 

20 . 0 •« 

25.0 « 

30.0 « 

35.0 « , X 

40.0 « ' X 


X 

X 

X 


X 

X 

X 

X 


GAIN 

- 198.83 

” 197,54 

« 196.53 

195.62 
” 192.36 

•' 190.02 

“ 190.83 

« 186.95 

5 ' 186.11 


POirjT-AflEAD ANGLE 20 .MI CRORAD TANS 
INC. 

.0 « 

5.0 « 


GAIN 

X« 200.48 

X « 199.53 


A-.18; 



X 


10.0 5 « 

15.0 « 
20.0 « 
25.0 « 
30.0 « 

35.0 « 

110.0 « 



POINT-AHEAD ANGLE 25.MICRORADIANS 
INC. 

.0 « 

5.0 » 

10.05' 

15.0 !« 

20.0 5 { 

25.0 5 : 

30.0 55 

35.0 55 

40 .0 5S X 


X 


X 

X 



X 


X 

X 


POINT-AHEAD ANGLE 30 .III CRORADI ANS 
INC. 

. 0 5 : 

5.0 5 ! 

5< - 



15.0 

20.0 
25.0 
30,0 

35.0 

40.0 


X 

X 

X 

X 

X 


«ur 


POI NT-AHEAD . ANGLE 35 .MI CRORADI ANS 
INC. 

.0 
5.0 
10.0 
15.0 
20.0 
25.0 
30.0 

35.0 

40.0 


X 

X 


X 


X 


X 


X 


POINT-AHEAD ANGLE 40 .M I CRORADI ANS 
IRC. 

. 0 5C 

5.0 5 ! 

10.0 5S 

15.0 ” 

20.0 *s 

25.0 5 : 

30.0 5 { 

35.0 5 : 

40.0 5{ 


X 
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196,82 
194,92 
192 . 81 
188,51 
189,38 
188,16 
186,50 


GAIN 

19 ’. 30 
199,16 

197.37 
194,39 
193.11 

190.38 
188.79 
188,68 
186.6B 


GAIN 

200.01 

198.87 

198»-14 

195.93 

193.58 

^Si-71 

188.77 

188.52 

186,60 


GAIN 

200,33 

109,04 

198.03 

196.75 

19^.09 

192.53 

189.39 

188.62 

186.49 


GAIN 

198.08 
197.88 
196.73 

106.09 
193,80 
192.35 
189.78 
188.57 
186.30 


X 



f’OTNT-AHEAD ANGLE i45.MT CRORAOI ANS 
INC. 

.0 

5.0 
10.0 
15.0 
20.0 
25.0 
30.0 

35.0 

40.0 ” X 

POINT-AHEAD ANGLE 50 ,MI CRORADI ANS 
INC. 

.0 

5.0 

10.0 
15.0 
20 . 0 
25.0 
30.0 

35.0 

40 .0 

"OINT-AHEAD ANGLE 55 . M I CROR ADI ANS 
INC. 

.0 
5.0 

10.0 
15.0 
20.0 
25.0 
30.0 

35.0 

40.0 

POINT-AHEAD ANGLE 60 .MI CRORADI ANS 
INC. 

.0 " X 

5.0 

10.0 
15.0 
20,0 
25.0 
30.0 

35.0 

40.0 


X 


X 


X 


X 


X 


X 


GAI 

N 

195, 

33 

195. 

63 

195. 

25 

194. 

48 

193. 

0 2 

191. 

00 

189. 

24 

187. 

98 

185. 

76 

GAI 

N 

193. 

.99 

193. 

33 

194, 

,03 

193. 

,55 

192, 

,74 

190, 

,72 

189. 

.32 

187, 

,81 

185. 

. 56 

GAIN 

190, 

.38 

190 

.68 

191 

.71 

192 

.37 

192 

.01 

190 

.0 4 

183 

.57 

187 

.01 

184 

.87 

GA 

TN 

185 

.13 

186 

.96 

188 

,36 

190 

.42 

191 

,03 

189 

.97 

188 

.30 

186 

.29 

134 

,20 


.06 MICROMETERS WAVELENGTH 

^01 NT-AHEAD ANGLE 5, MICRORADIANS 
INC. 

. 0 « 

5.0 , « 


GAIN 

108,09 

197.21 


A-20' 



X 


10.0 '! 
15.0 « 
20.0 « 
25.0 « 
30.0 « 
35.0 « 

Ho . 0 " 


X 

X 


X 


X 


X 


195.54 
« 192.96 
« 188.85 
« 182.62 
« 176.45 

« 178.83 

« 179.45 


POIfIT-AHEAD ANGLE 
INC, 

.0 « 

5.0 « 

10.0 « 

15.0 « 

20.0 « 

25.0 « 

30.0 « 

35.0 « 

40.0 5« 


10. MICRORAD IANS 


X 

X 

X 

X 



GAIN 

« 196,69 

» 196,04 

« 194.63 

« 192.35 

« 188.60 

« 183.00 

177.30 
178,68 
» 179.25 


'’OINT-/\HEAD ANGLE 15 , M T CROP ADI ANS 


INC. 

.0 •” 

5.0 « 

m.o ” 

15.0 « 

20.0 5! 

25.0 « 

30 . 0 

35.0 ” 

40.0 ” 


X 


X 


X 


X 


X 


X 

X 


GAIN 

194.81 

194.45 

193.40 
191.49 
188.14 
183,11 
177,77 

178.41 
178.96 


^0 ! N T- Ai \ E AD' AN G LE 2 0 . : 1 1 C RO R AD I AN5 


INC. 

.0 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

40.0 


X 


X 

X 

X 


X- 


X 

X 


X 


X 


GAIN 

192,60 

192.55 

191.98 

190.58 

187.79 

183.49 

178,57 

178,20 

178.62 


point-ahead ANGLE 25. MICRORAD I /^NS 
me... ■ - , . . _ . : 


,0 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

40.0 


X 

X 

X 

X 


X 




GAIN 


190.40 

190.36 


190.33 

189.62 


187,60 

184,18 


179.72 

178.09 

178.18 


X 

X 


X 


X 



POINT- AHEAD ANGLE 

30.M1CRORAOIANS 


INC. 




.0 

n 


X 

5.0 

«• 


X 

10.0 



X 

15.0 



X 

20.0 

55 


X 

25.0 


X 


30.9 

55 

X 


35.0 

55 

X 


40.0 

55 

X 



GAIN 

« 189.96 

" 189,05 

« 188^84 

« 188.58 

« 187.25 

« 184.57 

5« 180.51 

« 177.90 

« 177.56 


POINT-AHEAD 

INC. 

.0 » 

5.0 ” 

10.0 « 
15.0 » 
20 . 0 « 
25.0 " 
30.0 « 

35.0 » 

40.0 « 


ANGL E 35 . M I C ROR AD I AMS 


X 

X 

X 

X 


X 

X 

X 

X 

X 


point-ahead angle 4o.microradians 


INC. 




.0 ' 




5.0 • 




10.0 J 



X 

15.0 • 



X 

20.0 ' 



X 

25.0 « ' 


X 

30.0 5 

35.0 5 

X 

X 


40.0 -5 

X 




X 


X 


POINT -AHEAD 
INC. 

.0 « 

5.0 « 

10 . 0 « 
15.0 « 
20.0 « 
25.0 « 
30.0 « 

35.0 « 

40 . 0 ” 


ANGLE 45.MICRORADIANS 


X 


X 


X 


X 


X 

X 

X 


X 

X 


POINT- AHEAD 
INC. 

.0 

5.0 ” 

10.0 " 


ANGLE 50 .MICRORADIANS 


X 

X 


GAIN 

” 190.68 
” 189.26 
- 188.32 
“ 137.91 

” 186.96 

5^ 184.80 

« 181.06 

'* 177.78 

176.94 


GAIN 

« 190.98 
“ 189.64 
« 188.23 
« 187.49 
” 186.81 
« 185.12 
« 181.73 
« 177.92 
" 176.34 


GAIN 

” 190.35 

" 189.40 

” 187.94 

« 186.96 
« 186.58 
« 185.38 

« 182.40 

« 178.23 

” 175.78 


GAIN 

” 188.57 

5 ^ 188.28 

" IS7.I3 
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X 


X 


15.0 

20.0 

25.-^ 

30,0 

35.0 

40.0 



X 


X 

X 

X 


'^OIMT-AHtiAD ANGLE 55. MICRORADIANS 


I NC . 



.0 •' 


X 

5.0 » 



10.0 « 


X 

15.0 « 


X 

2 0 , 0 "• 


X 

25.0 


X 

30.0 « 


X 

35.0 


X 

40 . 0 ” 

X 


oOINT-AHEAD 

ANGLE 

6o,microradians 

INC. 



.0 « 


X 

5.0 " 


X 

10.0 « 


X 

15.0 " 


X 

20,0 •• 


X 

25.0 


X 

30.0 


X 

35.0 " 


X 

40.0 " 

X 



5'STOf’" 


« 185,98 

” 185.90 

« lB5.'3fl 
182,83 
» 178.52 

« 175.24 


GAIN 

» 185,77 

” 186 , 3 ^ 

« 186.07 

" 185.00 

« 184.52 

« 184.12 

5 « 182.07 

“ 177.90 

" 174.25 


GATH 

« 182.87 
» 183,75 
« 184.39 
” 183.61 
« 183,18 
5= 183.47 
« 182.11 
5 > 178.08 
5' 173.77 


A*-23 



APPENDIX B 


LAGEOS SIGNATURE DATA 


T.O INTRODUCTION 

The temporal response of the model LAGEOS target varies with point-ahead 
angle, wavelength, and illuminating pulse width. The data and the program 
used to generate this data is presented in this appendix. The gain data used 
in this analysis was generated by the program in Appendix A. 


2.0 SIGNATURE ANALYSIS PROGRAM 

The program used to compute the temporal response of the LAGEOS target 
is presented in the following pages. The program was written in Fortran IV 
for the MDAC Sigma 9 time-sharing computer system. 


couvv 09:16 MAR 06,'75 


100 DIMENSION AM CO , 2 4 0 , AVC9 ^ 0^0 > TR ANf 80 
110 C0M*1ON AM,AV 

120 DATA D2LT,C0M1/1.87794E-11,2.5'^B/ 

130 wAVf-=.53 
140 22 OUTPUT v/A'/E 
150 READ/7AVE 

160 call OPENFIC TEDAT' ^10 

170 IFO/AVE.LT.l. OOOTO 25 

130 0021 1=1,144 

190 READCl, 9 d1/d2,d3/d4 

200 21 CONTINUE 

210 25 CONTINUE 

220 DO 20 1=1,9 

230 READC1,>D1,D2,D3,d4 

240 PRINT,D1,D2,D3,d4 

250 DO 23 J=l,12 

260 READC1,0TRAN 

270 AMCl ,2-'.J-1') = TRANC3') 

280 AMCT;,2-N0=TRANC7) 

29 0 AV C T , 2 J -1 ■) = T R AN ( 4 0 

3 0 0 A V C I , 2 ‘N J ) = T R AN C 8 ) 

310 23 CONTINUE 

320 READC1,)DUM 
330 READCl, ODU*! 

340 RE AD Cl, ODUM 


B-1 



350 20 CONTIN'.J'; 

3 30 CALL CLOSFl 

370 1 ouTPijT"’UL0E '/idt;i' 

3o0 RZ/VD>TPW ■ - 

300 2 OUTPUT' TMIM,TMAX/TSTEP' 

400 OSAJ^TM IN,TMAX,T5TEP 

'HP D-0 5 1 = 1,12 

42 0 II-O^vT 

430 pa=-5.=H 

440 T = TMr-J 

450 4 TN=-.5”OELT 

460 SUM=0 

470 5UMVAf’=0. 

4 30 SU'11=0. 

490 V‘0l=O, 

500 DO 3 M=1,76 

510 TM=TU+DELT 

520 0 A:' 1=ACOS Cl . -TN”COUl ) 

530 I P(OAM .r,T .0 . 7854001=0, ;GOT010 
5 40 CALL eSTVALC Dl , VAR , 0 A.' 1 , MO 
350 S JMl=OUMl->-Dl 
360 VS1=VS1+VAR 

570 D'j;i = EX P (-2 . " C CT-TM ) / TP’ /O ””2 ) 

530 S JMV AR=S UnVAR +'/ AR-'DUfl" " 2 

590 10 COMTIUUE 

600 3 5U:i = SUM+Dl"DUM 

610 I F C CT- TH I r 1 ) . GT . TS TEp / 2 , >G 0T06 

■5 2 0 G A V E = 1 0 , A L OG 1 G C S U<M 1 0 

630 S r ! R = 2 . " G ;\y E - 1 0 . " A L OG 1 0 C '/ S 1 0 

6 40 PRI.’iT 101, PA, GAVE, SNR 

650 101 FORMATC/F7* 0 , f8 . 2, F7. 2) 

660 J CONTINUE 

670 GAVS=10 . ”ALOGl0CSU:n ’ 

630 SnR = 25''GAVE-10 . "ALOGlOCSUMVAPsO 

690 PRINT 10 0,T, GAVE, SNR 

700 100 FOR‘1ATC9pF7.2,pf3.2,f7.2) 

710 T=T+TSTEP 

720 IFCT. LE.TMAX)G0T04 

730 5 COfJTINUE 

740 OUTPUT* rlE' I VALUES OR MEU RUN ’ 

750 IN=0 

760 P>EAD,I‘1 
770 IFTTM.[rn^l)r,r)X02 
730 IFnM.E0.2)G0T0l 
790 I F C n . EO . 3 )G0T022 
300 STOP 
310 END 

820 SUBROUTINE ESTVALCO.NU, VAR, T IE, JpO 
330 DIMENSION A*1C9,24)^ AVC9,24) 

340 COMMON AM,AV 
350 GNU=0. 

360 VAR=G. 

S70 XI=THE”11.4592 
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380 IL=IMTCXI)+1 

390 ,2M)RETURM 

OOG ifcil.ot.ogreturn 

910 G2=0. 

920 V2=n, 

930 G1=AMCIL,.JR9 
940 V1=AVC IL^ J”') 

050 IFriL.En.QOGOTOl 
960 G2 = A‘iriL+l, J'^) 

970 V2= AyCT L+l, J'^') 

980 1 G=G1+CG2-G1)"CXT+1-IL) 
990 V=Vl+CV2-Vl)”CX I+l-IL) 
1000 G.9U=EXRCG) 

1010 VAR=FX'‘>C2”V') 

1020 RET'JR'J 
1030 E^|0 


3.0 SIGNATURE DATA 

The following pages present the temporal response characteristics of a 
LAGEOS type target, when illuminated by a 200 ps pulse at 0.53 pm and at 
1.06 pm wavelengths. 
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FIGURE B-2 
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FIGURE B-4 
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FIGURE B-7 
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FIGURE B-8 
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FIGURE B-9 
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FIGURE B-14 
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FIGURE B-19 


B-22 




TiM£ Cus') 


B-23 


FIGURE B-20 





































APPENDIX G 


LAGEOS SIGNATURE CENTROID CALCULATION 


1.0 INTRODUCTION 

The nonsymmetrical shape of the return pulse (Appendix B) from LAGEOS 
results in the pulse centroid varying considerably with conditions and pulse 
width. This program used to calculate the centroid and the resulting data 
are included in this appendix. The data used for these calculations were 
gerfe rated by the programs presented in Appendix A. 

2.0 BIAS ESTIMATION PROGRAM 

The program used to compute the centroid of the pulse was written in 
Fortran IV for MDAC's Sigma 9 time-sharing computer system, and was used in 
the on-line mode. The program is listed in the following pages. 


& I A S i S T 10:30 MAR OS, ' 75 


100 DIMENSION AMO, 24) ,AV(S, 24), TlUNC3),A(7S) 

110 COMMON AM,AV 

120 DATA DELT,C0N1/1.37794£-11,2.5E3/ 

1^0 DATA Al, A2,A3,A4/.2 78595 , .23 03 89 ,.0009 72 , .078108/,RTPI/1 . 77245 / 

1 40 inF (X ) = 1 . - 1 . / ( 1 .+ X* ( A 1 + X* ( A2-(- X* ( A3+X* A4 ) ) ) ) ** 4 

150 WAVE:. 53 

ISO 22 OUTPUT WAVE 

1 70 R 2A D, WA VE 

130 CALL OPENFK'TEDAT 1) 

ISO IFCWAVE.LT.l .) G0T025 

200 D021 I r 1,1 44 

210 READCl ,) D1 ,D2, D5,D4 

220 21 CONTINUE 

250 25 CONTINUE 

240 DO 20 1=1,9 

2 50 itEAOC 1 , ) D1 , D2,D3, D4 
2 SO PRINT, 01 ,D2, D3,D4 
270 DO 23 J=1 , 12 

280 HEAD(1,)TRAN 

290 AM(I,2*J-l)=TRAN(3) 

50 0 AH (I ,2*J) = TRAN(7) 

310 AV(I,2*J-l)=TRAiN(4) 
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i20 AVd ,2*J)-TftAN(8) 

550 25 CONTINUZ 

540 READ(l,)DUM 

550 READ(1,)DUM 

5S0 nZAD(l,)DUM 

5 70 20 CONTI 

580 CALL CLOSFl 

5S0 1 OUTPUT 'pulse WIDTH' 

400 READ, TP W 
410 D05I=1,12 
420 11-2*1 
450 PA:5.*I 
440 4 TN=-.5*D£LT 
450 DO 7 J=1 , 76 
460 7 A(J)=0. 

470 DO 10 N = l,76 

430 TN:TN+DELT 

490 GAMrACOSU .-TN*C0N1 ) 

500 CALL ESTVALCDl , VAR,GAM,1I) 

510 A(M) = D1 
520 10 CONTINUE 
550 SUMl-0. 

540 SUM2-0. 

550 DO 5 N:l,7S 

560 sum =5Um+D£LT*CN- 1)*A(N) 

570 5 SUM2=SUM2+ACN) 

580 SDLT-SU.V11 /SUM2 
590 SUMlrO. 

SO 0 DOS N- 1 , 76 
SIO NN = N 

S2 0 SUMl :SUm + ACN) 

S5 0 I P C S UM 1 . 3T .S UM2 /2 , ) GOT 09 
S4G 3 CONTINUE 

650 9 SDST-(FLOATC NN)- .5)*DELT 
SSO EHHrSDST 

670 IF ( TPW . GT .0 . 5E-9 ) ERR rSDLT 
680 11 SUMFrO. 

690 SUMDF-0. 

70G DO 12 N=l,76 
7 1 0 DUM r ( DELT* ( N - 1 ) - ERR ) /T? W 
720 DUi'l 1 : £ R F ( AB S ( DU K ) ) 

750 DUM2 : SI GN ( DUM 1 , DUM) 

74 0 S U M F =S UM F+ A ( N ) * DUM2 

750 12 SUMDF = SUKDF+ACN)*£XP("DU1*1**2) 

760 D EL ER R : S U MF * P W *R TP I / SU MD F / 2 . 

7 70 ERR-ERR+DELERR 

780 lF(ABS(DEL£Ra/£KR) .GT. 1 .E-4) GOTOll 
79 0 PR I N T 1 00 , PA , ER R , SD ST , S DLT 
800 100 FaRMAT(F7.0,9PF3.5,2F3.5> 

^10 5 CONTINUE 

82 0 OUTPUT 'new VALUES OR NEW RUN' 

"750 lf»):0 
840 READ.IM 
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S50 IF(IM.EQ.2)QOTOl 
IF(IM.Eq.3)G0T022 
8 70 STOP 
<»30 iND 

R50 SUBROUTI NE ESTVALC GNU, VAR, THE, JP) 
SOO DIMENSION AMC9,24) , AV(9 ,24) 

910 COMMON AM,AV 
920 GNUrO. 

93 0 VARrO, 

940 X1:TH£*1 1 .45’92 
95 0 IL=INT(XI) + 1 
9 SO lF(JP.GT.2 4)RETU,iN 
970 IFdL.GT .9) RETURN 
980 G2=0. 

99 0 \£=0. 

1000 G1=AM(1L,JP) 

1010 V1 = AV(IL,JP) 

1020 lF(IL.Ea.9)G0T01 
1050 G2: AM(IL+1 , JP) 

1040 V2:AV(IL-M , JP) 

10 50 1 G:G1+(G2-G1 )’k(Xl+l-IL) 

1 060 V : V 1 4- ( V2 - V 1 ) * ( X I + 1 - I L) 

1070 GNU=£XP(G) 
lOSO VAR-EXP<2*V) 

1090 RETURN 
1100 END 


3.0 DATA, 

The data resulting from running this program is presented in the follow- 
ing pages. The printout is in two distinctive parts. The first part of the 
printout contains 9 lines of data, read from the data file containing the 
LAGEOS cube-corner gain data, and is printed to ensure a correct data input 
has accured. Following entry of the (half) pulsewidth, the next twelve lines 
of data list the pointahead angle in microradians, the computed pulse centroid 
(nanoseconds) for the chosen pulsewidth, and the short and long pulsewidth 
approximations for the centroid (nanoseconds). 
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WAVE = .530000 


3.80000 

5.30000 

3.30000 

3. 30000 


.530000 

.550000 

.530000 

.530000 


.000000 

5.00000 

10.0000 

15.0000 


7.30000 

7.30000 

7.50000 

7.50000 
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5.50000 

.550000 

20.0000 

7,50 000 

5.50000 

.550000 

25 .0000 

7.50 000 

5.50000 

.550000 

50.0000 

7.50000 

5.50000 

.550000 

55.0000 

7.50000 

5.30000 .550000 

PULSii WIDTH 

40 .0000 

7.50000 

7.1E-9 

5 . 

.553 

.547 .555 


10 . 

.2 79 

.291 .5 08 


15. 

.180 

.178 .255 


20. 

. 1 52 

.141 .258 


25. 

. 1S6 

. 160 . 24 7 


50. 

. 1 55 

.1 60 .256 


55. 

. 170 

.160 .258 


40. 

.198 

.197 .262 


45. 

.222 

.216 .285 


50. 

.255 

.254 .511 


55 . 

.284 

.291 .555 


SO . 

.556 

.3 4 7 .5 72 


NZV VALUE 

S OR NEW 

RUN 


?2 

PULSE WIDTH 



?2 .i-9 

5 . 

.555 

.54 7 .555 


10 . 

.50 7 

.291 .508 


15.' 

.254 

.173 .255 


20. 

.256 

.141 .253 


25. 

.246 

.160 .247 


50. 

.255 

.160 .256 


55. 

.257 

.ISO . 258 


40. 

.2 51 

.197 .2 62 


45. 

.284 

.216 .285 


50. 

.510 

.254 .511 


55 . 

.555 

.291 .555 


SO. 

.5 71 

.547 .572 


VALUE 

S OR NEW 

RUN 


?5 

WAVE r .550000 



?1.05 

5 .50000 

1 .0 6000 

.000000 

7.30000 

5 .300 00 

1 .06000 

5 .00000 

7 .30000 

5 .30000 

1.06000 

10.0000 

7 .30000 

5 .30000 

1 .06000 

15.0000 

7.30000 

5 .30000 

1 .06000 

20.0000 

7 .50000 

5 .80000 

1 .06000 

25.0000 

7.30000 

5 .300 00 

1 .0 60 00 

50.0000 

7 .500 00 

5 .30000 

1 .0 6000 

55.0000 

7 .30000 

5 .80000 1 .0 SO 00 

PULSE WIDTH 

40.0000 

7.30000 

?.l £-9 

'5. 

.088 

.085 .125 


10 . 

.098 

.085 .158 
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15. 

.110 

.105 

20 . 

. 150 

. 122 

25. 

.161 

. 1 60 

50 . 

. 137 

. 19 7 

55. 

.196 

.197 

40. 

.204 

.216 

45 . 

.220 

.255 

50 . 

.248 

.25 4 

55. 

.258 

.2 72 

60 . 

.504 

.510 

NiW VALUES 

OR NEW 

RUN 

?2 

PULS? WIDTH 


72.E-9 

5. 

.124 

.085 

10 . 

.157 

.085 

15. 

.154 

.105 

20 . 

.179 

. 122 

25. 

.210 

.1 SO 

50 . 

.254 

.19 7 

55. 

.241 

.197 

40 . 

.246 

.21 6 

45. 

.258 

.255 

50 . 

.23 0 

.25 4 

55. 

.290 

.2 72 

SO . 

.52 6 

.510 

rJEW VALUES 

OR NEW 

RUN 

? 

. t 


.155 
.179 
.211 
.255 
.241 
.246 
.258 
.28 0 
.291 
.526 


.125 
.158 
.155 
. 1 79 
.21 1 
.255 
.241 
.246 
.258 
.280 
.291 
.52 6 
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APPENDIX D 


LINK MARGIN ANALYSIS 


1.0 INTRODUCTION 

The link margin analysis summaries and the program used to generate them 
are presented in this appendix. 

2.0 LINK MARGIN ANALYSIS PROGRAM 

The program used to compute the link margins for the various links and 
terminals was written in Fortran l]I , for MDAC's Sigma 9 time-sharing computer 
system, and may be used either on-line or in the batch mode. The program 
is listed in the following pages. 
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I;4THGEn. LISTC24n) . LMJ1;C36) • 



DlCiU^SlOt^ .VALC30 0 



120 

DATA '’I ,OLAMCK,SPnL/3.1' 

'11592653, 

6 , 625pE 

n 'JO 
-> j 

DATA LIST/ ’TRANSniTTER ' 

■’OWER 

DBW 

lAOS 

‘TRAiiSMITTCP. LOSSES 

DO 


150 s 

•TCAMSMIT ANTEfRJA GAIN 

DO 

/ 

160s 

'DOI.iTING LOSS 

DB 


170 c 

‘FREE SPACE LOSS 

OB 

f 

9 

i3o<; 

' ATnoSPIIERI C LOSSCDO'./iMO 

DB 

I 

/ 

1.10^, 

* T ARGET GAUL PRODUCT 

DO 

P 

j KT 

’ ATf'.OSP. IE RI C LOSS CU”^ ) 

DB 

1 

P 

210G 

SPACE LOSS 

DB 

P 

221;: 

*RECEr'/E ANTENNA GAIN 

DB 

t 

/ 

2 •Jno 

‘nECEI VER LOSSES 

DD 

1 

P 

2A 

* CE I VE D S I GN AL PO'. /ER 

DBW 

f 

2 30o 

’ TRANSHIT PULSE RATE 

DB/5EG 

t 

) ^ A 

‘ENERGY ^cR PULSE 

D3J 

t 

2?OG 

‘ENERGY ®ER PHOTON 

D3U 

f 

P 

2305 

' RECHI VED PilOTONS/^ULSE 

DB 

1 

P 

2105 

‘•TUANTUM EFFICIENCY 

DB 

t 

/ 

0 r> p 

'RECEiyED P-E/r’ULSE 

D3 

« 

/ 

3106 

‘ d ACnGROUHD rad 1 AN CE 

DBW 

P 

3205 

‘ RE C E 1 V c R FO 7 C S T E R RAD. ) 

D3 

1 

P 

390- 

‘OPTICAL FILTER B AND'. /I DTHDS 

% 

P 

34c 5 

‘RECEIVE ANTENNA AREA 

DBM2 

f 

P 

3305 

‘RECEIVER LOSSES 

OB 

P 

3b05 

‘RECEIVED BACK. POWER 

DBW 

t 

P 


3A,2,e97925E8/ 

t 



37ncMN£RGv oi£n t’loTon dbj S 

3 OOC 'RilCCIVED P'lOTOislS/GEC DB/SEC 

39 ns ’ HJANTUM EFFICIENCY D!3 

ilOOf; ‘BACKGROUND P-E/SEC Dii/SEC \ 

410(9 'RECEr/ER GATE UIDT.U OEJ~GEC 

429G 'l^ACKGROUND P-E/GATE DB '/ 

■ ••}3'^ data LTNK/’UATRS to LATSCOASGIVE) ■ 

44nc‘ lATRS TO LATSCACTTVE9 ’ 

4590, ‘EATS TO 'lATRS * 

4G9S , '. IATRS to GROUJIDCp ASS I VE9 » 

' 1700 'GUUTTLE TO GROUND ' 

4 80.9 'GROUND BEACON TO SNUTTLE'/ 

490 DATA SIGMA, PN0TN,1MAR/. 13333E-9, .2E-9/9.^ 

50 c DATA APT,ALT, JNDT;I,^OINT,v;AVE/.O01, . 8, 5 . E-6 , 1. F-6 , , 5?E-6/ 

5190 , RANGE, A ri tD, DTAR, ATMU/ 3. F7, 1. , 20 8 . 5, 1. / 

5 20 5 , 9 RI D I A, 5 E CDI A, ALR , ” R AT E, 0 E / 1 . , , 1 9 , . 4 , 1 . , , 2 5 ' 

5 30 s , BRAD, FOV, FI LTER, RGU, 5PG, LTROL, NEXT/ . 017^. 5E-4 , 5 . / 2 . E-9, 1 . , 1 , 1/ 
54 OR ,PND, iVFA, B LOO^/1. E- 6 , 1. E- 8 , 10 . / 

559 NANELIGT 

560 DBCXO-lO.i'ALOGlOCX) 

57'^ CALL OPENFIC ‘OARAM’,!') 

580 GO TO 21 

590 22 COiTTINUE 

50C I FCNEXT.E0.1)G0T02n 

610 CALL CLOSFl 

5.20 STOP 

6 30 20 CONTI i iUE . . - 

540C COMPUTE LINK PARAMETERS 

650 D'jn=S TGr-1A ’ 

550 IFCO!Jfl,GT.n, 2 "PUDTN)DUM= 0 , 2 ”Pl'/DTN 
u7^ ACO;j = CPI”OUf1/P!'/DT:i)F"2 

680 i.nux=0 

690 IFCLTROL,EO,2)INDEX=1 
7^0 I'-CLTROL.E0.3)IMDEX=1 
719 I FOLTROL. E9,6)INDEX=2 

729 D9T=DLiCA9T) 

730 0LT=DBCALT) 

740 DU,M=32./Bv/DTN”"2 
75 O GTX=DBCDUM) 

76 0 D L9 T G =- 8 . 6 85 8 9 •' C 2 P 0 1 )N T / B ' I D Ti O ” •'* 2 
779 DUM=\/AVE/4. /PI / range 
730 FSL=2.-DBCDUM) 

790 DAT?1D=D3rAT;iO) 

300 DATMU=DBCATMU) 

310 R A A = P I « c ”p. I D I A** « 2- S E C D I '■ 2 ')/ 4 . 

320 TRG=Pw^A”4.”PI/17AVE""2 
839 GTRG=UaCTRG) 

349 DLR=DBCALR) 

350 DPRX=DPT+DLT+GTX+DLPTG+2”FSL+DAT^1D+DTAR+OATMU+GTRG+nLP 
360 IFCINDEX. EO.nOGOTOlO 

870 nnRX=DPRX-'^SL-DATrOJ-nTAR 

830 I FCINDEX. E0.2')DPRX=DPRX+DAT.'1U-DAT;in 

399 i.n COiJTIN'JE 



9on dprate=dbcprate') . 

910 DE'’=DPRX-DPRATE 
920 EPP='’LAi--ICK«3PDL/'7AVE 
920 OEP”=DB(E'’P) 

'skn DRP'’=D£'’-nEPP 
950 r3'VE = DiiCOE) 

960 DR'’EP = DRPP+DOE 

970 RP€»«Mi"’KORPEP,A10, 9 

9SO DBRAD=DliCaRAD9 

990 DUri=PI*vCFOV/2,>”"2 

1000 DF0V=Dt3CDUfO 

1010 DFILT=D3CFILTER). 

1020 DRAA=DliCRAA) 

1030 DR3P=UURAD+DFOV+DFI LT+DRAA+DLR 
10^0 DRP5=DRL}P-DEPP 
1050 ■)3PES=DR'’S + DQE 
lOrO DRG,/ = D3CRG'./9 
1070 OiiPEG=nuPES+nRG'/ 

10 30 i3^EG = in.””CPBPEG/10.) 

1090 I FC I MAR. E'P . 1)G0T0250 

1110 ms I G= Cl . +S O RT Cl , +2 "U P EG” ACON99 / AGON 

1110 TFCIMAR. En.2)ROSIG=5”BLOOP”ri.+SORTCl. + , it”3'’EG"PRATE/3LOOP99/aDATE 
1115 rFCl."iAR,E0.39CALL T.1P.ESHCPFA,PND,3 oEG,ROSIG9 
1120 250 COrJTIfiuE 

1130 drosig=dbCrgsig) 

1140 AMARG=RPEP/R05IG 
1150 DMARG=DBCAMARG) 

1160c PUT OUTPUT PARAMETERS INTO VALUE LIST 
1170 valC1)=dpt 
1180 VALC2)=DLT 
1110 valC 39 =gtx 

1200 '/ALC4)=DLPTG 
1210 VALC5)=FSL 
1220 VALC6) = DATfin 

1230 valC 7 )=dtar 
12 40 VALC8)=DATf1U 

1250 valC 9 )=fsl 

1260 VALCIOO-GTRG 
1270 VALC11)=DLR 
1280 VALC129=DPRX 
1290 valC 139 =dprate 
1300 '/alC14)=dep 

1319 VALC15)=DEPP 

1329 VALC16)=DRPP 

1330 valC179=doe 
1340 VALC18)=DRPEP 
1350 valC 19 )=dbrad 
1360 VALC20)=DFOV 
1370 VALC219=DFILT 
1380 VALC22)=DRAA 
1390 valC 23 )=dlr 
1400 VALC24)=DRBP 
1410 VALC25)=DEPP 
1420 VALC26)=DRPS 
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1430 VAL(27)=D0E 
1440 VAL(28>=D3PES 
1450 VALC290=DRGV/ 

14 60 VAL(3O0=DiiPEG 
1470C P0If4T LINK SUflMARY 
14 30 J.JJ=6”LTR0L 


1490 JJ=JiJJ-5 

1500 PRINT200,CLINKC'J)/>J=JvJ,J-JJ) 

1510 200 FOPMATC//*Sa4) 

1520 PRINTIOO, WAVE, RANGE /PWDTH . , 

1530 100 FORMATC ’ Cl RCUI T MARGIN SUMMARY ' A* OPT I CAL WAVELEiiGTU ' , 
1540GGPF6. 2, ' MI CROMETERS, RANGE ’ 

1550G ‘r.|>LSEV/IDTH=%9PF4,2, ' NANOSECONDS »//T6, ‘PARAMETER’, 

15C)0G T 32, ’UNITS V/T44, ’VALUE ') 

157^ DO 3 1=1,30 

1530 I FCTNDuX. EO.OOGOTO ll 

1590 I FCCCI .LT.79.OR. Cl .GT. 99 9. AND. INDEX.E0,19G0T011 
l6oo ifCCCI .LT.59.0P- «CI .gt.799 .and . I NDEX* £0 * 2 9 GO TO 11 
1610 GOTO 12 


1620 11 continue 

16 30 1 1 1 = 8" I 

1640 ii=iii-7 

1650 IFCr .E0.1)PRINT105, I 

16 6 0 I F C I . EQ . 3 9 PR 1 NT 10 6 , 1 
1670 IFCI.E0.49PRINT106,I 
l680 IFCI .En.l09'‘’RTNT110, 
1690 IFCI .e6.139pRI<JT111, 
1700 IFCI .En,299PniNTll2, 
1710 IFCI .eo.179pRI'NT 107, 
1720 IFCI. EO . 1 8 9 p R I NT 10 3 , 
173^ I F c T . EP , 19 9 p R I N T 10 9 , 
1740 IFCI .E'^, 20 9:PPIMT106, 
175'^ IpCI, E'^,309PRInT108, 

1 7 6 p ” p. T r ! T 1 01 , 1 , c L 1 s T c . J 9 , 
17?P IFCCI . EO. 129 .OR. Cl .E 


, CltstCu9,j=ti,h 

,CLISTCJ9,U=ir,II 
,CLISTCU9,U = II,n 
I,CLISTCJ9,U=IT,I 
I,CLISTC.J9,U = II,I 
I, CLISTCj 9,J=II,I 
I ,CLISTCJ9,J=n,I 
I ,CLisTc.J9,<j=n ,i 
I,CLISTC.J9,J = II>I 
I ,CLISTCU9,U = n ,I 
i,Clistcj9,.j=ii,I 
.) = II,IIT9,VALCr9 
O. 249 9'’RINTin3 


l9,VALCn 

I9,VALCI9 

I>,VALCI9 

M9,VALCI 

II9,VALCI 

II9,VALCI 

1I9,VALCT 

in,VALCl 

in,vALCi 

II>,VALCI 

II9,VALCT 


, ApTtG0T012 
,B\IDTH;G0T012 
,POINT;GOTOl2 
9,PPID1 A; GOT 01 
9,PRATE;G0T012 
9, RG/; GO TO 12 
9,oE;GOT012 
9,RPEP;G0T012 
9, OR AD; GOT 01 2 

9 , FO'/; CO TO 12 

9,OPFG;GOTOlP 


1780 12 CONTINUE 

1790 IFCCI. EO . 189 . OR. c I . EQ . 30 9 9 PR I NT 103 
1800 3 CONTI riUE 


3.0 LINK MARGIN SUMMARIES 

The link margin summaries are listed in the following pages. 


1810 I fCIMAR. E0,1)PRINT251 

1320 251 FORMATCTo, UNLINK MARGIN9 

1830 I FCIMAR, EO, 09PRrNT252 

1840 252 FORMATCtG, ’LINK MARGIN-MAX LIKELYilOOD DETECTOR’) 
1850 IFCli‘:AR.EO.2)PRINT253/0Ln0P 



•1360 253 FORMATCTe, ’LINK MARGIN, LOOf’ BANO’/I DTi 1= ' , GB; 3, ' Nf;RT7’0 

1362 IFCIM.AR.EO, 3)RRIijT25^l,RFA,’’ND 

1864 254 FQRMAtCT6, ‘LINK MARGI N, PFA, PND ’ , 2010 . 2) 

1870 1=18 
1830 1 1 1 = 8” I 

1890 ii = rri-7 

1900 ” R I NTlO 8 , I , C L I S T C J ) , J = 1 1 , 1 1 1 ) , V ALC O , R^ P 
1910 PR1NT113,DR0SIG,R0SIG 
1920 '’RINT114,DMARG,AMARG 

1930 113 FORMATCT3, 'REOUIREO n^E/PlJLSE • ,T32, ' DS ’ , T42 /f8 . 2, 

19 405 F9. 2,’ P-E’9 

1950 114 FORMATCT3, *32’, T6, ’ilNK MARGIN’/T32, ’ DLi ' ,T42 , F8 . 2, F9 . 2 ) 

196c ■’RINT104 
1970 21 CONTINUE 
1930 MPUTCIO 
1)00 GOTO 22 

2000 101 FORMATCT3/I 2,T6,6 A4, T32,2A4/T42, F8. 2) 

2010 105 FORHATCT3,I2,T6,6a4,T32,2A4,T42/f 8,2,3”F9.1, * fll LLIV/ ATTS ’ 0 
20 20 106 FORflATCT 3, 1 2, T6 , 6 A4 , T 32 , 2 A4 , T 42 , FB , 2 ,6 p ^9 . 0 , ’ MI CRORADI ANS ’ 0 
2030 107 FORHATCT3,I2,T6,6a4/T32,2A4,t42,F8.2,2PF9,1, ’ PERCENT’) 

2040 108 FORMATCT3,I2,t6,6a4,T32,2A4,t42,f8.2,F9.2, ’ P-E’) 

2050 109 FORMATCT3,I2/T6,6a4,T32,2A4,T42,F8.2,F9. 4, ’ WATTS /M2-A-ST ’ ) 

2060 110 FORMATCT3,I2/T6,6a4,T32,2A4/T42/f8.2,F9,2, ’ METERS DIA. ’) 

2070 111 FORMATCT3,I2,T6,6A4,.T32,2A4,T42,F8.2,F9.n/’ r^ULSES/SEC. ’ ) 

2030 112 FORf1ATCT3,I2,T6,6A4,T32,2A4,T42,F8,2,9'’F9.2,’ NANOSECONDS') 
2090 102 FORMATCIM ) 

2100 103 FORMAT C/) 

2110 104 FORMATC//////) 

2120 END 

2200 SUOROOTINE THRESi IC'^FA, <^ND,B ACK,SI G) 

2210 AL'=FA=ALOGCPFA) 

2220 ALPND=ALOGCPND) 

2230 GAM=BACK+5.6”S0RTC3ACiO 
2240 1 DUM=AL0GCGAM/BACK) 

2250 F=GAM”C1.-DUM)-3ACK-ALPFA 
2260 DG=F/DUf1 
2270 GAM=GAM+DG 

2280 IFCABSCDG/GAM).GT.l. E~4)G0T01 
2290 S I G=GAM-BACK+5 . "SORTCGAM) 

2300 2 DUM=GAf1/CBACK+SIG) 

2310 F=GAM-3ACK-SIG-GAM”AL0GCDUM)-ALPND 

2320 FD=DUM-1, 

2330- DS«»FAi=D ' 

2340 SIG=SIG+DS 

2350 IFCABSCDS/5IG).GT<1.E-4)G0T02 
2360 RETURN 
2370 END 
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lATRS TO LATSCPASSIVE) 

CinCUIT MARGIM SUMMARY 

OPTICAL WAVcLEWGTH ,53 MI CP.OMETERS , RANGE 30000, KM 
PULSB/IDTH= .20 NANOSECONDS 



PARAMETER 

UNITS 

VALUE 



1 

TRANSMITTER POV/ER 

DBV/ 

-20,00 

10,0 

MI LLI'IATTS 

2 

TRANSMITTER L0SS£S 

DB 

-.97 



3 

TRANSMIT ANTENNA GAIN 

DB 

121.07 

5. 

MI CPORADI ANS 

'4 

POINTING LOSS 

03 

-1,39 

1. 

MICROPADT ANS 

5 

FREE SPACE LOSS 

D'i 

-297.04 



5 

AT MO S P H E R I C LOS S C DOl /N ) 

DB 

.00 



7 

TARGET GAIN PRODUCT 

DB 

209.00 



8 

ATMOS P 4 E R I C LOSS CUP ) 

DB 

.00 



9 

FREE SPACE LOSS 

DB 

-297,04 



10 

RECEIVE ANTENNA GAIN 

DB 

131. 00 

,61 

METERS DIA. 

11 

RE CE I VEP> LOSSES 

DJ 

-3.98 



12 

PvECEIVED SIGNAL pQUER 

D3U 

-159.35 



13 

TRANSMIT PULSE RATE 

DB/SEC 

10.00 

10. 

PULSE 3 /SEC, 

1^ 

EflERGY PER RULSE 

DBJ 

-159.35 



IJ 

ENERGY PER p.MOTON 

DB(J 

-134.20 



16 

RE CE I VE D PI 10 rOiMS /P'.JLSE 

DB 

14.91 



17 

OUAfJTUM EFFICIENCY 

' DB 

-6.02 

25.0 

PERCE rs!T 

13 

RECEIVED P-E/ PULSE 

DB 

8.89 

7. 75 

P-E 


19 

B ACKGROUND RADI ANCE 

DB\'/ 

-17.70 

.0170 

WATTS/ M2 -A- ST 

20 

RECEI VER 

FOVCSTERRAD. ) 

DB 

-87.0 7 

50. 

MICOORADI ANS 

21 

OPTICAL 

FILTER BAND’/IDT.i 

DB 

6,99 



22 

RECEIVE 

ANTENNA AREA 

D3M2 

-5.51 



23 

RECEIVER 

LOSSES 

DB 

-3.9 8 



24 

RECEI VED 

BACK, POV/ER 

.DBW 

-107.26 



25 

E 4ERGY PER PMOTON 

DBJ 

-184.26 



26 

RECEIVED 

PMOTONS/SEC 

DB/SEC 

77.00 



27 

OUANTUM 

EFFieiENCY 

DB 

-6.P2 



28 

BACKGROUND P-E /SEC 

DB/SEC 

70.9 8 



29 

RECEIVER 

GATE WIDTH 

DB-SEC 

-86.99 

2,00 

MANOSHCONDS 

30 

BACKGROUND P-E/GATE 

DB 

-16.01 

.03 

P-E 


LI NIC MARGIN 




18 RECEIVED P-E /PULSE 

DB 

8.89 

7.75 '’-E 

31 RED U I RED P-E/ PULSE 

DB 

.00 

1.00 P-E 

32 LINK MARGIN 

DB 

3.89 

7.75 



HATRS TO GROUND CP ASS I VE) 

CIRCUIT MARGIN SUMMARY 

OPTICAL WAVELENGTI-I .5? MI CROMETERS .RANGE 2^1000. KM 


PULSE\/IDTM= ,20 NANOSECONDS 
PARAMETER 

1 TRANSMITTER POV/ER 

2 TRANSMITTER LOSSES 

3 TRANSMIT ANTENNA GAIN 

4 POINTIi'JG LOSS 

5 FREE SPACE LOSS 

6 ATMOS PliERLC LOSS-CDOUN) 

7 TARGET gain PROnUCT 

8 ATMOSPilERIC LOSS CUP) 

9 FREE SPACE LOSS 

in RECEI'/i; ANTENNA GAIN 

11 RECEIVER LOSSES 

12 RECEIVED SIGNAL pOV/ER 

13 TRANSMIT PULSE RATE 

14 E JERGY PER PULSE 
19 ENERGY PER PilOTON 

I-') RfiCFilVED PHOTO.I5/PULSE 

17 OUANTUM EFflCIENCY 

18 RECEI'/ED P-E/PULSE 


19 liACKGROUND RADIANCE 

20 RECEIVER FOVCSTERRAD. ) 

21 OPTICAL FILTER BANnWIDTM 

22 RECEIVE ANTEN.JA AREA ' 

23 receiver losses 

24 RECEIVED BACK. POV/ER 

29 E./ERGY PER pmOTON 

26 RECEIVED PHOTOf^S/SEC 

27 OUANTUM EFFICIENCY 

28 BACKGROUND P-E/SEC 

29 RECEIVER GATE V/IDTN 

30 iiACKGROUilD P-E/GATE 


LINK MARGItl 
18 RECEIVED P-E/PULSE 

31 REOUIRED P-E/PULSe 

32 LI 'IK MARGIN 


UNITS 

VALUE 


DBW 

-20,00 

10,0 

DB 

-.97 


DB 

121.07 

5. 

DB 

-1.39 

1. 

DB 

-295.10 


PB- 

, -U65 


DB 

210.00 


DB 

-4.6G 


DB 

-295.10 


DB 

131.00 

,6l 

DB 

-3.98 


DBW 

-160.78 


DB/SEC 

10,00 

10. 

DBvl 

-170 .78 


OBJ 

-184.26 


DB 

13,4 8 


PB 

-6.0 2 

25,0 

DB 

7. 46 

5.57 


now 

-17.70 

.0170 

DB 

-87.07 

50. 

DB , 

6.09 ' 


DO M2 

-5.51 


DB 

-2.0 8 


DUV/ 

-107.26 


DON 

-184.26 


DB/SEC 

77.00 


DB 

-6.02 


DB/SEC 

70.98 


Db-SEC 

-86.99 

2.00 

DB 

-15.01 

■ .03 


DB 

7.46 

5,57 

DB 

.00 

1.00 

DB 

7.46 

5.57 


D-7 


MILLIWATTS 

MI CRORADT ANS 
MI CRO RADI AMS 


METERS DIA, 


PULSES /SEC. 


PERCENT 

P-E 


WATTS/ M2 -A- ST 
MI CROP AD I ANS 


NANOSECONDS 

P-E 


P-E 



iATRS TO LAT5(PASSIVt;> 

CIRCUIT MARGIN SUMMARY 

OPTI CAL 'UAVtiLtiNGTH .33 HI CROMETERS , RANGE 30000, KH 

pulse'./idtm=4,oo nanoseconds 



parameter 

UN I TS 

VALUE 



1 

transhitter pov/er 

DRW 

•10 ,00 

100. n 

MILLIV/ATT5 

2 

TRANSfilTTER LOSSES 

DB 

-,07 



3 

TRANSMIT ANTENNA GAIN 

DB 

121,07 

5. 

.MI CRORADT ANS 


P0I:JTING LOSS 

DQ 

-1, 39 

1, 

M.I CRORADT ANS 

5 

FREE S'^ACE LOSS 

DB 

-297, on 



6 

A Tf 1 0 S P: I E R I C loss C D O’. / N ) 

DQ 

. ,on 



7 

TARGET GAIN PRODUCT 

DB 

209,00 



P 

■j 

ATriOS ’’MERI C LOSS CUP) 

DB 

,00 



■1 

FREE SPACE LOSS 

DB 

-297, on 



10 

PECEIVd A;^TENNA GAIN 

DB 

135,30 

l.PO 

METERS PI A. 

11 

RECEIVER LOSSES 

DI3 

-3.9B 



12 

RECEIVED SIGNAL PO'/ER 

DBV/ 

-145.05 



1 7 

-V 

TRANSMIT PULSE RATE 

DB/SEC 

10,00 

■ 1-, 

PULSFS/SEC, 

in 

E’lERG'' PER PULSE 

DBvJ 

-155,05 



13 

ENERGY PER PHOTON 

DLi.l 

-ian,2c: 



15 

Rii CE I VE D PM OT 0:iS / UL S E 

• DB 

29,21 



17 

OUA.NTUM cFFICIE.ICY 

DB 

-5,02 

25.0 

PERCENT 

IS 

RECEI VED P-E/”ULSE 

no 

23,19 

203, 51 

P-E 


1 ) 

BACRGROUND RADIANCE 

DBW 

-17,7'- 

,0170 

'.7ATTS/H2-A-S 

2 -'^ 

RECEI'Offt FOVCSTERRAD. ) 

DB 

-87.07 

50. 

HI CRORADT ANS 

21 

OPTICAL FILTER BAND'./IDTM 

DB 

6,99 ' 



22 

RECEIVE ANTENNA AREA 

DBM2 

-1,21 



00 

j 

RECEIVER LOSSES 

DB 

-3.93 



2n 

RECEI VED BACi;, PO’VER 

■DBW 

-102,95 



25 

ENERGY PER PHOTON 

D B J 

-I8n,26 



25 

R E CE I VE ;D P ! \ OT ON S /SEC 

Du /SEC 

31.30 



27 

OUANTUM EFFICIENCY 

DB 

-5, 02 



28 

BACKGROUND p- E/SEC 

DB /SEC 

75, 28 



20 

RECEIVER GATE WIDTH 

D3-SEC 

-83.01 

5.00 

MANOSECOrJDS 

30 

B ACKGRO UN D P- E /GAT E 

DB 

-7.73 

■ ,17 

P-E 



LINK MARGIN 




18 

RECEIVED P-E/PULSE 

DB 

23.19 

20R.51 o-E 

31 

RED U I RE D P-E / PULS E 

DB 

in.77 

30,00 p-e 

32 

LINK MARGIN 

D3 

8.n2 

6,95 


D-8 



lATRS TO GROUNDCPASSI VE) 

CIRCUIT , MARGIN SUMMARY 

OPTICAL WAVELENGTH .53 MI CROMETERS , RANGE 2^000. KM 
'’ULSEV/IDTHsil.OO NANOSECONDS 


PARAMETER 

1 TRANSMITTER POWER 

2 TR ANS M I T TE R LOS S E S 

3 TRANSMIT ANTENNA GAIN 

4 POINTING LOSS 

5 FREE SPACE LOSS 

6 ATMOSPHERIC LOSSCDn'/N) 

7 TARGET GAIN PRODUCT 
3 ATMOSPHERIC LOSSCUP) 

9 FREE SPACE LOSS 

10 RE CE I VE ANT ENN A G A I 

11 RECEIVER LOSSES 

12 RECEIVED SIGNAL POWER 

13 TRANSMIT pulse RATE 

14 ENERGY HER ^ULSE 
Ij ENERGY PER PHOTON 

16 RE GE I VE D PHOTOiNS /• P UL S £ 

17 OUANTUM efficiency 
13 RECEIVED P-E/'’ULSE 


19 I3ACKGROUND RADIANCE 

■ 20 RECEIVEP. FO VC ST ERR AD. ) 

21 OPTICAL' FILTER BAtID’IIDTH 
3? r>EC£IVE antenna AREA 

23 RECE I VER ' losses 

24 RECEIVED BACK. ROWEP. 

25 ENERGY PER PHOTON 

26 RECEIVED PHOTONS /SEC 

27 OUANTUM EFFICIENCY 

20 BACKGROUND P-E/SEC 
29 RECEIVER GATE WIDTH 
3'J' BACKGROUND P-E/ GATE 


LINK MARGIN 
IS RECEIVED P-E/PULSE 

31 REOUIRED P-E/PULSE 

32 Lirk MARGIiJ 


UNI TS 

VALUE 



DBW 

-10.00 

100 . 0 

MILLIWATTS 

DB 

-,97- 



DB 

121.07 

5. 

Ml CRORADI ANS 

DB 

-1. 39 

1. 

MI CRORADI ANS 

DB 

-295.10 



03 

-1.65 

• 


DB 

210 . 00 



DB 

-4.66 



DB 

-295.10 



DB 

135.30 

1.00 

meters dia. 

DB 

-3.98 



DBW . ■ 

-146.43 

/ 


DB/SEC 

10.00 

10. 

"ULSES/SEC. 

DBU 

-156.43 



OBJ 

-134.26 



DB 

2 7,. 7 3 



DB 

-6.0 2 

25.0 

"ERCENT 

DB 

21.76 

149.92 

P-E 


DBW 

-17.70 

.0170 WATTS /M2-A-ST 

DB 

-37,07 

50, MI CRORADI AMS 

DB 

6,0 0 


DBM2 

-1.21 


DB 

-3.93 


DBW 

-102.96 


D3U 

-184.26 


DB/SEC 

31,30 


DB 

-6.0? 


DB/SEC 

75.23 


DB-SEC 

-83,01 

5,00 NANOSECONDS 

DB 

-7,73 

.17 p-E 


03 

21.76 

149.92 P-E 

DB 

1^.77 

30 ,00 r>_5 

DB 

6,99 

'5.00 


D-9 


tr 



SWUTTLi-; TO 0R01J(JD 
CIRCUIT MARGIN SUMMARY 


0^ TI CAL V/ AVE LENGTH . 5 3 M I C ROME TERS . RANGE 

350. 

KM 

P'JLSEV/IDTM= .20 NANOSECONDS 




PARAMETER 

U.NITS 

VALUE 


1 TRANSMITTER POV/ER 

D3V/ 

- 30.00 

1.0 MILLIWATTS 

2 TRANSMITTER LOSSES 

Dii 

-.97 


3 TRANSfMT ANTENNA GAIN 

DB 

83,57 

375, MICRO RAD IANS 

4 POINTING LOSS 

DU 

-.62 

50 , MICRORADIANS 

5 FREE SPACE LOSS 

03 

-258.38 


6 ATMOSPHERIC LOSS CDO‘,/^J ) 

03 

- 1.65 


7 TARGET GAIfJ PRODUCT 

DB 

196.00 


3 ATMOSPHERIC LOSSCUP) 

DB 

-4.66 


9 free SPACE LOSS 

DB 

- 258 . p 8 


10 RECEIVE ANTENNA GAIN 

DB 

125.01 

.30 METERS DI A. 

11 RECEIVER LOSSES 

DB • 

- 3.98 


12 RECEIVED SIGNAL POWER 

DBW 

-154.05 


13 TRANSMIT PULSE RATE 

DB/SEC 

.00 

1. PULS ES /SEC. 

14 EiJERGY PER PULSE 

DUJ 

- 154,05 


15 ENERGY PER PHOTON 

DBU 

- 184.26 


16 RECEIVED PHOTONS /PULSE 

. 03 

30 , 21 


17 QUANTUM efficiency 

DB 

-6.02 

25.0 PERCENT 

18 RECEIVED P-E/PULSE 

DB 

24.19 

262.24 P-E 


19 

B A CKG RO UN D RAD I AN C E 

DB>/ 

-17,70 

.0170 

./ATTS/M2-A-ST 

20 

R E C E I V EP FO V C S T E R R A D . ) 

DB 

- 69.50 

375 . 

MI CRORADT AMS 

21 

OPTICAL FILTER BANDWIDT'1 

DB 

6.99 ' 



22 

RECEIVE ANTE.'INA AREA 

D3M2 

-11.50 



23 

RECEIVER LOSSES 

DB 

- 3.98 



24 

RECEIVED BACK. POWER 

D3W 

-95.75 



25 

E:1ERGY PER PHOTON 

OBU 

-184.26 



26 

RECEIVED PHOTGN5/SEC 

DU /SEC 

88.51 



-Z 

OUANTUM EFFICIENCY 

DB 




23 

3 ACKG R 0 UNO P- E / S E C 

DB/SEC 

82,49 



29 

RECEIVER GATE WIDTH 

DB-SEC 

- 96.02 

.25 

^•!A^^OSECONDS 

30 

BACKGROUND P-E /GATE 

DB 

-13.53 

,04 

P-E 


LINK MARGIN-MAX LIKELYHOOD DETECTOR 


- 


18 

RECEIVED P-E/PULSE 

03 

24.19 

062.24 

P-E 

31 

R E 0 U I R CD P - £ / PI ) L S £ 

DB 

7.08 

5.11 

P-E 

32 

LINK MARGII4 

DB 

17.10 

51.-32 



D-10 



^.lUTTLE TO GROUND 
CIRCUIT MARGIN SUMMARY 

ORTICAL V/AVe LENGTH 1,06 Ml CROMETERS , RANGE ^5^, !>1 

PULScWIOTHs ,20 NANOSECONDS 


PARAMETER 

UNITS 

VALUE 



1 TRANSMITTEr; PO'./ER 

DEW 

-23,08 

^,0 

■u LLrnTTS 

2 TRANSMITTER LOSSES 

D8 




3 TRANSMIT ANTENNA GAIN 

D3 

77,55 

750 . 

MICRO RAD IANS 

n POINTING LOSS 

DE 

-,15 

50 , 

MI CRORADI AMS 

5 FREE SPACE LOSS 

Dd 

-252,36 



6 ATMOSPHERIC LOSS C DOWN) 

DB 

-,56 



7 TARGET GAI(>1 PRODUCT 

DB 

190.00 



8 ATMOSPHERIC LOSS C UP) 

OB 

-3.57 



') FREE SPACE LOSS 

DB 

- 252,36 



10 RflCEIVE ANTENNA GAIN 

DB 

118.99 

.30 

METERS DIA, 

11 PECEIVER LOSSES 

DB 

- 3,98 



12 RECEIVED SIGNAL POWER 

DOW 

-151,38 



13 TRANSMIT PULSE RATE 

D3/SEC 

.00 

1, 

"ULSES/SEC, 

14 ENERGY PER PULSE 

DBJ 

-151.3 8 



15 energy per PHOTON 

DBv) 

-187.27 



16 RECEIVED PHOTONS/PULSE 

DB 

35.89 



17 OUANTUM efficiency 

' DB 

-20,00 

I.O 

PERCENT 

18 RECEIVED P-E/ PULSE 

DB 

15.89 

33,82 

P-E 


19 

BACKGROUND RADIANCE 

DBW 

- 20.58 

.0087 

’,/ATTS/M2-A-ST 

20 

RECEIVEB FOVCSTERRAD, ) 

DB 

- 63.55 

750 . 

MI CRORADI AMS 

21 

OPTICAL FILTER BANDWIDT'1 

DB 

6,99 ■ 



22 

RECEIVE ANTENNA AREA 

DBM2 

-11,50 



23 

RECEIVER LOSSES 

D3 

- 3,98 



24 

RECEIVED BACK, POWER 

■DBW 

-92.61 



25 

ENERGY PER PHOTON 

DBU 

-187.27 



26 

RECEIVED PilOTONS/SEC 

DB/SEC 

94.66 



27 

OUANTUM EFFICIENCY 

DB 

-20,00 



28 

3 ACKG ROUrj D P- E / S EC 

DB/SEC 

74.66 



29 

RECEIVER GATE WIDTH 

DB-SEC 

- 96.02 

.25 

f'MNOSECONDS 

30 

BACKGROUND P-E /GATE 

DB 

- 21,36 

,01 

— F 


18 

LINK MARGIN-MAX LIKELYHOOD DETECTOR 
RECEIVED P-E/PULSE DB 

15 . 89 

38. 82 P-E 

31 

REOUIRED P-E/PULSE DB 

7.05 

5,07 P-E 

32 

LINK MARGIN ’ DB 

8. 84 

7.65 


D-n: 


L<) 



S:iUTTLE TO GROUND 
CIRCUIT MARGIN SUMMARY 

OPTICAL WAVELENGTH ,53 MI CROMETERS/RANGE 350, KM 
P'JLSE'/IDTH=i|,00 NANOSECONDS 


PARAMETER 

UNI TS 

VALUE 



1 TRANSMITTER "OWER 

DBW 

-20,00 

10.0 

fllLLIW.ATTS 

2 TRANSMITTER LOSSES 

DLi 

-.97 



3 TRANSMIT ANTEl'INA GAIN 

DR 

83.57 

375. 

MICRORADI ANS 

4 POniTING LOSS 

DB 

-.62 

50. 

fll CRORADT ANS 

5 FREE S'^ACE LOSS 

OB 

-258.38 



6 ATMOSPHERIC LOSSCDOI/N) 

DB 

-1.65 



7 TARGET GAIN PRODUCT 

DB 

196,00 



8 ATMOSPHERIC LOSSCUP) 

03 

-4,66 



9 FREE SPACE LOSS 

DB 

-258,38 



10 RECEIVE ANTENNA GAIN 

DB 

128,53 

. 46 

f-IETEPS DIA. 

11 RECEIVER LOSSES 

DB 

-3.98 



12 RECEIVED SIGNAL POWER 

DBW ■ 

-140.53 



13 TRANSMIT PULSE RATE 

DB/SEC 

.00 

1. 

PULSES /SEC. 

14 ENERGY PER PULSE 

D3J 

-140,53 



15 ENERGY PEP. PHOTON 

DB.J 

-184,26 



16 RECEIVED PHOTONS/^ULSE 

DB 

43,73 



17 OUANTUM EFFICIENCY 

■DB 

-6,02 

25.0 

PERCENT 

18 RECEIVED p-i-E/pULSE 

DB 

37.71 

5900,40 

p-E 


19 

B ACKGROUND RADI ANCE 

DBW 

-17.70 

.017^ 

WATTS /f12-A-ST 

20 

RECEI VER FOVCSTERRAD. ) 

DB 

-69.57 

375, 

MICRORADI AMS 

21 

OPT! CAL ' FILTER BANDWIDTH 

DB 

6.90 



22 

RECEIVE ANTENNA AREA 

DBM2 

-7.98 



23 

RECEIVER LOSSES 

DB 

-3.98 



24 

DECEIVED BACK. POWER 

DBW 

-92.23 



25 

EHEPGY P£R PHOTON 

OBJ 

-184.26 



26 

RECEI '.'ED PHOTONS /SEC 

DB/SEC 

92.03 



27 

OUANTUM EFPICIENCY 

DB 

-6.02 



28 

aACKGR{)UND P-E/SEC 

DB/SEC 

86.01 



29 

RECEI '/ER GATE WI DTH 

DB-SEC 

-83.01 


NANOSECONDS 

30 

B ACI:GR0U.ID p-E/GAT e 

DB 

3.00 

2 ,PP 



LT.NK MARGIN-MAX LIKELYHOOD DETECTOR 


- 


18 

RECEIVED P-E /PULSE 

DB 

37,71 

c^O P O . 4 p 

n^rr 

31 

PEOUI-RED P-E/PULSE 

DB . 

22i,66 

'i84.3^ 


32 

LINK MARGIN 

DB 

15.05 

32.00 



v/) ij'i 



SHUTTLu TQ (GROUND 
CIRC'JTT margin 5'JflMARY 


PTICAL WAVELENGTH 1.0 6 MICR 

OMETERS, RANGE 

350. 

KM 


ulse'Mdt'I-=4,oo nanoseconds 






parameter 

UNITS 

VALUE 



1 

transmitter power 

DBW 

“13,98 

.tr 

0 

0 

MILLIWATTS 

2 

TRANSMITTER LOSSES 

DB 

-,97 



3 

TRANSMIT ANTENNA GAIN 

DB 

77.55 

750 , 

MICRORADIAMS 

i| 

POINTING LOSS 

DB 

-.15 

50 , 

MICRORADIANS 

5 

FREE S'^ACE LOSS 

DB 

- 252,36 



6 

ATMOS PHERI C LOSS COOWN ) 

D3 

-,56 



7 

TARGET GAIN PRODUCT 

DB 

190,00 



o 

' i 

AT lOSPHERI C LOSSCU'^O 

DB 

-3.57 



9 

FPEE SPACE LOSS 

DB 

- 252.36 



10 

PJICEIVE ANTENNA GAIN 

DB 

122,51 

.46 

METERS DTA. 

11 

RECEIVER LOSSES 

DB 

- 3,93 



12 

RECEIVED SIGNAL POWER 

DBW 

- 137.86 



13 

TRANSMIT PULSE RATE 

DB/SEC 

,00 

1. 

PULSES /SEC. 

14 

ENERGY PER PULSE 

DBaI 

-137,86 



15 

ENERGY PER PHOTON 

DBJ 

-187.27 



16 

RECEIVED PHOTONS/PULSE 

•DS 

40.41 



17 

O'JANTUM EFFICIENCY 

DB 

-20,00 

1,0 

PERCENT. 

18 

RECEIVED P-E/PULSE 

DB 

29.41 

873,46 

P-E 

19 

BACKGROUND RADIANCE 

DBW 

- 20.58 

.0087 

WATTS /M2-A-ST 

20 

RECEIVER FOVCSTERRAD. ) 

DB 

- 63.55 

750 . 

MICRORADIANS 

21 

OPTICAL FILTER UANDWIDTH 

D3 

6.99 



22 

RECEIVE ANTENNA AREA 

DBM2 

-7,98 



23 

RECEIVER LOSSES 

DB 

- 3,93 



24 

RECEIVED BACK. POWER 

DBW 

- 89 , 09 



25 

ErJERGY PER PHOTON 

DBJ 

-187,27 



26 

RECEIVED PHOTONS /SEC 

DB/SEC 

98.18 



27 

QUANTUM EFFICIENCY 

DB 

-20,00 



28 

BACKGROUND P-E/SEC 

DB/SEC 

78.18 



29 

RECEIVER GATE WIDTH 

DB-SEC 

- 83,01 

5.00 

r,'ANOSECONDS 

30 

BACKGROUND P-E/GATE 

DB 

- 4,80 

.33 

p— E 


LINK MARGIN-MAX LIKELYHOOD DETECTOR 




18 

RECEIVED P-E/PULSE 

DB 

29.41 

873.46 

P—5 

31 

REOIJIRED P-E/PULSE 

DB 

22.62 

182.72 

p-E 

32 

LIfjK MARGirj 

DB 

6.79 

4.78 



D-13 



GROUrJD BEACON TO SHUTTLE 
CinCUIT MARGIN SUMMARY 


OPTICAL V/AVELENGTH .88 MICRO'^ 

lETEPS^ R ANGE 

^6 It, 

UM 


’’ULGE'/I DTHS5, on NANOSECOMOS 






PARAMETER 

UNITS 

VALUE 



1 

TR ANS M I T TER P OV/E R 

PSV/ 

-3’.ni 

.5 

MILLI/ATTS 

p 

TPANS.MITTEr LOSSES 
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DB 
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7.85 
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APPENDIX E 


PRECISE OPTICAL PULSE TIMING 


1.0 INTRODUCTION 

The data in this Appendix presents fundamental results of investigations 
into optimal range estimation techniques, and was partially funded under this 
Study. This material was presented in IEEE Southeastcon '75, in April, 1975, 
and published in the proceedings. The reprint is included in this Appendix 
for completeness. 
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PRECISE OPTICAL PULSE TIMING* 
by Gary Lee and George Schroeder ^ 
McDonnelLDouglas Astronautics Company 
St. Louis, Missouri 


ABSTRACT 


A(t, cx) - S(t, a) + 


ai 


This paper evaluates the performance of the ML and several 
sub-optimum arrival time estimators for optical pulses. The per- 
formance analysis of a “sliding window'* counter followed by a 
threshold detector is achieved by deriving an expression for the 
first crossing density of the output of the counter. A sub-optimum 
arrival time estimator using a threshold rather than a maximum 
detector is suggested and evaluated by simulation. It appears to 
prtn ide near optimum performance. 


INTRODUCTION 


Lasers permit the use of very narrow pulse widths for both 
communication and ranging. The combination of very narrow 
pulses and narrow beam widths offers the promise of heretofore 
unachievable efficiencies in terms of the power required to achieve 
a given data rate or ranging accuracy. 

Pulse position modulation allows the conversion of precise pulse 
arrival time resolution into high data rates with minimal 
transmitted power per pulse. If reasonable data rates are to be 
achieved, it is necessary that the pulse detection technique not be 
so complex that it limits the rate at which pulses can be detected. 
Thus it is desired to achieve the best possible pulse resolution 
without limiting the usable piulse rate. The number of received 
signal photoelectrons per pulse will always be small (i.e. on the 
order of 50 or less) so that analyses of achievable pulse resolution 
with non-coherent detection depend on the properties of non-sta- 
tionary Poisson processes. 

Ranging requires lower pulse rates than communication and 
thus may, depending on the application, use more complex pulse 
detection techniques. Some programs (such as the NASA EOPAP 
program] require satellite borne ranging systems with very precise 
(centimeter] range resolution. These systems must be small, light 
and reliable but yet capable of less than a hundred picosecond 
pulse arrival time resolution with only a few received signal photo- 
electrons per pulse. 

This paper defines achievable optical pulse resolution as a 
function of the pulse shape, received signal and background power, 
and the pulse detection technique. Optimum (maximum 
likelihood) detectors, simple counters and weighted counters are 
evaluated both analytically and by simulation. A novel 
approximate solution for the first crossing density of the output of 
a “finite window” counter or integrator is used to define both false 
alarm rates and pulse resolution properties for the “finite window” 
counter case. 


ANALYSES 

It is desired to determine the achievable pulse arrival time 
resolution for a single 200 picosecond optical pulse. Two hundred 
picoseconds is about the minimum presently achievable pulse 
width from a Nd: YAG laser. Since non-coherent detection must be 
used, the output of the optical detector will be assumed to be a 
non-stationary Poisson process. Thus the output of the optical 
detector consists of a Poisson process whose mean rate A(t, a ) is 
given by 


where S(t, a ) represents the contribution of the signal pulse and 
u\y represents the uniform background level, a is the unknown 
pulse arrival time and (0, Tjnaxl is the total observation interval. 

Hoversten ct al 1 , 2 obtained a partial differential equation for 
the conditional density of a using non-linear filter theory. These 
results, although useful for developing insight into optimal 
estimator structures, are not directly implemen table into hardware 
without a considerable number of approximations. 

Bar-David 3 derived an expression for the maximum likelihood 
estimate of a conditioned upon the observations obtained during a 
fixed observation interval. Let jtj, t 2 , . . . t^} = {Tj^i f be the set of 
observed photoelectron arrival times. Bar-David showed that the 
probability of jT]^}, given a , is 

M 

( 2 ) 

P {T^} - exp[-Q(a)] n x(t , a) 
where 

T 

/ max 

A(t, a) dt 

0 


The maximum likelihood estimate of a is the value which 
maximizes the value of P|Tm} or equivalently log PjT|^|. 
Bar-David also d eterm ined the following approximate expression 
for the variance, £a^;Of the estimate of a for the special case of 
S(t,a I everywhere differentiable and S(t, -T/2] = S(t) = 0 for \t\ 
>T/2. 




S(t) + Oj, 


1-1 


dt 


(3] 


Figure 1 shows a block diagram for generating log P jTjvi} and 
the ML (maximum likelihood) estimates of a . It is probably 
possible to implement the ML detector digitally for the case where 
the average number of received photoelectrons per pulse is small, 
say less than 20. This could be accomplished by detecting the 
arrival time of individual photoelectrons (i.e. {ti. t 2 » . . t^jj), 

digitizing the arrival times, and using a computer to determine the 
value of a which maximizes P{Tm}* Individual arrival times can 
be resolved as close as 5 picoseconds using an optical detector with 
an image converter. An image converter is a rather complex device 
which first uses a cathode ray tube face to map and store an image 
of a received pulse. The tube face is then scanned and the output is 
sampled and A/D converted with a spatial sampling rate which 
corresponds to a 5 picosecond sampling rate across the original 
pulse. The tube has sufficient gain to detect the arrival of an 
individual photon within a resolution element. Note that this 
approach for implementing the ML detector is neither simple nor 

tGeorge Schroeder is also with the Department of Engineering, 
Southern Illinois University at Edwardsville. 

♦This work was partially supported by NAS5-20646. 
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capable of being performed rapidly enough not to be the limiting 
factor in determining data rate for an optical PPM communication 
system. 



FIGURE 1 MAXIMUM LIKELIHOOD ESTIMATOR 


where E ( • ) represents the expectation with respect to the time and 
number of occurrences to the input process. The characteristic 
functions of I''('Qq) and Qq) ate found using the expectation of 
product functions with X(t) » S(t — a ) + njj 

tfexp [ix F(a)]J - Efexp [-ixJ^rCt. -a)]} 

j ^ 19) 

- EflT exp [-ix r(t - a)]} 

J ^ 


Analog maximum detectors, though conceptually simple, are 
seldom practical. The basic approach is to first differentiate the 
signal representing log P|T]vi|and then use a threshold detector to 
determine when the derivative passes through zero. The value of 
the signal and the time of occurrence of these relative extrema are 
sampled and stored so that the time of the absolute maximum can 
be determined at the end of the measurement period. Note that for 
a small number of photoelectrons per pulse there will be several 
relative maxima per pulse. Thus the sampler must be capable of 
taking samples on the order of a few picoseconds apart in order 
that the detector look like an “ideal” maximum detector. This, if 
possible at all. would be very difficult. 

'Phe preceding discussion has been mainly concerned with the 
practical difficulties involved in building maximum detectors. 
Another area of difficulty with the ML detector is the uncertainty 
in knowledge of S{t). Laser amplitude instability and multi-path 
(for the ranging case) limit the ability of the receiver designer to 
specify S(t). Note also that even if S(t) is known exactly it may be 
\ ery difficult to build a filter with several gigahertz bandwidth and 
an impulse response matched to log [S|t| + niJ. 

Thus, in summary, although it is probably possible to 
implement the ML detector at some cost in size, weight, and pulse 
ra te, it is clear that this should only be attempted if the achievable 
gain in resolution of the ML detector over other simpler 
approaches is substantial. The goal of this paper is to define this 
gain so that the tradeoff can be made by each designer for his own 
application. 


This expectation can be evaluated in terms of 


and 


/max . MO) 

G “ / fS(t - Oq) + n^] exp {-ixr(t - a^) ] dt 

0 X 

^raax 

Q - / x(t)dt Ml) 

^ 0 


Neglecting end effects, let, Qq = Tmax/2. Then, with 


r — t - Tjnax/2 


and 


G = 


max 

2 


i 

/’ 

-T 


[S(t) + exp [-ixr(x)]dT 


( 12 ) 


max 

2 


X 

^ max 

J ^ [S(t) + 


n^]dT 


(13) 


max 

2 


EFFECTS OF FILTER MISMATCH 

First, the effect on resolution of mismatch between the signal 
and the filter response in the ML estimator will be determined. 
This is accomplished by deriving an approximate expression for 
the variance in the estimate of a assuming the filter response is 
designed for a signal S’ ( t) when, in fact, the signal is S(t). Since a 
signal pulse shape of S’ (t) is used in the filter design 


and 


r<t) ^ log [S»(t) -I- ri^l 
F(a) » r(t. - a) 


<41 

(5) 


Tj r(t 
j ^ 

F( a ) is the input to the maximum detector and the tj*s are the 
times of occurrence of the received photoelectrons. Freezing in a 
manner similar to Bar-David, an approximation to the variance of 
the estimate^is obtained. 

Expanding F( a) about the actual delay, a q, and using the 
definition of 'ct 


It can be shown that 

E{exp [txF(a)]} * exp (-Q + G) 


exp f“Q + 


T 

I 


(14) 


max 

2 


[S(t) + X ] exp [-ixr(T)]dTl 
n ‘ 


max 

2 


In the same manner, letting 
G 


/ 


max 

2 


[S(t) + n^l exp [ix r(T)]di 


max 

2 


F(a) = 0 =* Ffa^) (a - a^) FCa^) + 


This yields 


Thus 


“O ^ ® , 

FCa^) 


e2^E{[F 


[E{F (a^)}]^ 


( 6 ) 


(7) 


(8) 


E{exp[ixF(a) ] } 


exp 


-Q + 


exp(-Q + G) 
T 

« max 
2 


I ^ [S(t) + n^] exp [ixr(T)]dT 


max 

2 


( 15 ) 


Using these results 
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2 T ^ 

^ ^max 


t 

r 2 

1 [S(T)+n^]r(T)dT 

I {S(x)+t^][r(T)]^dT 

,p(t) -.a(t) exp 

1 a(x)dx 
^ o 

-T 

max 

1 ’’^max 

where 


2 i 

I 2 ^ 




/ T_ 


max 

2 


^2 


/ iUO 

-T 


[S(t) + n^] r(T)dT 


max 

2 


(16) 


U8) 


a(t) = limi Prob. {crossing in (t, t + A ) 1 no crossing before t} 
A-OA 


a(t) can be approxunated by 
a(t) - 


L-1 


k-0 


■“ vk 


kl 


(19) 


Assuming that 


where L is the threshold and 


T \ . / T 

I , max I I , max ■ « 

r ±— 1= r l± — 1= 0 


r T_. 


/ 

“T 


max 

2 


r(T) S(x)dT 


max 

2 


i 

y. max 

-T 


[S(t) 4- n^] [r(T)]^dT 


max 

2 


T 

r 




rCt) S(x)dT 


max 
2 

and if r ( t ) = 0 for t t J > t/ 2 the limits in the integrals of (17) can 
be replaced with ± T/2. 

SLIDING WINDOW COUNTER 

By replacing the maximum detector with a threshold and the 
matched filter with a “sliding window” integrator or counter, a 
simple analog detection system is obtained (see Figure 2). This 
system estimates arrival time by observing the number of photo- 
electrons received during the previous T seconds and thresholding 
this value. Determination of the timing resolution of this system is 
not simple since it requires the calculation of the “first crossing” 
density of the output of a “finite window” counter. To the best of 
the authors knowledge, there are no other results for this problem 
(either exact or approximate). 



FIGURE 2 "SLIDING WINDOW" COUNTER ESTIMATOR 

Lee and Fogle ^ obtained an expression for the “first crossing” 
density of a random process in terms of a function a(t) somewhat 
analogous to the infinitesimal transition probability for Markov 
processes. Appendix A summarizes this development and shows 
how a(t) can be approximated for a Poisson process and a “finite 
window” counter. Note that the optimum threshold value is 
determined by numerical optimization. 

To calculate the variance in threshold crossing time, it is 
necessary to determine the probability density hmction, p(t), of 
the time of the first threshold crossing t . From Appendix A. 


nT 


t 

L 


X(x, a)dx 


( 20 ) 


For the problem under consideration 

0 < t < a - T/2 

+ S(t-a) a - T/2 < t < a + T/2 

% 


( 21 ) 


a + T/2 < t 


117) 


And, ignoring end effects for 0 < t < T 

^ n^T 0 < t < a “ T/2 


uX + 


1. 

/ 


t - a 

S(x)dx a - T/2 < t < a + T/2 
-t/2 


( 22 ) 


nT« 


T/2 


, T + 


/ 




S(x)dx a + T/2 < t < a + 3T/2 
t-a-T 

a + 3T/2 < t 


With a finite observation interval there is a non-zero probability 
of no threshold crossing. Therefore, for comparison purpose we 
determine the mean square error given that there has been a 
threshold crossing in (0, Tinax)* 


Prob {t < t I T < 


Prob . {x < t} 


Prob. {(x < t)n (i < X^^^} 

____ . Tx "<~T T" 

max 


Prob.Cx < T ) 
max 


t < T 


t - T 


(23) 


Thus the conditional density used in finding the expected values 
given a threshold crossing before Tmax is 


Pj Ct)- 
max 




Prob.fx < T } t < T 

max max 


t > T 


(24) 
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Thus the mean threshold crossing time, given there has been a 
crossing before Tmax» is 

T 


® * Prob { 


j f 

t < T } J 

max •'o 


max 

ta(t) exp 


i 


t 

a(x)dx 


(25) 


dt 


and the mean square value of crossing time is 




Prob ii 


Integrating (25) and (26) by parts 
E(t) - 1 _ 

^ W 


T 

ftaax - 

t a(t)exp 

’ t 

- j a(x)dx 

f 

o 

•'o 


dt 


(26) 


(27) 


"T exp 
max 


T 1 r ^ 

- [ ”:<x)dx + j - p<x)dx 

Jo o L 


dt) 


and 


e(t^) « 


J 

O 1 

J 

1 

2 

exp 

T 

rmax 

T } ' 

max 1 

- 1 a(x)dx 

1 max 

Jo 



m 




•max 
2 I t exp 


i: 


a(x)dx 


dt 


where 


Prob {t < T } ^ 
max 



dt 


max 

1 - exp |- / a(x)dx| 


Note that by replacing Tmax signal pulse be^s 

one obtains the probability of an early detection, can easily 
be obtained from equation (27) ~ (29). 

COMPARISON OF ESTIMATION TECHNIQUES 

Comparing the optimum (ML) arrival time estimator shown in 
Figure 1 with the more practical estimator configuration shown in 
Figure 2, two differences become obvious. The first difference is 
that the ML estimator uses a filter that is in some sense matched 
to the signal and background conditions whereas the “sliding 
window” counter configuration does not. The other difference is 
that the ML estimator requires a maximum detector where as the 
“sliding window” counter uses only a simple threshold detector. 

The variance in arrival time estimation will now be computed 
assuming a raised cosine optical pulse shape. Note that this pulse 
shape was chosen in order to meet the conditions on S(t) required 
to derive equation ( 3 ). It is not required for evaluating variance for 
the “sliding window” counter. 


S(t) 


s ^ 2 tt 
— (1 -I- cos — t) 


|t| < T/2 
|t| > T/2 


(30) 


ns is the mean number of received signal photoelectrons per p 
From (3) the mean square timing error using the ML estimaV' 


a * 4 tt 2 [iig + -Vn T (2 + n^^T) ] 


( 31 ) 


To develop an intuitive feeling for how much of the performance 
difference between the “sliding window” counter and the ML 
estimator is due to the difference between a threshold and 
maximum detector and how much is due to mismatch between the 
integrator impulse response and the pulse shape, the performance 
of a mismatched ML estimator was also evaluated. Figure 3 
compares the rms error for the ML estimator configuration just 
evaluated 




r(t)-l«g [(X +COS + ii|^] and S(t) =^(l + cos 

with the rms error for a mismatched ML estimator configuration 
with 


2 ti. 


r(t) 


[1 + cos t] and S(t) = ^(1 + cos |^t) 


(28) 


(29) 


The rms error for the latter case using equation (17) is given by 
n_ + 

(32) 


2 

e = 
a 


o 2 2 

27t n 

s 


"b%2 


Figure (3) also shows some computer simulation points for the 
mismatched case. The simulation is discussed later. 



Note that for this particular case mismatch does not cause much 
degradation from the optimum. 

The “sliding window” counter and ML estimator will now be 
compared. 

For the sliding window counter and a raised cosine pulse shape, 
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X(t, a) x.-/ 


0 < t < a - T/2 

^[i+cos |S- (t-a)]a-T/2<t< a+ T/2 

a + T/2 < t (33) 


% 

V_ 


and 




0<t<«-T/2 


nT=J 


h <,-T/2<t<«+T/2 

(34 

%T+a^[|-t^-|-sin(2*£fa)j «+T/2<t«.+3T/2 




a + 3T/2 < t 


Thus, equations (27) and (28) can now be evaluated numerically 
using equations (19), (29), (33), and (34). 

Figure (4) compares the rms error in estimating arrival time 
(normalized by the pulse width) for the ML estimator and the 
“sliding window” counter for a background level of 0.1 and 1.0 
photoelectrons per pulse width. Note that the ML estimator is 
roughly a factor of 2 superior to the counter in resolution. This 
corresponds to an extra 6 dB of required laser power which is 
definitely significant. 



RECEIVED SIGNAL PHOTOELECTRONS/PULSE (Hj) 

FIGURE 4 COMPARISON OF M.L. AND “SLIDING WINDOW” 
COUNTER ESTIMATORS 


Based on a previous result for threshold detection ® , the 


optimum ranging threshold level for the counter is expected to be 
approximately given by 


L» 


InT 


n. 

s 


-. n T -I- log (T /T) 
D tnax 

Iog?nJ7n^T) 


(35) 


where InT(x) is the smallest integer greater than or equal to x. 

Figure (6) compares L' with the optimum threshold determined 
numerically using equations (27) and (28). Note that L’ appears to 
be a good approximation for the optimum threshold. 


VALIDATION BY SIMULATION 

Since the derivation of both the expressions for resolution of the 
ML estimator and the ‘‘sliding window” counter contain approxi- 
mations, a computer simulation was performed to verify the 
analytical results. 

The first step in the simulation is to generate the sequence of 
photoelectron arrival times. This is done by first dividing t he pulse 
width, T, into 500 intervals each of width A — 3V500. The 
probability of receiving n photoelectrons in an interval centered 
about time t is approximately 


n 


= 0.1 


= 1.0 
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FIGURE 5 COMPARISON OF EXACT AND APPROXIMATE THRESHOLDS 


P(n) ~ exp [-AA(t)] • [AA(t)]^/n! (36) 


(The approximation is due to the fact that A(t) is not constant 
over A ,) For each interval a uniform random number, X, is 

n~l n 


A value n is then determined such that 


<Yj 


,p(j) 


j=0 j=0 


n is then the number of photoelectrons received in the interval. For 
our case with A = T/500, n is never greater than 2 and is usually 
1 or 0. When n = 2 occurs it is treated as a single occurrence of 
double strength. 

The signal pulse is centered in an interval of lengt h 4T. Since the 
effect of false alarms or anomalous errors on rms timing error is 
not considered in the analysis of the ML estimator performance, it 
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was desired to eliminate the effect of these errors on rms timing 
error of the ‘‘sliding window” counter so that a direct comparison 
could be made. The choice of Tmax ^ ^ compromise 

between minimizing these errors and eliminating end effects. 

In the simulation of the “sliding window” counter, the number 
of photoelectrons received over the last T seconds are added and 
when the sum exceeds a set threshold, the threshold crossing time 
is used as an estimate of the pulse delay. To save computer time, 
several thresholds are considered simultaneously, and when the 
highest has been exceeded the generation of a new process is 

started. The minimum value of from the values obtain^i 

for the different thresholds is plotted. 

The ML estimator applies a weighting function obtained from 
the desired filter response, sums the values obtained, and selects 
the time of the absolute maximum of this sum as the estimate of 
the delay. The weighting function corresponds to a filter matched 
to log |S(t) + n^l 

Figure (6) compares the simulation and analytical results for 
rms timing error for both the ML and “sliding window” counter 
estimators. Note that even though both analytical results involve 
approximations, the simulation results for both are identical with 
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FIGURE 6 VALIDATION BY SIMULATION 


from the optical detector. This is equivalent to replacing the 
“sliding window” integrator with a filter. The output of the 
weighted counter is applied to a threshold detector. This design 
comes about either by analogy with the matched filter result from 
Gaussian detection theory or by mdoifying the ML estimator by 
replacing the maximum detector with a threshold detector. 

Performance of this estimator is difficult to evaluate 
analytically because the instantaneous distribution of the counter 
output is a complex shot noise process. Figure (7) compares the 
rms estimation error of two weighted counters with the “sliding 
window” counter. One counter uses the weighting function log 
[S(t| + n]^] while the other uses S(t). Use of S(t) as a weighting 
function appears to be superior to use of log fS(t) + nj^,] and 
roughly 20% superior in rms error to the “finite window” counter. 
Thus, using a “matched” filter rather than the “sliding window” 
integrators (i.e. a filter with a square impulse response) before the 
threshold offers some performance improvement but is not as good 
as the ML estimator. 

Comparison of the previous results indicate, but definitely does 
not prove, that the use of a maximum detector rather than a 
simple threshold provides significant performance gains. Some of 
the difference between the ML estimator and the “sliding window” 
counter can be eliminated by replacing the integrator with a 
carefully designed filter but the magnitude of the improvement is 
limited. 
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the analytical results to within the expected error of the 
simulation. 

WEIGHTED COUNTER SIMULATION 

An intuitive design of a pulse arrival time estimator might use a 
weighted photoelectron counter with the weighting function 
matched to the received pulse shape. Instead of simply integrating 
the output of the optical detector over the previous pulse width, 
the weighted counter includes a weighting function, w<t) in the 
integration. Thus the output of the weighted counter is given by 

r ^ x( t) w( f — t + T/2 )dt’, where x( t ) is the sequence of impulses 
t-T 


FIGURE 7 

One way to approximate a maximum detector is to use an 
adaptive threshold or equivalently an instantaneous AGC. One 
solution to this problem is to delay the estimation of a by at least 
Tmax- This allows an estimate of the total energy in the pulse 
which is then used to normalize the received signal before 
thresholding. This system is shown in Figure 8. Figure 9 shows the 
performance of this system using both a “sliding window” and 
weighted counter (weighting function is S(t ) )with AGC. Note that 
the weighted counter with AGC is as good as the ML estimator for 
this case. The main difficulty with this approach is the realization 
of wideband ideal delay lines. 





OPTICAL 
DETECTOR I 


-ST 


MAX 


WEIGHTED 

COUNTER 




J 

^0 



FIGURE 8 ADAPTIVE THRESHOLD TECHNIQUE 


for a two state Markov process with one absorbing state. 

P( t+ A ) = Prob [no crossing before time t } ( 38 ) 

•Prob I crossing in the interval (t, t+A) 

I no crossing before time 1 1 + 

Prob { first crossing occurs before time t } 

Prob { first crossing before time t } == P( t ) ( 39 ) 

Prob I no crossing bhfore time 1 1 = 1 — P(t) (40) 

P( t+ A ) = [1-P( 1 1 ] Prob {crossing in the interval (41) 

(t,t+A) I no crossing before time t} +P (t) 


SUMMARY 

Tlio timing resolution of various sub-optimum estimators were 
compared witli the resolution of the ML estimator using both 
analvt icnl and simulation techniques. For t3q>ical signal and back- 
ground parameters the ML estimator is approximately a factor of 
2 superior in timing resolution to a simple ‘‘sliding window’* 
integrator followed by a threshold detector. It appears that most 
oi this difference is due to the use of a maximum rather than a 
tlireshold detector. This difference can be made up by using an 
adaptive threshold detector. 


Thus dividing equation (41) by A and taking the limit as A 0 
equation (42) is obtained. 

Urn P(t+A) - P(t) [1-P(t)l lim 
A-^0 A ” A-i-0 (42) 

Prob {crossing in the inverval (t, t + A ) | no crossing before 
time tf 

~ ^ - 
^ = ll-P(t)l a(t) (43) 
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FIGURE 9 ADAPTIVE THRESHOLD SYSTEM PERFORMANCE 


APPENDIX A 

Let S(t, a) be a function of time and the random vector process 
a At is desired to find P (t) where P (t) is defined to be the proba- 
bility that the first t for which equation (37) holds is less than t. 


where a( t) = lim 1 Prob { crossing in the interval 

A-*-0a (t,t+A) Ino crossing before time t| 

Thus, 

P(t) ■= 1 - (l-Pg) exp a(t')dt'] (44) 


where P(0) =* P^. 

Equation (44) is a formal solution to the first level crossing 
problem. The remaining problem is the determination of the 
transition probability a(t). 

Let S(t + A, a) - S(t, a) - AS (45) 


If S(t, a.) is assumed to be mathematically well behaved, then in 
the limit as A -► 0 , S(t, a) can have only one zero crossing in the 
interval (t, t + A). 

a(t) - lim i Prob {AS>L-X | S(t,a^) = X, S(tV, a) <L, 
A-»0 

, (46) 

0 1 t* 1 Prob {S(t, a) = X I S(t', a) <L, 


0 £ t* £ t> dX 

' Assuming S(t, a) is a suitable smooth function, AS can be 
determined from equation (47). 


AS 


dS(t,a) 


At 


dt 


i=l 


dS(t,^ Aa. 


da . 


(47) 


If S (t, ^) is not a continuous function of t, P ( t) must be updated at 
each discontinuity using the laws of conditional probabilities. 
Equation (47) shows that the distribution of A S can be found if the 
joint distribution of the is known. 

For many cases of practical interesta or S(t, a) is represented as 
a vector Markov process. When S(t, £) is a vector Markov 
process, 


S (t, a) > L ( 37 ) 

P(t) can now be determined by the same reasoning as that used 


Prob {AS>L-X |S(t, a^) = X, S(t', a) <L, 0£t*<t} ( 43 ) 

« Prob {AS>L-X |S(t, a) = X} 
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Note that the right side of equation (48) can usually be evaluated 
without difficulty. Using the Markov assumption, Prdb { s (t , ^ 
“ X I S(t, ^ <L, 0 £ T ^ t}can be determine from, 

p <x,t) 

Prob{S(t,^- x|S(t’,a)<L,0<t'<t} ■ . — (49) 

/ p. (y.t) dy 

0 • 

where Pg (y,t) is the solution to the diffusion equation for S(t, a) 
with an absorbing boundary at S(t, a) = L. 

Using the idea of an absorbing boundary, several authors ^ i ^ 
have considered the special case where S(t, a) is a process with 
independent increments since the first crossing density can then 
be determined directly by solving a Fokker-Planck equation. For 
this case the derivative of S(t, a) with respect to time exists almost 
nowhere and it can be argued intuitively 
lim Prob{S(t,a)*x|S(t*,a)<L,0<t:’<t}* ? (X, t) * 0 

X->L X^L 

X£L X£L 

Unfortunately, the zero boundary condition comes directly from 
nonphysical attributes of independent increment processes (for 
instance, an infinite number of level crossings in zero time). There 
appears to be no general technique for evaluating the boundary 
condition for Prob S (t,^) = X I S(t', a)<L, 05t’< t| at X = L for 
physically realizeable processes. Yang and Shinozuka ^ have 
suggested a way around the unknown boundary value problem by 
formulating the solution in terms of the Kolmogorov backward 
equation rather than the forward (Fokker-Planck) equation. 

Perhaps the most valuable attribute of this whole approach is 
not that it leads to exact solutions of previously intractable prob- 
lems (which it usually does not), but that it expresses the first 
crossing distribution in terms of quantities that can often be 
accurately approximated. 

APPROXIMATION FOR THE CONDITIONAL DENSITY 


The conditional density of A S can often be exactly determined, 
however general solutions for the conditional density of S(t, a) are 
not known. One useful approximation is given by equation (50). 

Prob{S(t, a)«x|s(t*, a)<L ,0£t * <t}^ Prob{ S(t, a)»x| S (t ,a)<L} (50) 


By definition, 


10 


for X > L 


Prob {S (t , a)*X | S (t , a) < L}* 


llXotJ 

f(X,t)dX 


(51) 

for X < L 


where f(S,t) is the distribution of S(t, a). 

The use of this approximation and this general approach in 
obtaining closed form expressions for the first crossing density for 
various types of random processes is discussed in reference 4 . 
The approximation makes uses of the fact that for many processes 
of practical interest specifying that the process has not crossed 
some level in the past does not significantly alter its present 
distribution. 

We now consider the optical pulse detection case where photo- 
electrons are being emitted from an optical detector according to a 
Poisson distribution with mean rate X(a , t). A “sliding window” 
counter, that is a counter which determines the number of photo- 
electrons received during the previous T seconds, is used to filter 
the optical detector output. The output of the counter is then 
applied to a threshold detector which estimates pulse arrival time. 


LetN(t) = the number of photoelectrons received in the interval 
[t — T, t] (i.e., signal and background). 

Define, 


P(t) * Prob (N(s) <L, 0<s<t}. 
Thus P(t) is determined by equation (44) where 


a(t) - lim i Prob {N(t+A)>L|N(t ’) <L, O^t'^t) 
A-K) 4 
L-1 

- Urn i ^ Prob {N(t+A)^L|N(t)-M,N(t’)<L, 
A-K) “ M-1 

•Prob {N(t) - M |N(t’) <L, 


Prob {N(t)-M|N(t)<U= 


/-..sM 

(nT) 

m: 


. exp(-nT) 


L-1,-„.K 

E (nT) _ 

— e xp(-nT) 

K=o k: 


for M <L (62) 


where n is the effective mean rate for the interval [t — T, t]. 


In the limit as A ->■ 0, j Prob {N(t+A) >_ L | N(t)=M, N(t) 
£L, 0 £ t £t } approaches zero unless M = L — 1. For M = 

L— 1, the probability of a threshold crossing due to one or more 
counts dropping out and two or more counts occuring in time 

A also goes to zero faster than L 

A" 


A^ f- Prob {N(t+A)^L|N(t)=L-l, N(t’)<L,0<t’<t}- n (53) 

where n is the instantaneous mean rate at time t. For the case 

where L — 1, a(t) can be evaluated exactly since 

Prob {N(t)-M|N(T) <l,0£T£t} = Prob {N (t)-M| N (t) 0 (54) 

. Thus for the case of background alone ( n = ), the false alarm 

probability in a time period t is exactly given by 

P^a ” l-exp(-nj^t) (55) 

For the general case 


a ( t) :fcn( t ) 



(56) 


where n(t)« X(t ,a) is the instantaneous mean rate at time t. Note 
that the density of the first crossing time can easily be expressed 


P(t) 


dP£^ 

dt 


a(t)exp[-J* a(t')dt*] 


(57) 
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APPENDIX F 


TARGET SIGNATURE ANALYSIS 

1.0 CUBE-CORNER GAIN AND ANTENNA PATTERN 

For this study, the LAGEOS type of optical retro-reflector was chosen as 
the model retro-reflector. The LAGEOS retro- reflector is a 3.8-cm diameter, 
fused silica cube-corner, with an intentional dihedral angle offset of 1.5 
+^0.5 arcsec. This specific design was chosen to optimize the return at 0.6943 
ym wavelength for ranging from the ground to a low to medium altitude satellite. 
Over the broader range of conditions encountered in this Study we find the 
cube-corner is usable but not necessarily optimum. 

The cube-corner is modeled as a receiving/transmitting aperture to which 
an incident phase wave is applied. As the wave passes through the cube-corner, 
the phase is perturbed, both by systematic effects, such as dihedral angle 
offsets, and by nonsystematic effects, such as temperature gradients. The 
amplitude of the wave is also affected. The LAGEOS cube-corners use total 
internal reflection, which cause a linearly polarized incoming wave to be 
decomposed into two, nonuniformly^^ ^ distributed, orthogonally polarized 
reflected waves, and some net energy loss is encountered. For the purposes of 
this study, both temperature gradient and polarization state transfer effects 
were ignored. 

We define the gain product of an ideal cube-corner retro-reflector as the 
product of the effective receive aperture antenna gain and the effective 
transmit aperture antenna gain. At normal incidence, the effective aperture 
is the face diameter, and the on-axis gain, with respect to isotropic radiators, 
is simply, 

G(o) = (ttD/x)‘^ (F-1) 

D = face diameter, m 
A = wavelength, m 


^ ^R. F. Chang, et al , "Farfield Diffraction Pattern for Corner Reflectors with 
complex reflection coefficients," Journal of the Optical Society of America, 
Volume 61, Number 4, April 1971, ps 431-438. 
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The far-field antenna pattern is then given by. 


G(e) = (ttD/X)^ 


2 J^(it9D/A) 
[ (TreD/A) 


12 


(F-2) 


This is the well known (J^(x)/x)^ pattern from a circular aperture vn'th uniform 
amplitude and phase distribution within the aperture. 


In order to compute the antenna pattern when the effective aperture is 

noncircular, and the phase is nonuniform, we must resort to numerical techniques. 

The starting point is the Fraunhoffer approximation to the Fresnel -Kirckhoff 
(21 

integral . ' 


E(x2,y2.L) = if 


i KL - 


2L^^2 


4 ) rr . 




E(x-| ,y-| ,o) exp|^(xi ,X£ ’’’ ^1^2^} 


— CO 


E is the complex spacially dependent "wave function," and x-| , y-j are dimensions 
in the plane containing the transmit aperture, X 2 , y£ dimensions in a 
parallel plane, spaced a distance L from the, aperture plane, which describes 
the location of the point detector, and K is the wave number, 2ir/A. The 
exponential external to the integral is simply a phase term, and may be dropped 
from this expression for our purposes. 

The power within the aperture is simply, 

00 

’^1 *o) dXydy^ (p_4 j 


When the amplitude is uniform over the aperture, as we have assumed. 


E(x^.y^ ,o) = Eq exp |if(x^,y^)|, 

within the aperture and zero elsewhere 

P = E^ A 
^A “"o ^ 

where A is the aperture area. 


(F-5) 


( 2 ) 


Klein, Miles V.; 


"OPTICS," John Wiley & Sons, New York (1970), pg 421. 
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The field from an isotropic radiator with power 1s, 




I aeI 

»y 2 »L) = "w ■ A 
▼ 47tL^ 


(F-6) 


The cube-corner receive gain, compared to an Isotropic radiator, Is 
2 

4ttA/x . Thus, the normalized field at the detector 1s, 


^ E(x«,y«,L) r 

E^(x2,y2.L) = t^(x2.y2,L) 


4ttA 

,2 


(F-7) 


C^(x2-y2-L) = % JJ exp ^|f(x,,y,) + f<x,x2 + 


If 


dx^dXg 


Next, let (|)j^ = X 2 /L and a = y 2 /L 


Then, 

^l\|(<f>x><i>y) = -y ff 

^ A 

Equation CP-8) reduces to Equation CP -2) for uniform phase tf Cx-j ,y-])=o) . 

In this form, we clearly see that E^ Is simply the two-dimensional Fourier 

•j f (xi ,yi ) 

transform of ^ ' , taken within the aperture. Since we wish to describe 

the entire far-field pattern over the region of interest, it is most economical, 
numerically, to use the fast Fourier transform (FFT) technique. This requires 
representing the function to be transformed as a set of samples. 



K(xi<l)x + yi<|) ) 


dx.|dyi 


(P-8) 


P{xi,yi,o) 

-£ Ej^^^ 6(x-| - mAx,y^ - nay) (F-9) 
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where 


^m,n ~ , mAx, nAy within the aperture. 


and 


Then 


= 0 elsewhere, 

- 00 . 

E|Nj(<l'x*<J>y) = %AxAy)^ ®^P |iK(mAX(|)j^ + nAy(|)y)j^ 

— 00 

Gf^(<l>X»<l>y) = EN(<|)j^v<(.y) • (<|)x><!>y) 


(F-IO) 


A short digression on the FFT is in order. The basic purpose is to com- 
pute the Fourier transform of a periodic process, represented by a set of 
samples. There are two significant requirements for reliable application of 
this technique. First, the sample spacing must satisfy the Nyquist criteria. 
Second, the sample set must be sufficiently long to provide the desired output 
resolution. 

In spatial terms, a one-dimensional transform, with sample spacing Ax, and 
a length Nj (samples), results in a sample spacing after transform, of A<i) = 
X/NjAX. 

The domain of the sample set, after the transform, is +(|)p(<f)p = A/2Ax). In 
this context, (f)p can be considered as the folding frequency. If we choose Ax 
sufficiently small, and Nj sufficiently large, the transformed sample set is 
a good estimate of the Fourier transform of the nonperiodic process. 

The two-dimensional transform process is accomplished by first transform- 
ing each row of a sample matrix, S(I,J), and then transforming each column 
(or conversely). Finally, we note that the sample matrix (MxM) may be smaller 
than the transform set (Nj), provided only that the transformed sample matrix 
covers the domain, in transform space, that is of interest. 

The output antenna pattern in units of receive/transmit gain product is 
obtained by processing the transformed sample matrix, 

G(lA(}>^,jA(|)y) = ^4TrAXAy/X^^ . S(I,J) • S*(I,J) (F-11) 

The next step is to establish the process for obtaining the sample matrix. 
This would be a trivial operation for an open cube-corner with no dihedral 
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angle offsets, since the output aperture is simply the input aperture folded 
successively about each axis, and a ray which intersects the face plane within 
both the entrance and exit aperture is a valid (nonzero) sample ray. However, 
for the solid fused silica cube-corner, with dihedral angle offsets, it is 
necessary to account for both refraction and phase taper. 

The technique we employed was to construct a rectangular sampling matrix 
normal to the input wave direction vector. The intercept of each ray with 
the face plane was computed, and the ray was refracted into the cube. If the 
ray was within the entrance aperture, the reflection surface contact sequence 
and exit coordinates were computed. If after reflection the ray was within 
the exit aperture, the phase was computed based on the contact sequence and 
the magnitude of the dihedral angle offset, A complex variable with unit 
amplitude and the computed phase was then stored in the sample matrix. 

After the sample matrix was completed, the two-dimensional FFT was 
performed, followed by the scaling conversion shown in Equation (F-11). 

The resultant data matrix was processed by several data presentation 
subroutines. The subroutine used for most subsequent data processing con- 
verted the data from the rectangular sample spacing into the mean gain in 
annular slots, typically 2.5 urad wide. For the final production runs, the 
annular slot gain data was stored in a file for subsequent access by other 
programs. 

The complete data file includes the annular slot mean (and variance) 
gain data for the nominal cube-corners, at incidence angles from zero to 40° 
in 5° steps, for both 0.53 urn and 1.06 ym illumination, in 2.5 yrad bins 
from 2.5 yrad (center) to 72.5 yrad. During execution of this program, 

180.000 gain points were generated, which were compressed to 540 values 
for the mean and variance of the gain product in the annular slots. This 
program, and two of the data conversion programs used to process the gain 
data, are presented in Appendix A. 

2.0 LAGEOS SIGNATURE DATA 

Conceptually, the LAGEOS type of target presents a uniform signature 
regardless of the direction of the incident radiation, resulting in a target 
which can be used for ranging purposes both from the ground and from a space- 
craft in a high altitude orbit. For this reason, the LAGEOS was chosen as a 
prototype for the passive low altitude target (LATS). 



The LAGEOS configuration, which was used in this 'study, is shown in Fig- 
ure F-1 . Retro-reflectors are installed in rings (constant latitude), and slightly 
recessed into the surface. Each retro-reflector is identified by a location 
vector in polar coordinates, summarized in Table F-1. When illuminated by an 
impinging laser pulse, with a direction vector v, each retro-reflector, whose 
position (in terms of central angle between the retro-reflector location 
vector and the complement of the pulse direction vector) is less than some 
critical angle, will contribute to the return pulse. The angle between the 
first contact point and the retro-reflector location vector is given by 
Equation (F-1 2), 

a = ARCCOS^COS(LATp) COS(LATp) + SIN(LATp) SIN(LATp) COS(LONGp-LONG^)^ (F-12) 
This is shown, schematically, in Figure F-2A. 

LAGEOS CONFIGURATION 



32 
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FIGURE F-1 



ILLUMINATION GEOMETRY 



FIGURE F-2 
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TABLE F-1:LAGE0S RETRO- RE ELECTOR LOCATIONS 

(UNIFORMLY SPACED WITHIN RINGS) 
TYPICAL, EACH HEMISPHERE 


RING 

NO. OF 
RETROS 

RING 

LATITUDE 

(DEGREES) 

LONGITUDE 
OF FIRST RETRO 
(DEGREES) 

1 

32 

4.87° 

5.625° 

2 

32 

13.25° 

0° 

3 

31 

22.98° 

5.806° 

4 

31 

31.23° 

0° 

5 I 

27 

40.96° 

0° 

6 

23 

50.69° 

0° 

7 

18 

60.42° 

0° 

8 

12 

70.15° 

0° 

9 

6 

79.88° 

0° 

10 

1 

90° 



Figure F-2B shows this geometry in the plane containing the two vectors. 

The pulse reflected from this retro-reflector must travel 2d further than if 
it were aligned with the pulse vector. Thus, for the nth retro- reflector, therj 
is a characteristic delay time, t . 

n , 

2r 

\ = “n) 

r^ = radius of the sphere (m) 
c = speed of light (2.997925 x 10® m/s) 

Now, if the laser pulse has an amplitude time envelope, g(t), g(o)=l, 
the reflected pulse has a mean amplitude envelope, S(t). 


N 

S(t) = £ g(t - T ) a 
n=l ^ 


(F-14) 
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where 


is the gain product constant for the n^^ 
retroref lector, with an incidence angle, a^. 

The spherical surface area of the annular ring defined by a, a+Aa, is 
simply 2Trrg(Ad). Thus, if the spherical surface area is divided into N such 
rings, we would find Ad = 2Trrg/N. Thus, a LAGEOS type of target with a 
"uniform" retro-reflector distribution would have one retro- reflector in each 
annular ring, resulting in a target impulse response consisting of a series 
of equally spaced impulses (At = STTr^/CN), whose amplitudes nominally decrease 
with increasing delay. 

Now, if the incident laser pulse is temporarily Gaussian, the incident 
pulse g(t), is defined by Equation (F-15). 


g(t) = exp 



(F-15) 


Then, by substitution in Equation F-14 S(t) is given by Equation (F-lb) 


N (t-T 

S(t) = E a exp 5 — 

n=l ^ r ^ ^ ” 


(F-16) 


The exp{jw^t} term is simply the laser carrier frequency term. 

The = (je^) term is a phase term which is, effectively, a random variable, 

uniformly distributed in 2 tt radians. When the return pulse is detected, the 
detected signal, X(t), is given by Equation (F-17). 


N N 

X(t) = S(t) S*(t) = £ £ a^a^ cos(e^-e^) exp 

n=l m=l 


(t-T„)^ + (t-T 


(F-17) 


Now, the a^ are random variables whose randomness is principally attributable 
to angular distribution variations. The are random variables determined by 
®n " ^'^c^n^ large, and are therefore relatively 

independant. Finally, the x^ are random variables, where the randomness is 
attributable to the randomness in the latitude and longitude of the incident 
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wave direction vector, in terms of the target coordinate system. The expected 
value of X(t) can be obtained by first taking the expectation with respect to 
0^, then with respect to a^. This results in Equation (F-18) 




(F-18) 


The term, '•s simply the mean power gain of the retro-reflector at the 

incidence angle corresponding to 

The next step in the process is to determine the relative magnitude of 
the various contributers to the randomness of X(t). The first step is to 
integrate X(t) over all time, resulting in Y, the total energy in the pulse. 


00 

Y = / X(t) dt 

— CO 



a a cos (9 -e ) exp - 
n m ■ ■ ' n m' ^ 


2T^ 


(F-19) 


Then, let. 


Then, 



(F-20) 


(F-21) 


(F-22) 


(F-23) 


F-10 



2 

Now, if we consider E{Y'} as the "signal", and 0 as the noise, due to gain 
variations, the target signature signal-to-noise power ratio (SNRq) is given by 
Equation {F-24). 



(F-24) 


Thus, we conclude that the mean, or many pulse, integrated average detected 
power SNR is dependant only on the gain variation noise, assuming variations 
in location cause negligible effect on the mean gain. Figure F-3 shows the 
dependence of SNRq to point-ahead angle magnitude at both 0.53 ym and 1.06 ym 
wavelengths. These curves were computed using the data from the cube-corner 
gain analysis discussed in Section 1.0 of this Appendix, using the program, CONVV, 
shown in Appendix B. The data shows that the effective target SNR, although dependant 
on point-ahead angle and wavelength, is on the order of 15 to 20 dB, which is 
considerably greater than the expected shot-noise signal-to-noise ratio for 
narrow pulse ranging. 

In order to evaluate the effect of the nonuniform distribution of retro- 
reflectors over the surface of the LAGEOS sphere, we resorted to a simulation. 

For each time (t), an incident wave direction vector was chosen randomly (over 
4 tt steradians), and E{X(t) x^} was computed. The variance of E{X(t) x^} was also 
computed. This process was repeated on the order of 50 times for each sample 
time, and the mean and variance of the estimates of the mean gain, due to incident 
wave direction variations were calculated. The variance of the estimate of the 
mean gain, due to direction variations, was uniformly small compared to the 
variance of the gain due to gain variations. These data were generated for a 0.2 
ns laser pulsewidth. The results, therefore, allowed us to consider, in all 
subsequent work, that the LAGEOS was well modeled as having a uniform distribution 
of retro-reflectors. For a less symmetrical target, this assumption cannot be 
made, and target signature variations with incident wave direction vector must be 
considered. 

The effect of point-ahead angle and wavelength on E{X(t)} as defined by 
Equation F-18 using the cube-corner gain data, discussed in Section 1, of this 
Appendix, was evaluated for a temporally Gaussian incident pulse of 0.2 ns width 
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TARGET SIGNATURE SIGNAL-TO-NOISE RATIO DUE TO GAIN VARIATIONS 




p 

at the e" power point. This data, and the generating program, are presented in 
Appendix B. A considerable variation in mean return pulse amplitude, shape, 
and width was observed. 

We can examine the effect of the variation in mean return pulse amplitude, shape, 
and width on a split gate range tracking loop by considering the function f{e). 


f(e) = f X(t) dt - 
e 


e 


J X(t) dt 


(F-25) 


where x(t) is given by Equation (F-17). 


With E{X(t)} defined as in Equation (F-18), this results in an expected 
value of f(e) given by Equation (F-26). 

g(e) = E|f(e)} = ^ Sp' E ^{®n} (F-26) 


where. 


erf(x) = — 



If we expand g(e) in a Taylor series, and ignore the second and higher 
ordered terms, we find the Eq, such that qIeq) = 0, is given by Equation (F-27). 


tV7 

r M 
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n=l 


erf- 

T /T^ 

i V j 
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N 

Ln=l 


[ exp ■ 

{- 


(F-27) 


Equation F-27 is evaluated in Figure F-4 as a function of point-ahead angle, 
for two laser pulsewidths and at both 0,53 and 1.Q6 ym wavelengtiis . The program 
used to compute this data is shown in Appendix C, 
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STATIC BIAS - LAGEOS VS POINT-AHEAD 
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When the incident pulse is very short (T « x.,), and x = nAx, e is within the 

IN n 0 

range of values which satisfy Equation (F-29). 

INT(£q/At) N 

^ \ = T. % (F-29) 

n=l n=l+INT(ep/Ax) 


In addition to computing e^, we also computed the two approximations, and 
found close agreement between computed using Equation (F-27), and the 
approximations, where we compared the 200 ps pulse to the very short pulse 
approximation and a 4 ns pulse to the very long pulse approximation. The 
calculated pulse centroids, for the 200 ps pulsewidth, were also shown 
on the signature curves in Appendix B. We conclude that is the cube-corner 
response is uniformly close to the "ideal," which was used to compute the 
mean gains, we can easily compensate for the static bias, since the point- 
ahead angle is well known in operation. 

Note, however, that manufacturing tolerances and thermal gradients cause 
variations in the far-field response of the individual cube-corner. These 
variations are most like to significantly affect the location and depth of deep 
nulls, which are only observed at near normal incidence angles. Thus, we 
would nominally expect to find manufacturing tolerances and thermal gradients 
having maximum influence on the leading edge of the pulse. The nominal, 
near normal incidence angle antenna pattern at 0.53 ym (Appendix A, page A-9) 
exhibits a significant reduced gain region at small point-ahead angles. This 
reduced gain is responsible for the significant rise in static bias at small 
point-ahead angles as shown in Figure F-4. This effect disappears at 
1.06 ym, since the aperture is only half the number of wavelengths in diameter. 

The effects of thermal gradients on the far field antenna pattern are, 
of course, a function of the type and direction of the gradient. The axial 
gradient, i.e., where the temperature within the corner-cube varies linearly 
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from apex to face, and is constant in planes parallel to the face is a good 
example. This type of thermal gradient produces a parabolic phase shift in 
the reflected energy, which may either enhance or combat the phase shift 
caused by dihedral angle offset. Depending on the direction and magnitude 
of the gradient, then, the far-field pattern can vary from near diffraction 
limited performance for zero phase error, to many interference rings within 
the nominal range of point-ahead angles. Thus, the extremely important 
leading edge of the return pulse is strongly influenced by local conditions 
with the target, and we could hypothesize that the static bias could easily 
very, at small point-ahead angles, between the 80 ps value calculated for 
0.53 pm and the 340 ps value calculated for 1.06 pm. Gradients of 1 to 2°C 
are sufficient to cause this magnitude of effect. This should be much less 
severe at larger point-ahead angles. Note that for ranging between two space- 
craft, small point-ahead angles are encountered only for nearly coplanar space- 
craft, and even then, only for short periods. For a HATRS/LATS link, with 
coplanar spacecraft at 22000 km and 1000 km altitudes, respectively, the point- 
ahead angle is less than 15 prad for about 36% of a viewing opportunity, roughly 
centered between acquisition (at horizon) and closest approach. This effect 
should be virtually unobservable from the ground, since the point-ahead angle 
nominally remains large over a viewing opportunity. Finally, we note that the 
potential variation in static bias is relatively unaffected by laser 
pulsewidth. 

3.0 COHERENT INTERFERENCE EFFECTS 

In the previous section, the LAGEOS signature was computed, assuming 
noncoherent (power) summing of the returns from the individual retro-reflectors. 
In this section, we address the more general situation, and determine the 
statistics of the processes as fully aspossible. Figure F-6 shows the detection 
process in block diagram form. Previously, we noted that the return pulse, S(t), 
was simply the convolution of the laser pulse with the impulse response of the 
LAGEOS. Since these processes are linear, the order is immaterial, and is shown 
in Figure F-5 to emphasize the contribution of laser pulsewidth on the overall 
process. The responses are summarized below. 

\(t) = l; s(t-T ) (F-30) 

*■ n=l " " 
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(F-31) 
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S(t) = hL(t)*hp(t) = E ^ ^ 

n=l ^ 


(F-32) 

x(t) = S(t) . S*(t) = E E 

n=l m=l ^ ^ ^ 

j2 jZ 

(F-33) 
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FIGURE F-5 


Equations (F-30) through (F-33) describe the processes through the detection 

block. Note that a„ and e„ are randotii variables. Y{t) can be described as a 

n n 

non-stationary Poisson process, whose mean rate varies in a random manner. 
Equation (F-32) is the same as Equation (F-16) with the carrier frequency term 
dropped, and Equation (F-33) is the same. as Equation (F-17). We would like 
to determine the statistics of the output of a perfect integrator, Z, on 
a pulse by pulse basis. 

When the laser pulse is narrow compared to the mean temporal spacing 

between individual impulse responses of the LAGEOS retro-reflectors, the phase 

terms drop out, and power summing is exactly correct. Then, the probability 

distribution of Z can be determined by repeated convolutions of the probability 

distribution of Zn, n = 1,2 •.*, N, where Z is the incremental increase in Z 

n 

due to the nth impulse. 
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Thus, we write, Pr{Z = K C„} = (C )VK!e 

n n _ n 


-c_ 


where c is a factor which 
n 


relates the retro- reflector gain, a^ , to detected photoelectrons. Then, 
Pr{z„ = K|= / Pr{z„ = K|cJf,(C„)dC„ 


Since we will be convolving a number of these individual functions, it is 
convenient to work in the transform domain, i.e., to compute the eharacteristic 
function, 


Cw) = f: Pr{z„ - k} 

^n K=0 ^ ^ 

(^2 (w) = / expic^fe'''^ - 1)} ^c^^n^ 
no' 


In order to proceed to the solution, we need to assume a model for the 
probability distribution of C. One simple model is the exponential dis- 
tribution, which has equal mean and root variance, which is a reasonable 
match to the cube-corner gain statistics. Thus, assuming f^iO) = 
results in Equation (F-35). 

♦z (w) = [%(«’*' - V) -l]’’ 

(This is the characteristic function of the Geometric distribution.) Then, 

Z has a characteristic function which is the product of the characteristic 
function of each term. 

N r . i-1 

♦z(w) = n - 1) -ij (F-36) 


Then, solving for the mean and variance of 2, we find 
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Thus, if the are small («1), the distribution will retain a Poisson 
type of characteristic; if they are large, the distribution will tend toward 
the exponential distribution. 

Other distribution models for f^(C) will alter these results somewhat, 
however we expect the small m^ cases to retain the Poisson characteristic 
virtually independant of the actual probability distribution of cube-corner 
gain. 

For the long laser pulse case, the results are similar. In this case, 
T»Nat thus, equation F-33 reduces to, 

o.2/t2 N N 

X(t) = e 

n=i 1^1 " 

which can be written as 
X(t) . 

Thus, X(t) is the product of a non-time-dependent random variable and the non- 
random pulse amplitude time function. 

The statistics of the random variable can be deduced readily. Let 
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A and B are then orthogonal random variables with zero-mean and equal variances. By 
invoking the Central limit theorem, we hypothesize that A and B tend to normal random 

variables, thus the sum of A^+B^ will have an exponential distribution. 

Then, the statistics of the output of the integrator (Z) can be deduced 

from, 

oo 

Pr{z=K[ = f Pr{z=K|s} f (s) ds 

•^0 


1 

m-K! 


/ 


0 





Hence, Z has a geometric distribution, with n2=m> 

For the case where the pulse is larger than the cube-corner spacing but shorter 
than the overall impulse response of the target, the preceding approximations 
cannot be relied upon. A simulation was conducted to determine the pulse shape 
variations that could be expected for the intermediate pulse width condition. 

Figure F-6 shows the mean pulse shape and six return pulses obtained in successive 
trials, where the optical phase terms in equation F-33 were chosen randomly for 
each pulse. A laser pulsewidth of 0.15 ns (FWHM) was used. 


As a first order approximation, then, we can divide the return signal into 

several segments, each of which is essentially uncorrelated with the others. 

This results in the probability distribution of Z determined by convolving the 

probability distributions of the individual segments. 

Thus, we can see that independant of the pulsewidth, the mean value of Z is, 
N 

2 ■ 23 '^n* variance of Z differs, however, for differing laser pulsewidths. 
Th,is is shown in the next three equations. 


N N , 

E "In + 23 > T<At (F-40) 
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FIGURE F-6 
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In conclusion, we see that as the pulsewidth increases from very small to 
very large, the mean value of the integrator output is unchanged, assuming equal 
energy per pulse, but the variance increases with the pulsewidth in the transi- 
tion region. Thus, assuming that the mean detected signal has an exponential 
distribution is a conservative estimate with the degree of conservatism 


decreasing as the pulsewidth increases to values larger than the maximum target 


delay. 

If the target were a planar array of retro-reflectors at near normal inci- 
dence, the coefficients (a^) would be the same, i.e., the randomness is signifi- 
cantly reduced. For large numbers of retro-reflectors, the orthogonality argument 
is still potent, and the probability distribution of Z should tend toward 
exponential. This was verified by simulation, and Figure F-7 shows the cumulative 
probability distribution for four planar arrays, containing 4, 9, 25, and 36 
retro-reflectors. All four cases represent the histogram results of 1000 trials, 
and confirm the reasonableness of the exponential distribution model . Thus, for 
all subsequent work, we elected to use the exponential distribution as the limit- 
ing case, regardless of target geometry or laser pulsewidth, when the target is 
composed of several individual retro-reflectors. 


4.0 GROUND TARGETS 

Ground targets may be single retro-reflectors or arrays of retro-reflectros. 
Both types of targets were considered during the study. The basic problem is 
obtaining enough signal return to satisfy the ranging function. Since the 
ground targets do not exist, at this time, a modest amount of latitude in 
selecting cube-corner characteristics is assumed, i.e., we can tailor the cube- 
eorner physical characteristics to optimize the return. 

The major factor in determining the fundamental limits on the achievable 
cube-corner gain is the restrictions which must be imposed on the far- field 
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CUMULATIVE PROBABILITY DISTRIBUTION FUNCTION - SIMULATION DATA 



Z/n (TYPICAL, ALL GRAPHS) 






antenna pattern, due to velocity aberration (point-ahead angle) effects. When 
viewed from the ground, the maximum point-ahead angle occurs when the host 
spacecraft passes directly overhead. Assuming a circular orbit, the spacecraft 
tangential velocity component is the orbital velocity at this point, and ranges 
from approximately 3.75 to 7.73 km/s for altitudes of 22000 km to 300 km, 
respectively. The maximum point-ahead angle (2Vj/C) is therefore between 25 
and 51.5 yrad, assuming negligible effect from Earth rotations (Vjg 0.074 
km/sec, PAg < 0.5 yrad). The minimum point-ahead angle occurs at minimum 
elevation angle, when the spacecraft orbit plane passes through the ground 
site. The variation in point-ahead angle with elevation angle under these 
conditions is. 


2Vc /r. cos Ef 

PA ^ y 1 - j ( F-43 ) 

where is the spacecraft orbital velocity 

C = 3 X 10® km/s 

r = earth radius ~ 6378.16 km 
e 

r^ - s/C radius = r^ + h^ 

E = elevation angle 

For the high altitude S/C, h^ = 22000 km, Vg = 3.75 km/sec, and the point- 
ahead angle varies only slightly with elevation angle. (from 24.4 to 25.0). For 
the low altitude spacecraft, h - 300 km, V = 7.73 km/sec, the point-ahead 

Is O 

angle varies considerably with elevation angle, as shown in Figure F-8. For a 
minimum usable elevation angle of 30°, the PA varies from 29 yrad to 51.5 yrad. 


For the low altitude S/C-ground link, when the target latitude is com- 
parable to the orbit inclination, the angular orientation of the tangential 
velocity vector may be in almost any direction, thus we require the "optimum" 
cube-corner transmit antenna pattern to be circularly symmetric, with nearly 
uniform gain in the annular ring between ^35 and 55 yrad. The gain product 
for this coverage is, therefore,' 


G = 


4(TrD/x)^ 

( 4 - 4 ) 


(F-44) 
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POINT-AHEAD ANGLE VARIATION FOR 300 KM ALTITUDE S/C - GROUND TARGET 
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FIGURE F-8 


Evaluating equation F-44 results in. 


where 

and 


Gpg = 227.676 + 20 log (D) 

D =? face diameter (m) 

A = 0.53 ym 


Table 2 summarizes the gain product for several cube-corner sizes. 


(F-45) 


FACE DIAMETER 

GAIN PRODUCT 

(cm) 

(dB) 

3.8 

199.3 

5.0 

201.7 

10.0 

207.7 


For the high altitude S/C ground link, the point-ahead angle range is much 
smaller, permitting substantially larger net gain products. Further, since 
the high altitude spacecraft orbit is near polar inclination, the range of 
tangential velocity vector orientations for most target latitudes is substan- 
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tially reduced, permitting some additional freedom in shaping the cube-corner 
transmit antenna pattern. 

The next step in this process is to determine which cube-corner character- 
istics can be used to optimize the far-field pattern. One possible alternative 

is to grind the face to provide a lens effect, i.e., a physical radius of 
curvature which provides a parabolic phase shift in the return wave front. For 
relatively modest values of phase shift, this process should provide an other- 
wi s e un di s to rte d wave f ron t . 

At normal incidence, for a circular aperture, the far-field pattern of 
the cube-corner with a parabolic phase shift is shown in Appendix G to he. 


where 







(- liKri/ili) 

r • *1 
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(F-46) 


(F-47) 


The second form is obtained by expanding the Bessel functions in series form 
and rearranging the terms. This expression is convenient for numerical evaluation. 
The gain is G(4>) = E^((j>) • E^*(<j)). 

Figure F-9 shows the gain, vs point ahead angle (({>) computed from equation 
F-47 for a 10-cm diameter cube-corner, at normal incidence, as the parabolic 
phase shift, measured at the edge of the aperture, is increased from ir/2 to 13 
ir/2 in steps of ir radians. The shaded portions of the curves denote the regions 
where the gain product is greater than 200 dB. In the nominal region of interest, 
the shading is deeper. The shift of energy from the center out to 50 yrad is 
clearly evident. The ripples in the gain curves result from the physical size of 
the aperture and exhibit a period of ^<6.5 yrad, which is the nominal period of the 
diffraction pattern for the uniform phase front case. 
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FIGURE F-9 


POINT-AHEAD ANGLE (ji RAD) 



Figure F-10 shows the gain product vs point ahead angle for a 5 cm diameter 
cube-corner. In this case, the region above 194 dB gain product is shaded. We 
note, when comparing this figure with the preceding figure (for a 10-cm diameter 
cube-corner), that analagous performance characteristics are obtained with 
approximately one-half the phase shift at the edge. 


CUBE CORNER ANTENNA PATTERN VARIATION WITH 
PARABOLIC PHASE TAPER 



Thus, based on these two figures, we conclude that the desired cube- 
corner pattern characteristics, for normal incidence, can probably be achieved 
for any reasonably sized device, and that the net gain-product in the region 
of interest increases essentially in proportion to the cube-corner face area. 

Also shown in Figures F-9 and F-lO are the gains (dashed lines on the last 
curve), for the low altitude S/C computed from equation F-44. The computed 
gain produced is within 3 dB of this gain over most of the region of interest. 

At nonnormal incidence angles, the effect of a parabolic phase shift is 
less clear. The changes in gain at nonnormal incidence angles results from two 
distinct effects. First, as the incidence angle increases, a "cos (e)" scaling 
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term begins to be effective. This scaling term, for an otherwise unperturbed 
aperture shape, results in a relatively minor beam spread. 

For a cube-corner, however, as the incidence angle increases, the effective 
aperture decreases rather significantly. The basic concept is that the path thru 
the cube-corner can be shown as a straight line, which is analagous to folding 
the cube-corner about the three reflecting planes, in the ray contact sequence. 
Then, the exit aperture is in a plane parallel to the entrance aperture, and 
spaced a distance of tan e^p. Thus, the projection of the entrance and exit 
apertures shows an offset as a function of incidence angle. When the aperture 
face has a large radius of curvature, resulting in a small parabolic phase shift, 
the net phase shift is the sum of the excess phase at the entrance and exit points. 
The resulting phase shift is also a paraboloid of revolution, whose center is 
exactly half-way between the two centers. Thus, the dashed lines in Figure F-11 , 
represent lines of constant phase. 

FACE-PLANE PROJECTION - 15° INCIDENCE ANGLE 



FIGURE F-11 

Since the paraboloid of revolution is unchanged, the defocussing effect is 
also unchanged. Thus, as the incidence angle increases, the first effect to be 
observed is simply the area decrease, and the transmit beamwidth is virtually 
unaffected. However, as the incidence angle continues to increase, the nominal 
diffraction limited beamwidth of the cube-corner increases, so that the effective 
beamwidth begins to increase, resulting in an increasing steeper rate of reduction 
of gain. 
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We find that an approximate gain law vs incidence angle is G a 10 ' , 

where 0 is the incidence angle in degrees, is a quite reasonable approxiraation 
in the region of interest, and seems to hold fairly well provided the beam is 
spoiled at normal incidence angles. 

Thus, we can now project the expected performance capability for optimized 
cube-corner retro-reflectors which are used singly or in arrays for ground 
targets. Equation (F-48) presents this projection for the low altitude S/G target 
configuration. 

Gpg - 225 + 20 log 0 - 6/2.5 (F-48) 

D = diameter (m) 
e = incidence angle, degrees 

For the HATRS/ground link, the target is assumed to consist of a planar 
array of cube-corners, in a 4 x 4 arrangement, which is used only for local 
elevation angles of 60° or above. The individual cube-corners are 6- cm diameter, 
and spoiled to maximize the return at 25 yrad point-ahead. The expected gain 
for each cube-corner is '^'200 dB at 30° incidence angle (60° elevation), yielding 
a mean array gain of ^^212 dB. At 60° elevation angle, the physical dimensions 
of the array result in a significant amount of pulse spreading. Alternatively, 
the array could be mechanically gimballed to reduce the incidence angle. Several 
techniques seem reasonable for this purpose, and could be used in conjunction 
with environmental protection devices to reduce periodic maintenance requirements. 

In any event, it is reasonable to assume a ground target gain product of 
210 dB for use by the HATRS radar, which results in comparable link margins for 
the HATRS/ground and the HATRS/LATS links. For the Shuttle or Space! ab experiment, 
however, the target viewing opportunity is rather short, and it would be pre- 
ferable to choose a ground target configuration which resulted in a single cube- 
corner response, to facilitate dual frequency (1.06 ym and 0.53 ym) ranging for 
atmospheric delay evaluation. 

Thus, for reasonably sized cube -come rs , the maximum usable mean target gain 
product is considerably less, more on the order of 192 to 196 dB at 0.53 ym, and 
186 to 190 dB at 1.06 ym wavelength, for incidence angles up to 30°. Environ- 
mental protection is desirable, although possibly not required during a Shuttle 
mission assuming favorable locations and weather conditions. 
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APPENDIX G 


FAR FIELD RADIATION PATTERN 
TRUNCATED GAUSSIAN APERTURE 


1. INTRODUCTION. This appendix presents a solution of the integral describing 
the far field radiation pattern of a finite aperture illuminated with a 
Gaussian intensity distribution. 


2. ANALYSIS. The far field radiation pattern of circularly symmetric plane 
radiation from a finite circular aperture is described by, 

a 



f{r)JQ(Kr sine) rdr]^ 


(G-1 ) 


where 


0 is the angle from the optical axis 
a is the aperture radius 
K is the wave number (2n/x) 

A(L) is a propagation constant 

is the zero order Bessel function 


and 


f(r) is the electric field distribution in the aperture 


For small angles, and neglecting unimportant constants, this reduces to. 


a 

I-| = if f(r)jQ(Ker) rdr]' 


( 6 - 2 ) 


A Gaussian intensity distribution in the aperture is described by, 


f(r) = exp {-r^/R^} 


(6-3) 


where R is the radius of the e" power point of the intensity distribution. 

Substituting equation (G-3) into equation (6-2) results in, 
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Next, we substitute t = resulting in, 

i \ 

= I y* e ^ J^j(Ke/t)dt 
D 

Now, the Laplace transform of J (2/bt) 
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Then, let. 
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We find, therefore. 
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Next, we find, 
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This is the result obtained by solving equation (G-2) with e=0. In order 
to establish the limit of as a «>, we use the identity, 
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z (2a/K0R^)^ J.(Kea)e ^ = e - z (||^) J„(K6a)e ^ 

i=l ’ n=0 " 


Whence, 


1 


2 12 

Lim Ii = ^ R e ^ 2 ' 

a oo 


(G-15) 


This is the expected Gaussian intensity far field pattern for a nontrun- 
cated Gaussian plane radiation. 

Last, we note that. 


? 2 Ji(K0a) 


Lim I-i = a ( 


K0a 


) 


(G-16) 


Thus, as R the aperture illumination becomes uniform, and the 

p 

expected (Ji(x)/x) intensity pattern is obtained. Note that Equation (G-12) 
also applies when a defocus condition is postulated. In this application 
we simply replace with the complex quantity. 


Y = 


_ 1 + j 


..x2 
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